-  AD-A236  743 

REPORT  "  WMIlilir . . 


Ohm  No.  o7o*4ita 


i 


catKOaiar 

Oa 


_  .  inpiiH  iiHi«  ■»  igr  maun  mn  aural,  io  MatiHnatax 'waMuanltl^ 
HO*.AHiima«.<A  22202-AMl.  ana  taut*  OMt««a«Mai<aqa<w<  an*  tuaqpt.  a; 


w  >aai— nn  ingfmnuw^  laarcnnn  aniana  aata  lOMwai. 

”**  wpawa  o>  p^aglar  aiaaci  at  nwa 


Ovacnraa  asr  Miona 

A 


ana  RaaaaOk  Utl  ii 


2.  report  date 

3.  REPORT  TYPE 

And  OATES  covered 

June  1,  1991 

Final 

8/1/84  -  3/31/89 

I.  AaiNCY  USI  OfN.r  (Lm««  M«iurj 


4.  TITU  ANO  SUBTTTU 

Physics  and  Chemistry  of  Small  Scale  Structures  for 
Modern  Devices. 


«.  AUTMOmS) 

T.  C.  McGill 


S.  PUNOWO  NUMSIRS 


Contract  N00014-84-K- 
0501 


7.  PEMOIUMINO  ORGANIZATION  NAMI(S)  ANO  AOORfSStiS) 

Dept,  of  Applied  Physics  128-95 
California  Institute  of  Technology 
Pasadena,  CA  91125 


t.  PERfORMING  ORGANIZATION 
RfRORT  NUMMR 


s )  I 


V..  ^ 


N/A 


9.  SPONSORING/ MONITORING  AGENCY  NAMC(S)  ANO  , 


Office  of  Naval  Research 


IJ 


10.  SPONSORING/  MONITORING 
AGINCr  REPORT  NUMBER 


11.  suppumentart  notes 

N/A 


1ZA.  OISTRlBUTION/AVABJIBaiTY  STATEMENT 


j  izbl  distribution  COOE 


Available  to  public 


13.  ABSTRACT  (MuamomiCOworOtt 

Under  the  support  of  this  contract  a  number  of  major  contributions  were  made  to 
the  development  of  nanostructures  for  microelectronics.  Fundamental  studies  of 
growth  have  led  to  a  better  understanding  of  the  role  of  strain  in  epitaxial 
layers  grown  on  semiconductor  substrates.  The  possibilities  for  growing  super¬ 
conductors  on  semiconductors  and  vice-versa  have  taken  a  major  step  forward  with 
the  successful  growth  of  vanadium  III  silicon  on  silicon.  The  basic  theory  of 
tunnel  devices  including  the  role  of  phonons  and  self-consistency  hive  been 
explored.  New  device  structures  based  on  the  combination  of  III-V  and  II-VI 
semiconductors  have  also  been  examined. 


14.  SUBlEa  TERMS 

Growth  of  epitaxial  layers,  MBE,  MOCVD,  high-speed  tunnel 
devices,  electronic  wave  function,  gallium  arsenlde/galllum 
aluminum  arsenide,  vanadium  III  silicon. 

IS.  NUMBER  OP  PAGES 

55 

18.  PRKE  COOE 

17.  SECURITY  CLASSDICATION 
OP  REPORT 

Unclassified 

IB.  SECURITY  CLASSifICATION 
OP  THIS  PAGE 

•Unclassified 

IB.  security  classipication 
OP  ABSTRAa 

Unclassified 

ZO.  UMITATION  OP  ABSTRACT 

SAR 

NSN  734041 -JEO-SSOO  SUnojrd  POfiw  (R4v.  2-89) 


a*  AMD  tra  £3*.'« 


THE  OFFICE  OF  NAVAL  RESEARCH 
CONTRACT  NO.  N00014-84-K-0501 


PHYSICS  AND  CHEMISTRY  OF  SMALL  SCALE 
STRUCTURES  FOR  MODERN  DEVICES 


by 

T.C.  McGILL 

THOMAS  J.  WATSON  SR  LABORATORY  OF  APPLIED  PHYSICS 
CALIFORNIA  INSTITUTE  OF  TECHNOLOGY 
PASADENA,  CALIFORNIA  91125 
(818)  356-4849 


91-01894 

llllllllillllllllllllllillli 


ft  6  11  172 


1 


INTRODUCTION 

The  development  of  new  materials  preparation  techniques  including 
molecular  beam  epitaxy  and  metal-organic  chemical  vapor  deposition  have 
made  possible  the  design  of  materials  at  the  atomic  scale.  These 
technological  developments  are  making  possible  a  whole  new  set  of 
electronic  concepts  that  are  likely  to  be  the  key  technologies  for  the  1990’s 
and  beyond.  Research  supported  under  this  contract  was  aimed  at 
addressing  some  of  the  key  questions  in  this  field. 

One  of  the  major  applications  of  these  new  device  preparation 
technologies  is  the  production  of  electronic  devices  that  are  likely  to  be 
replacements  for  the  standard  transistors  now  employed  in  the  period  well 
beyond  the  turn  of  the  century.  One  version  of  these  is  resonant  tunneling 
devices  in  which  electron  waves  incident  on  a  double  barrier  structure  are 
resonantly  transmitted  through  a  central  quantum  well.  This  phenomena 
leads  to  a  negative  resistance  and  is  capable  of  producing  bi-state  electronic 
devices.  Part  of  the  research  program  here  has  been  aimed  at  achieving  that 
subject. 

A  second  major  issue  that  is  addressed  in  this  research  program  is  the 
possibility  for  growing  various  materials  combinations  in  molecular  beam 
epitaxy  or  metal-organic  chemical  vapor  deposition.  Notably  in  this  work 
we  attack  two  major  problems.  One  of  them  was  a  stu^  of  the  role  of  strain 
in  determining  the  quali^  of  epitaxial  layers  grown  in  standard 
semiconductor  systems.  Historically,  one  of  the  major  issues  has  been  the 
role  of  strain  in  determining  the  compositional  clustering  of  atoms  in  thin 
films  grown  by  MBE  or  MOCVD.  Other  areas  of  intense  interest  in  growth 
are  the  possibilities  of  growing  heterostructures  involving  materials  with 
very  dissimilar  electronic  properties.  One  of  the  major  efforts  of  this 
program  was  to  attempt  to  combine  superconductors  and  semiconductors. 

Finally  ^  one  of  the  major  applications  for  tunnel  structures  is  in  high 
speed.  Hence  one  of  the  principal  activities  of  any  research  program  in  this 
area  is  to  attempt  to  make  measurements  of  the  high  speed  transient  behavior 
of  tunnel  devices. 


The  major  contributions  in  each  one  of  these  areas  is  outlined  in  the 
following  sections.  For  completeness,  copies  of  the  publications  are 
included  in  the  appendices.  One  of  the  major  outputs  of  any  research 
program  is  the  training  of  graduate  students.  A  list  of  the  graduate  students 
trained  into  this  program  is  included  in  a  section  just  prior  m  the  major 
concluding  section. 

MAJOR  CONTRIBUTIONS 
Tunnel  Devices 

Three  major  contributions  were  made  under  this  program  in  the 
area  of  tunnel  devices.  We  made  the  first  self-consistent  calculation  of  the 
electronic  wave  functions  in  gallium  arsenide/gallium  aluminum  arsenide 
interfaces  with  relevance  to  tunnel  devices.  The  major  issue  addressed  in 
this  effort  was  to  obtain  die  proper  wave  functions  at  the  interface  for 
developing  a  better  understanding  of  the  potential  profile  and  the  character 
of  the  wave  function  just  prior  to  tunneling  in  the  turmeling  device.  None  of 
the  previous  workers  had  actually  carried  out  self-consistent  calculations. 

We  find  that  there  are  major  modifications  in  the  electronic  wave  function 
and  in  the  potential  due  to  self-consistency.  The  results  of  this  research  are 
presented  in  the  paper  included  in  Appendix  1. 

One  of  the  major  issues  in  the  theory  of  tuimel  devices  is  proper 
prediction  of  the  peak-to-valley  ratio.  Much  of  the  theory  to  date  predicts 
peak  currents  correctly  but  has  great  difficulty  in  predicting  valley  currents. 
Under  the  support  of  this  prograin  we  carried  out  a  study  of  the  role  of 
inelastic  phonon  scattering  in  determining  valley  currents.  Basically, 
phonon  scattering  processes  can  effectively  broaden  the  resonance  in  the 
tunnel  stracture  and  lead  to  an  increase  in  the  valley  current  Theoretical 
studies  based  on  the  inclusion  of  optical  phonons  as  weU  as  acoustic 
phonons  in  this  process  have  indicated  that  the  phonon  scattering  mechanism 
can  change  the  valley  current  by  literally  many  orders  of  magnitude.  The 


3 


details  ofthis  research  are  contained  in  the  publication  included  as  Appendix 

2. 


Another  major  issue  in  the  development  of  tunnel  devices  is  to  come 
up  with  new  materials  combinations  which  can  lead  to  better  tunnel  device 
performance.  One  possibility  explored  under  the  support  of  this  program  is 
the  combination  of  ni-V/II-VI  heterojunctions.  Notably  systems  such  as 
gallium  arsenide/zinc  selenide  and  indium  arsenide/zinc  telluiide  are 
attractive  for  application  in  tunnel  devices.  In  these  systems  the  wide 
bandgap  E-VI  semiconductor  zinc  telluride  or  zinc  selenide  acts  as  a  barrier 
layer  replacing  the  IH-V  semiconductor  that  has  been  historically  used  as  the 
barrier.  The  results  of  this  theoretical  study  are  contained  in  detail  in  the 
paper  included  as  Appendix  3. 

One  of  the  major  issues  in  the  determination  of  the  performance  of 
tunnel  devices  and  many  other  heterojunction-based  electronic  device 
structures  is  the  value  of  the  band  offsets.  Under  the  support  of  this  program 
we  carried  out  some  of  the  first  measurements  in  our  new  MBE  esca  system 
of  the  band  offsets  for  the  extremely  well-studied  gallium  arsenide/gaUium 
aluminum  arsenide  system.  Because  of  the  unique  features  of  having  an 
MBE  connected  to  an  esca  facility  we  were  able  to  study  the  commutivity  of 
the  band  offset  for  the  (100)  structure  which  is  the  one  fliat  is  predominantly 
used  in  noiaking  devices.  The  results  of  this  effort  are  included  in  the 
publication  in  Appendix  4. 

Another  area  of  intense  importance  in  the  develq>ment  of  small 
device  structures  is  the  role  of  fluctuations.  The  atomistic  nature  of  all  the 
materials  used  in  making  these  devices  plus  the  importance  of  doping  leads 
to  an  inherent  role  of  fluctuations  in  all  electronic  devices  of  this  type.  One 
of  the  problems  being  treated  under  the  support  of  this  research  activity  has 
been  the  role  of  fluctuations.  Notably  we  were  interested  in  treating  the 
percolation  problem  under  finite  size  considerations.  The  results  of  this 
research  effort  are  included  in  Appendix  5. 

Yet  another  of  the  other  major  areas  of  interest  in  tunnel  devices  is 
their  high-speed  performance.  Tunneling  times  are  thought  to  be  measured 
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in  the  picosecond  range  for  these  device  structures.  Under  the  support  of 
this  program  we  carried  out  some  of  the  first  measurements  of  the  tunneling 
times  for  the  high-speed  tunnel  devices.  Using  a  novel  photoluminescence 
technique  which  measures  the  tunneling  of  electrons  and  holes  out  of  die 
central  well,  we  were  able  to  explore  the  transient  behavior  of  tunnel 
devices.  The  major  results  of  this  study  are  included  in  a  publication  which 
has  been  attached  as  Appendix  6. 

Growth 

Under  the  support  provided  by  this  contract  we  carried  out  two  major 
studies  of  growth.  First,  we  carried  out  a  theoretical  study  of  the  role  of 
strain  in  producing  phase  separation  in  epitaxial  layers.  In  this  pioneering 
effort  on  our  part  which  has  now  grown  into  a  full  blown  activity  by  a 
number  of  key  investigators,  we  examined  the  possibilities  for  strain  induced 
phase  separation  in  alloy  layers  grown  by  the  modem  growth  techniques. 
While  we  found  that  such  a  driving  mechanism  would  not  be  important  in 
the  gallium  arsenide/gaUium  aluminum  arsenide  system,  for  other  systems  it 
could  lead  to  substantial  phase  separation.  The  details  of  this  theoretical 
work  are  included  in  Appendix  7,  which  includes  the  published  paper. 

One  of  the  major  areas  of  intense  interest  is  die  possibility  for  growing 
superconducting  materials  on  semiconductors  and  then  growing 
semiconductor  materials  back.  Successful  development  of  this  materials 
technology  could  lead  to  a  whole  host  of  new  device  possibilities.  For 
example,  highly  controlled  Josephson  junctions  could  be  grown  in  this 
method.  The  possibility  of  exploring  semiconductor-enhanced  transition 
temperatures  in  superconductors  would  also  be  explorable  with  this  method. 
The  incorporation  of  superconductors  in  a  natural  epitaxial  way  with 
semiconductors  could  lead  to  a  whole  new  set  of  device  concepts  combining 
the  importance  of  semiconductor  electronic  properties  with  some  of  the 
novel  aspects  of  superconductors.  One  of  the  most  attractive  possibilities  for 
this  is  the  growth  of  vanadium  III  silicon  on  silicon.  Theoretical  studies 
previously  supported  by  the  Office  of  Naval  Research  have  indicated  that 
this  particular  system  has  a  good  lattice  match  between  vanadium  m  silicon 
and  the  silicon  substrates.  Hence  one  might  hope  to  be  able  to  grow  this 
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system  in  an  epitaxial  fashion.  However,  there  are  some  inherent 
difficulties.  Notably  the  desire  to  grow  a  metal-rich,  that  is  vanadium-rich, 
phase  on  top  of  a  semiconductor  or  vice-versa  in  a  very  thin  layer  suggests 
some  difficulties. 

Under  research  provided  by  this  contract  we  were  able  to  demonstrate 
the  growth  of  vanadium  m  silicon  on  silicon  substrates.  During  this  study 
we  were  able  to  obtain  thin  layers  of  vanadium  m  silicon  on  silicon 
substrates  for  superconducting  transition  temperatures  close  to  those  of  bulk 
vanadium  BI  silicon.  The  various  parameters  governing  the  growth  of  these 
superconducting  layers  on  semiconductors  were  explored.  While  the 
materials  to  date  have  not  been  single  crystal  but  have  consisted  of  relatively 
small  grain  polycrystal  in  vandium  lU  silicon,  the  success  of  growing  a 
superconducting  metal-rich  on  top  of  a  semiconductor  is  a  major  step  in  the 
right  direction.  This  set  of  experimental  results  is  included  in  the  papers  in 
Appendices  8  and  9. 

GRADUATE  STUDENTS  TRAINED  INTO  THIS  RESEARCH 
PROGRAM 

The  support  provided  by  this  program  made  substantial  contributions 
to  the  training  of  five  Ph.D.  students.  Their  names  are  Alice  Bonnefoi, 
Wesley  Boudeville,  Matthew  Johnson,  Yasantha  Rajakarunanayake  and 
George  Wu.  Each  of  these  students  is  continuing  to  make  contributions  in 
their  own  scientific  environment 

CONCLUSIONS 

Under  the  support  of  this  contract  a  number  of  major  contributions 
were  made  to  the  development  of  nanostructures  for  microelectronics. 
Fundamental  studies  of  growth  have  led  to  a  better  understanding  of  the  role 
of  strain  in  epitaxial  layers  grown  on  semiconductor  substrates.  The 
possibilities  for  growing  superconductors  on  semiconductors  and  vice-versa 
have  taken  a  major  step  forward  with  the  successful  growth  of  vanadium  m 
silicon  on  silicon.  The  basic  theory  of  tunnel  devices  including  the  role  of 
phonons  and  self-consistency  has  been  explored.  New  device  structures 
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based  on  the  combination  of  m-V  and  II-VI  semiconductors  have  also  been 
examined. 

Research  provided  under  the  support  of  this  contract  has  made  a 
substantial  contribution  to  the  development  of  microelectronics  into  die  turn 
of  the  next  century. 


APPENDIX  1 


Self-consistent  solutions  of  electronic  wave  functions  at  GaAs-Alj^Ga^.^As 
interfaces 

Y.  Rajakarunanayake  and  T.  C.  McGill 

California  Instivute  of  Technology,  Pasadena,  California  91125 

(Received  23  March  1987;  accepted  23  April  1987) 

We  report  the  first  study  of  the  self-consistent  electronic  wave  functions  at  a  heavily  doped 
GaAs-AJ^Ga,_^As  interface.  The  doping  daisities  in  the  GaAs  are  between  10*’  and  10*® 
cm“^.  These  doping  densities  are  characteristic  of  tunnel  structure  devices.  The  Schrodinger 
equation  is  solved  self-consistently  in  the  Hartree  approximation.  Some  of  the  more  important 
results  of  our  calculations  include  the  fact  that  there  is  only  a  single  subband  level  in  the  well  for  a 
wide  range  of  biases.  This  level  is  also  fairly  loosely  bound  (bound  by  a  few  meV)  in  spite  of  the 
fact  that  the  notch  at  the  interface  is  on  the  order  of  a  100  meV.  In  accumulation  layers,  the 
potential  at  the  interface  is  somewhat  similar  to  the  non-self-consistent  one.  However,  in 
depletion  layers  the  self-consistent  potential  can  be  substantially  different  horn  the  one  obtained 
in  the  Thomas-Fermi  approximation. 


I.  INTRODUCTION 


In  this  paper  we  have  studied  the  self-consistent  potential, 
and  the  associated  electronic  wave  functions  that  exist  at 
heavily  doped  GaAs-Al^Ga,  As  interfaces.  These  inter¬ 
faces  occur  very  commonly  in  single-barrier,  and  double¬ 
barrier  heterostructures  made  of  GaAs  and 
AljGa,  _;,As.'~^  For  example,  let  us  consider  a  single  bar¬ 
rier  structure  which  consists  of  GaAs  electrodes  and  a 
AJ,  Ga,  _ As  barrier.  To  describe  the  tunneling  of  electrons 
through  this  barrier,  we  have  to  study  the  response  of  GaAs 
electrons  to  an  applied  voltage  bias.  When  the  applied-elec¬ 
tric  field  is  perpendicular  to  the  interface,  accumulation  and 
depletion  layers  are  formed  at  the  interfaces  of  the  barrier  on 
either  side.  With  the  self-consistent  theory  developed  here, 
we  are  able  to  account  for  band  bending  and  the  details  of  the 
screening  charge  near  the  interfaces.  These  self-consistent 
solutions  given  here  are  of  major  interest  because  band  bend¬ 
ing  plays  a  major  role  in  determining  the  tuimeling  anrent  . 
and  other  device  properties.’^*’ 

The  first  part  of  this  paper  discusses  the  details  of  the 
Hartree  model,  and  the  important  elements  of  self-conastent 
calculations.  In  the  next  part  we  discuss  the  results  obtained 
for  accumulation  layers,  and  depletion  layers  separately. 
Two  representative  accumulation  layers,  and  two  depletion 
layers  are  studied  in  detail  to  explain  the  general  features  of 
the  solutions.  For  each  of  these  cases,  the  corresponding 
wave  functions,  carrier  densities,  and  the  potentials  are  in¬ 
cluded  in  the  figures.  One  of  the  aims  of  this  paper  is  to 
compare  the  merits  of  the  Thomas-Fermi  theory  versus  the 
ielf-consistent  theory.  Thus  the  corresponding  Thomas-  ■ 
-,,j^crmi  results  are  also  discussed  along  with  those  of  the  setf- 
'  consistent  theory.  The  last  section  of  this  paper  constitutes  a 
brief  summary  of  the  major  results.  ???  * 

II.  DESCRIPTION  OF  THE  MODEL  ' 

In  our  investigations,  we  have  used  the  self-consistent 
Hartree  model’”’  to  treat  the  accumulation  and  depletion 
layers  that  are  formed  m  GaAs-Al,Ga,  _j,As  heterostruc¬ 


tures  when  a  potential  bias  is  applied.  Right  firan  the  begin¬ 
ning,  we  can  isolate  the  two  interfaces  of  the  barrier  and 
study  each  side  separately,  since  the  coupling  between  the 
two  sides  is  exponentially  small  as  far  as  the  electrostatics  of 
the  problem  is  concerned.  Here  the  conduction  band  elec¬ 
trons  of  the  heavily  doped  GaAs  at  the  F  point  are  allowed  to 
rearrange  themselves  to  provide  the  screening  of  the  external 
electric  field.  The  motion  of  these  electrons  are  assumed  to 
be  described  by  a  simple  effective  mass  Schrodinger  Hamil¬ 
tonian.  Smee  the  F  point  of  GaAs  is  isotropic,  we  can  char¬ 
acterize  the  effective  mass  m*  by  a  simple  scalar  quantity. 

In  the  Hartree  approximation,  each  electron  moves  in  the 
average  electrostatiepotential  produced  by  all  the  other  con¬ 
duction  electrons,  and  the  uniform  positive  background  due 
to  the  ionized  donors.  However,  the  forces  acting  on  the 
electrons  can  be  separated  into  two  categories.  The  first  cate¬ 
gory  is  the  coulombic  electrostatic  forces  that  are  accounted 
correctly  by  the  Hartree  term.  The  other  category  is  the  ex¬ 
change  and  correlation  forces  that  arise  due  to  the  statistics 
of  spin  \  Fermions,  hi  principle  these  exchange  terms  have  to 
be  included  in  the  Hamiltonian  to  correctly  describe  the  mo¬ 
tion  of  the  electrons.  However  inclusion  of  these  terms  usual¬ 
ly  constitutes  a  many-body  problem,  and  cannot  be  ade¬ 
quately  described  by  an  mdependent  electron  picture.  The 
strength  of  the  exchange-correlation  force  compared  to  the 
electrostatic  force  depends  on  the  parameter  of  the  medi¬ 
um.  Here  r,  is  the  interelectronic  spacing  ea^iressed  in  the 
units  of  the  effective  Bohr  radius.  Thus  for  heavfly  doped 
materials  with  high  dielectric  constants  and  low  effective 
masses,  we  can  argue  that  the  exchange  terms  are  subdomin- 
ant.  Fortunately,  heavily  doped  GaAs  has  a  small  effective 
mass  and  a  large  cBelectric  constant,  and  thus  we  can  be 
contented  with  only  including  the  Hartree  term  of  the  poten¬ 
tial  in  the  Hamiltoinan.  The  various  approximate  methods 
available  for  estimating  these  exchange  and  correlation  po¬ 
tentials  as  functionals  of  the  local-electron  density,  that  go 
under  the  name  of ‘'density  functional  methods’*  will  not  be 
discussed  here.*’  We  have  also  ignored  the  electrostatic  im¬ 
age  forces  that  arise  near  the  interface  due  to  the  difference  in 
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the  dielectric  constants  of  the  two  materials. 

In  our  model  the  GaAs-Al^Ga,  _;,As  interface  is  as¬ 
sumed  to  be  sharp.  This  enables  us  to  treat  the  potential  seen 
by  the  conduction  band  electrons  of  GaAs  near  the  interface 
as  a  simple  potential  step  with  the  appropriate  magnitude  of 
the  conduction-band  offset.  The  band  offsets  that  occur  in 
our  structures  are  usually  about  200  meV.  This  corresponds 
to  an  AJ^Ga,  As  layer  with  a  value  of  x=:0.3.  We  have 
also  limited  our  attention  to  degenerately  doped  n-type 
GaAs  electrons,  at  zero  temperature.  However  a  similar 
analysis  can  be  performed  on  p-type  semiconductors,  and 
other  direct  gap  materials,  with  only  minor  modifications  to 
account  for  the  finite  temperatures.  The  effective  mass  of 
GaAs  used  is  m*  =:  0.065m, ,  where  m,  is  the  bare  electronic 
mass.  The  effective  mass  of  Al^Ga,_^As  used  is 
m*=:0.092m,.  The  background  material  (fiUed-valence 
electrons  and  ion  cores)  is  assumed  to  be  passive  e.xcept  to 
provide  the  uniform  dielectric  constant  of  the  medium.  In 
GaAs  f  =;  13.  Ifo,  where  is  the  permittivity  of  the  vacuum. 
In  Al;,Ga,  _,Ase=:12.2£o. 

The  potential  V  is  chosen  to  be  zero  far  away  from  the 
interface,  on  the  GaAs  side.  When  the  GaAs  is  attached  to 
the  negative  electrode,  electrons  see  an  attractive  potential, 
near  the  barrier.  Thus  an  accumulation  layer  forms,  and  the 
potential  is  lowered  below  zero  at  the  interface.  However, 
when  the  GaAs  is  attached  to  the  positive  electrode,  then  a 
repulsive  potential  is  formed  at  the  barrier.  In  this  case  a 
depletion  layer  is  formed  and  the  potential  is  raised  above 
zero  near  the  barrier.  We  have  chosen  the  z  axis  to  be  perpen¬ 
dicular  to  the  interface.  Thus,  taking  into  account  that  there 
is  no  variation  of  the  potential  in  the  x~y  plane,  one  can 
readily  reduce  this  problem  into  a  one-dimensional  one. 
However,  in  matching  the  boundary  conditions  at  the  inter¬ 
face,  the  parallel  component  of  the  wave  vector  of  each  elec¬ 
tron  k  II  is  assumed  to  be  conserved  across  the  interface. 
When  solving  the  self-consistent  problem  in  the  accumula¬ 
tion  region,  one  also  has  to  consider  the  bound  states  that  can 
occur  in  the  potential  welL 

In  a  typical  self-consistent  iteration  scheme,  we  start  with 
a  guess  for  the  form  of  the  potential,  and  solve  the  one-di¬ 
mensional  Schrodinger  equation.  From  these  solutions,  we 
then  compute  the  charge  density  due  to  all  the  other  conduc¬ 
tion  electrons,  and  the  positive  background.  Using  the 
charge  density  computed  in  the  last  step,  we  next  solve  the 
Poisson  equation  to  obtain  anew  potential.  The  details  of  the 
various  schemes  used  to  achieve  the  convergence  will  not^ 
discussed  here.'®  A  simple-minded  scheme  that  works  with 
considerable  success  is  to  modify  the  input  potential  by  feed¬ 
ing  back  a  fraction  of  the  difference  between  the  output  and 
the  input  potentials  of  the  last  iteration.*  We  have  found  that 
an  exponentially  decaying  potential  with  the  appropriate 
Thomas-Fermi  decay  length  works  as  a  very  good  initial 
guess  for  starting  the  self-consistent  iteration  scheme, 

In  this  section,  we  illustrate  the  solution  of  the  Schid- 
dinger  equation  in  the  effective  mass  approximation,  for  a 
guess  Hartree  potential  Kfz).  With  the  above  assumptions, 
we  can  write  the  wave  functions  for  this  problem  as 
The  resulting  solutions,  rf>{2)  are  the  set  of  en¬ 
velope  wave  functions,  in  the  z  direction.  The  actual  wave 


functions  are  a  produa  of  these  envelope  functions,  and  the 
appropriate  microscopic  Bloch  functions.  We  have  mea¬ 
sured  the  energy  of  the  electron  also  with  respect  to  the 
zero  of  the  potential  If  we  introduce  q  =  k^  to  simplify  the 
formulas,  then  we  can  write 

2m* 

where 


2m*  2m*  ' 

Substituting  these  relations  into  the  three-dimensional 
Schrodinger  equation  gives  us  the  familiar  one-dimensional 
Schrodinger  equation 

When  we  solve  this  equation  for  positive  energy  by  the 
above-described  method,  we  consider  incident  waves  im¬ 
pinging  on  the  barrier  from  infinity,  that  reflect  back  from 
the  barrier  with  only  a  negligible  transmitted  component. 
The  boundary  conditions  are  such  that  the  exponentially  in¬ 
creasing  solution  in  the  barrier  is  discarded.  At  the  interface, 
we  have  also  matched  the  boundary  conditions  that  the  cur¬ 
rent  and  the  amplitude  of  each  wave  function  be  the  same  on 
both  sides  of  the  interfree.  The  incident  and  the  reflected 
waves  interact  to  give  solutions  of  the  form  which  have  the 
asymptotic  forms,  xf>^i‘*~sxrr{qz  —  'q),  far  from  the  bar¬ 
rier.  The  normalization  of  these  continuum  states  have  to  be 
done  with  considerable  care.  Thus  we  introduce  a  macro¬ 
scopic  normalization  length  L,  and  a  cross-sectional  area  A 
for  the  bulk  of  GaAs  in  the  tunnel  structure,  and  let  £ — oo  at 
a  later  point  of  the  calculation.  Then  we  choose 


^,(r)~V27Z"  sinfgr— T/) 

as  z—  00.  Each  state  we  obtain  is  the  only  allowed 
linear  combination  of  -r  ?  and  —  q  states.  So  we  have  to  sum 
only  over  one-half  of  the  Fermi  sphere  with  ^  >  0  to  obtain 
the  charge  distribution.  When  there  are  bound  states,  we  also 
normalize  their  wave  functions  (z)  by  the  condition 


(z)fl!r  =  l 


for  each  bound  state.  Defining 

(z)  =  Vl/£  (z)  ...  .. 

we  can  write  the  contribndon  to  the  charge  density  from  the 
mobile  states  as  -  : -  -  r 

>n»bU)  =  “7^  f  -q^)<p*{z)<p^iz)dq  jh; 
(2ir)*a. 

after  dividing  out  by  the  factor  A  (cross-sectional  area).  In-  ; 
troducing  u  =  q/kp,  weean  write  this  expression  as 

Pmob  jf  ■ 


Thus  in  the  large  z  limit  since  the  average  value 
*  {z)q>^  (z) )  =  1,  we  immediately  recover  the  negative  of 
the  background  doping  density  .  giving  a  net  charge  den¬ 

sity  of  zero  at  large  distances  away  from  frie  barrier.  Here 


i 

I 

I 
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for  uniform  doping.  Similarly  for  bound  states, 

^bou  (2)  =  i  (1  -  (^) . 

where  n  is  the  number  of  bound  states  in  the  well.  Here£,  is 
the  energy  of  the  rth  bound  state,  measured  with  respect  to 
the  zero  of  the  potential  Thus  the  total  charge  on  one  side 
per  unit  area  is 

2  =  f  [PboaC^l  +P„ob(2) 


Jo 

Notice  that  it  is  proper  to  let  Z,  —  00  here  because  at  this  point 
of  the  calculation  all  our  formulas  are  expressed  in  terms  of 
quantities  independent  of  L.  This  excess  charge  Q  per  unit 
area  is  negative  in  accumulation  layers,  and  positive  in  deple¬ 
tion  regions,  and  is  in  fact  a  boundary  condition  of  the  prob¬ 
lem.  Now  in  order  to  find  the  potential  the  total  charge 
density  (z)  is  plugged  into  the  Poisson  equation  and  inte¬ 
grated, 

d-  e 

—  v(z)=-p^(2). 

dr  £ 

Using  v(z)  is  different  from  ^{z)  if  a  self-consistent  solution 
is  not  found.  Thus  we  iterate  this  scheme  with  different 
guesses  V{z),  until  the  self-consistency  condition 
y(z)  s=v(z)  is  satisfied  to  the  desired  numerica]  accuracy. 

III.  RESULTS 
A.  Accumulation  regions 

One  important  feature  that  we  discovered  in  our  calcula¬ 
tions  was  that  there  is  only  a  single-bound  state  in  the  GaAs- 
Al;,  Ga,  _ As  accumulation  well  for  a  wide  range  of  applied 
biases.  This  observation  was  confirmed  for  doping  densities 
ranging  from  10‘^  to  10”  cm“^,  when  the  applied  potential 
drop  in  the  layer  was  between  50-100  meV.  We  have  also 
found  that  this  bound  state  is  very  loosely  bound.  The  typical 
binding  energies  that  we  observed  usually  amounted  to  <  20 
meV,  and  were  about  5%-l5%  of  the  well  depth.  In  retro¬ 
spect,  the  evidence  that  there  is  only  a  single-bound  state 
makes  the  application  of  the  semiclassical  theory  to  this 
problem  invalid.  TheThomas-Fermi  semiclassical  theory  of 
screening  is  valid  only  in  a  regime  where  there  are  a  large 
number  of  bound  states  in  the  well.  This  is  the  assa^tfoh. 
that  enables  one  to  incorporate  the  extra  bound  electrons” 
into  the  Fermi  sea.  The  Wentzel-Kramers-Brillouin 
(WKB)  approximation  for  eigenvalues  predicts  that  the." 
number  of  bound  states  n  in  a  semiclassical  well  is  given  by* 


-  2m*  Viz) 


■■Ti'T'.v  •• 


Thus  we  see  that  the  number  of  bound  states  n-~-Jm*.  We 
could  then  argue  that,  in  the  GaAs  system  with  a  very  small 
effective  mass,  one  has  to  use  a  self-consistent  model  instead 
of  the  Thomas-Femri  model  since  the  number  of  expected 
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Fig.  1.  Representative  wave  functions  of  accumulation  layeis.  is  the 
potential  drop  in  the  GaAs  layer.  Note  the  decrease  in  theampHtnde  of  ^ 
near  the  barrier.  The  respective  energies  of  the  electrons  are  indicated  in 
meV.  These  curves  are  normalized  to  have  an  asymptotic  form  of 
2  sin^fqr  ~rj)  Tar  away  ftom  the  barrier.  The  energies  are  measured  with 
respect  to  the  conduction  band  edge  at  infinity 


bound  states  in  the  well  would  be  small.  The  discovery  that 
«  =  1  verifies  this  hypothesis  and  further  indicates  the  im- 
portpice  of  using  a  self-consistent  model. 

In  Figs.  1(a)  and  1(b),  we  have  shown  the  accumulation 
layer  self-consistent  wave  functions  of  two  typical  unbound 
states.  The  doping  density  is  10”  cm~^  in  Fig.  1  (a)  and  10” 
cm“^  in  Fig.  1  (b).  These  are  the  wave  functions  of  the  mo¬ 
bile  states  when  a  potential  difference  of  Vo  is  allowed  to 
drop  in  the  GaAs  layer.  In  the  graphs,  the  normalization  of 
these  states  is  done  in  such  a  way  that  far  away  from  the 
barrier  the  value  of  approaches  2  sin*(^ — 7).  This 
normalization  scheme  aUows  us  to  compare  the  rdative  dis¬ 
tortions  that  occur  in  each  wave  function  as  a  function  of 
energy  and  distance.  In  all  these  graphs,  we  definitely  see  a 
lack  of  negative  charge  near  the  barrier.  The  decaying  expo¬ 
nential  tail  of  these  states  gives  rise  to  the  turmeling  current 
on  the  other  side  of  the  barrier.  Thus  it  is  very  important  to 
know  the  amplitude  of  these  wave  functions  at  the  interface 
to  accurately  account  for  the  tunneling  current.  We  also  see 
that,  as  the  energy  of  the  electrons  is  increased,  their  distor¬ 
tions  from  plane  wave  states  is  reduced.  Thus  the  low  lying 
conduction-band  electronic  states  are  the  ones  that  are  most 
distorted  by  the  potential.  Each  solution  here  is  constructed 
with  an  incoming  wave  and  a  reflected  wave  that  is  almost 
exactly  in  anti-phase  with  each  other.  (The  phase  difference 
would  be  exactly  a-  if  the  barrier  was  infinite.)  The  definite 
phase  difference  between  these  two  waves  plays  an  impor¬ 
tant  role  m  giving  rise  to  a  self-consistent  charge  density  that 
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Fig.  2.  Corresponding  carrier  densities  obtained  from  the  wave  functions  of 
Fig.  1 .  The  self-consistent  (SC)  carrier  density  is  the  sum  of  the  two  curves 
mobile  +  background,  and  bound  carrier  densities.  Notice  that  the  self- 
consistent  curve  is  very  different  &om  the  Thomas-Fermi  prettiction.  There 
is  an  excess  of  negative  carriers,  S.7  X 10”  cm“*  in  case  (a)  and  4.9X  10" 
cm"*  in  case  (b). 


has  Friedel  type  osciHatioDS.  However,  when  a  scattering 
event  occurs,  the  phase  difference  between  these  two  waves 
would  be  randomized.  Thus  scattering  would  lead  to  an  ex¬ 
ponential  decay  in  the  Friedel  oscillations  of  charge  density. 
The  length  scale  of  this  phenomenon  can  be  estimated  by 
assuming  that  scattering  at  low  temperatures  is  predomi¬ 
nantly  charged  impurity  scattering.  Realistically  we  can  rep¬ 
resent  ^,^  =  1  -  cos(2^z  -  277)6-*''^  in  this  case.  Thus 
the  definite  phase  relationship  due  to  the  reflection  firem  the 
barrier  would  decay  beyond  a  typical  scattering  length  A 
away  from  the  barrier.  A  calculation  we  did  with  a  crude 
model  suggests  that  .  As  these  figure  also  fflus- 

trate,  higher  energy  electrons  penetrate  more  into  the  bar¬ 
rier,  and  thus  would  play  a  more  important  role  in  current 

'  conduction.  •  - .  ' 

*  In  Figs.  2(a)  and  2(b)  the  dashed  curve  denotwl  by 
“bound”  shows  the  bound  state  contribution  to  the  charge 
density.  All  curve  shown  here  correpond  to  the  wave  func¬ 
tions  of  Fig.  1.  This  bound  charge  arise  from  the  two-di¬ 
mensional  subband  of  electrons  that  is  bound  in  the  GaAs 
'  accumulation  wclL  This  curve  also  resemble  the  shape  of 
the  wave  function  of  the  bound  state,  since  these  two 

quantitie  are  proportional  to  each  other.  The  average  ^a- 
tial  extent  of  these  bound  state  is  about  2(X)-300  A,  and  the 


peak  of  the  bound  state  charge  density  occurs  arotmd  50-100 
A.  These  state  show  the  expected  trend  that  smaller  binding 


energie  reults  in  state  that  have  a  larger  spatial  etent.  The 
mobile  state  charge  density  plus  the  positive  background 
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charge  density  is  shown  as  the  dashed  curve  denoted 
“mobile  +  background.”  Thus  when  we  are  sufficiently  far 
away  from  the  barrier,  the  mobile  state  charge  density  and 
the  positive  background  cancels  each  other  out  to  give  a  neu¬ 
tral  medium.  One  noticeable  feature  of  the  mobile  state 
charge  density  is  that  it  is  depleted  near  the  interface.  Thus 
the  dashed  curve  has  a  net  positive  charge  density  due  to  the 
dominant  positive  background  near  the  interface.  The  rea¬ 
son  for  this  lack  of  charge  very  near  the  barrier  can  be  under¬ 
stood  in  terms  of  the  boundary  conditions  that  dictate  expo¬ 
nential  decay  of  the  wave  function  inside  the  barrier.  Thus, 
in  accumulation  layers  as  well  as  in  depletion  layers,  the 
interference  of  the  incoming  and  the  reflected  waves  pro¬ 
duced  only  a  small  amplitude  of  the  wave  fimetion  at  the 
barrier.  This  is  one  of  the  major  contrasts  to  the  Thomas- 
Fermi  theory  of  screening  which  predicts  a  charge  density 
that  also  dips  down  to  a  maTitniiTn  negative  value  at  the 
interface.  Figures  2(a)  and  2(b)  clearly  show  this  difference 
between  the  charge  densities  of  the  Thomas-Fermi  and  self- 
consistent  models.  However,  another  reason  for  the  lack  of 
charge  is  that  the  mobile  states  have  to  be  orthogonal  to  the 
bound  state  in  this  region.  Thus,  the  presence  of  the  bound 
state  in  the  well  expels  the  mobile  state  electrons  further 
away  from  the  weD,  causing  a  lack  of  charge  (a  hole-like 
state)  in  the  vicinity  of  the  boimd  state.  This  is  purely  an 
effect  due  to  the  orthonormalization  of  the  mobile  states 
with  respect  to  the  bound  state.  Both  these  effects  gives  rise 
to  depletion  of  the  mobile  electrons  away  from  the  well  re¬ 
gion  even  in  the  accumulation  layer.  These  curves  also  un¬ 
dergo  abrupt  jumps  as  they  cross  the  barrier  and  enter  the 
undoped  AlAs  side.  This  is  merely  an  artifact  of  the  disap¬ 
pearance  of  the  background  positive  charge  as  we  go  from 
the  doped-GaAs  side  to  the  undoped-Al,Ga,  _^As  side.  We 
are  also  able  to  observe  the  fraction  of  the  electron  charge 
density  that  has  penetrated  inside  the  barrier  in  these  curves. 

The  solid  curve  denoted  SC  (self-consistent)  is  the  sum  of 
the  charge  densities  due  to  the  positive  background,  bound 
electrons,  and  the  mobile  electrons,  and  is  proportional  to 
/7u,(  (z) .  This  is  the  charge  density  that  is  used  in  the  Poisson 
equation  to  determine  the  self-consistent  potentiaL  Another 
interesting  featureofthe  self-consistent  charge  density  curve 
in  Fig.  2(a)  is  that  there  is  a  dipole  layer  at  the  interface  that 
has  an  extent  of  about  20  A  on  each  side.  However,  in  the  low 
doping  case  of  Hg.  2(b),  we  do  not  see  the  formation  of  such 
a  dipole  layer.  The  excess  negative  charge  that  produces  the 
dip  in  the  potential  mainly  comes  from  the  bound  state.  The 
bound  state  charge  density  in  this  region  is  2-S  times  the 
magnitude  of  the  background  doping  in  Fig.  2(b).  By  com¬ 
paring  Figs.  2(a)  and  2(b)  we  also  see  that  the  rdative  im¬ 
portance  of  the  bound  state  in  accounting  for  the  charge 
density  is  greater  in  the  low  density  10’’  cm“^  case  than  in 
the  high  density  10**  cm  case. 

In  Figs.  3(a)  and  3(b)  we  see  the  potentials  correspond¬ 
ing  to  the  wave  functions  of  Fig.  1,  and  charge  densities  of 
Fig.  2.  The  exceptional  closeness  of  the  self-consistent  poten¬ 
tial  and  the  Thomas-Fermi  potential  is  indeed  surprising, 
considering  that  the  resulting  charge  densities  are  quite  dif¬ 
ferent.  One  way  of  understanding  this  is  to  say  that,  the  pres¬ 
ence  of  the  bound  state  in  the  well  alters  the  unbound  states 
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Fig.  3.  Corresponding  poteanal  for  Figs,  land  2  above.  Nonce  the  sunilari-  Fio  1  ^representative  wave  functions  of  depletion  layers.  K,  isthepoten- 
ties  between  the  Thomas-Fami  (TF)  and  self^^nsistent  (SC)  ptkentials,  liai  crtop  in  the  GaAs  layer.  Note  that  these  states  resemble  plane  waves  in 

in  the  accumulation  layers.  The  bound  states  are  very  loosely  bound,  with  the  depletion  region.  unBlte  in  the  accumulation  regioa  of  Fig.  1.  The  re¬ 
binding  energies  of  5.8  meVin  (a)  and  11.6  meV  in  (b).  Theapplied  poten-  spective  energies  of  the  electrons  are  inriirftfrf  in  meV.  These  carves  are 

dal  drop  is  50  mV  in  (a)  and  (b).  i>o'™>li2a<ltohaveanasymptoticformof2sin’(4Z  — i7}.&raway&omthe 

barrier.  The  energies  atemeasured  with  respect  to  the  conduction  band  edge 

.at  infinity. 


wave  functions  so  as  to  produce  an  orthonormalization  hole 
in  the  conduction  band.’’*  *  This  orthonormal  hole-like  state 
that  exists  in  the  conduction  band  moves  in  correlation  with 
the  bound  state  in  response  to  external  biases  so  as  to  make 
the  total  charge  in  the  layer  a  continuous  function  of  the 
applied  bias.  An  analogous  result  has  been  proved  by  Kohn 
and  Majumdar  even  in  the  event  when  the  number  of  bound 
states  in  the  well  is  variable.  Thus,  if  a  bound  state  vanishes 
.  at  a  certain  strength  of  the  applied  bias,  the  significant  frac¬ 
tion  of  charge  that  vanishes  with  it  reappears  in  the  conduc- 
'■y;  tion  band  as  unbound  charge.  This  corresponds  to  the  rear- 
Vj  rangement  of  charge  that  takes  place  due  to  the  vanishing  of  ' 

the  orthonormal  hole.  This  enables  one  to  lump  the  behavior 
,  of  the  bound  state  and  the  orthonormal  holelike  state  from  -. 

'  the  continuum  as  a  simple  excess  of  charge  thaTmonotoni- 
.  cally  increases  with  the  applied  bias.  This  reduces  to  the  ' 
'  same  type  of  behavior  that  the  Thomas-Fermi  theory  exh^  :  - 
.  its.  This  form  of  behavior,  added  to  the  constraint  that 
derivatives  of  the  self-consistent  and  Thomas-Fermi  pbtai-j!;- 
'•■Jtials  be  equal  at  the  interface  (condition  of  same  exc« 
^’'•'charge),  would  explain  why  the  potentials  look  very  similar.__'l 
^;^!5j-The  position  of  the  bound  state  that  occurs  in  these  wells*is 
shown  as  the  horizontal  line  segment  in  the  figure. 


;  B.  Depletion  regions  " 

■  .  '  In  Figs.  4(a)  and  4(b}  we  see  the  wave  functions  of  typic^ 
depletion  layers.  These  are  also  normalized  in  the  same. way 
as  in  Figs.  1(a)  and  1  (b),  to  obtain  their  asymptotic  forms  of 
2  sin’(92  — 17)  to  the  far  left-hand  side  of  the  barrier.  Results 


are  presented  for  a  doping  density  of  10”  cm*’  in  Fig.  4(a) 
and  10”  cm~’  in  Fig.  4(b).  These  electrons  are  reflected 
before  they  hit  the  barrier  due  to  the  repulsive  self-consistent 
potential  that  bends  upwards  in  these  layers.  Thus  the  com¬ 
ponent  of  their  wave  functions  that  penetrate  into  the  barrier 
is  vastly  reduced  fimm  those  in  the  armimiilatirm  layers.  We 
also  see  that  these  solutions  are  much  less  cfistorted  from 

plane  waves  than  in  the  case  of  the  accumulation  layen.  This 

lack  of  distortion  can  be  attributed  to  the  larfr  of  the  bound 
state  on  this  side  of  the  barrier.  Thus  there  is  no  distortion 
due  to  the  orthomamalization  requirement.  One  ntkw  no¬ 
ticeable  result  is  that  these  waves  actually  undergo  a  slight 
accumulation  in  the  first  half  wavelength  away  from  the  bar- 
tier.  Again  this  leads  to  Friedel  oscillations  in  the  charge 
density.  - 

Figures  5(a)  and  5(b)  show  the  correspon fling  charge 
,  densities  of  the  wave  functions  of  Hg.  4.  The  solid  line  SC  is 
'the  self-consistent  diarge  density  which  is  the  rf  the 
r  mobile  electron  contribution  and  the  backgronml  positive 
ion  contribution.  In  Fig.  5(a)  a  very  small  potential  drop 
,(20  meV)  is  appEol  with  a  deficient  charge  dasity  of 
9.0x  10”  cm*’.  Ob  the  other  hand,  in  Fig.  5(b)  a  largin' 
.potential  drop  (40  meV)  compared  to  the  Fermi  energy 
—  12meV)  isqiplied  with  a  deficient  charge  density  of 
2.2  X 10'*  cm  ’.  Here  we  clearly  see  the  long-range  Friedel 
oscillations  in  the  charge  density  that  are  characteristic  of 
self-consistent  calculations.  They  do  not  decay  with  an  expo- 
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Fig.  5.  Corresponding  carrier  dmsities  obtained  from  the  wave  functions  of 
Ftg.  4.  The  selfconsistcnt  (SC)  carrier  density,  and  the  Thomas-Fenni 
(TF)  carrier  density  are  quite  different.  Notice  also  the  character  of  the 
long  range  oscillations  of  the  seif-consistent  charge  density.  There  is  a  defi¬ 
ciency  of  negative  carriers,  9.0x10"  cm"’ incase  (a)  and 2.2x10“ cm"’ 
incase  (b). 

nential  decay  envelope  although  their  amplitudes  are  small. 
These  oscillations  are  dominant  at  high-doping  density, 
email  applied  positive  bias,  and  low  temperatures.  From  a 
physical  point  of  view  the  long-range  oscillations  in  Ae 
screening  charge  arises  from  a  sharp  Fermi  surface.  The  rea¬ 
son  is  that  it  is  not  possible  to  construct  a  smooth  function 
free  of  oscillations  with  the  restricted  set  of  waves  g<*jr. 
Thus,  the  screening  charge  will  always  exhibit  osdDations 
approximately  at  2kj,,  twice  the  cutoff  frequency.  The  s^-  . 
consistent  screening  charge  also  has  a  decay  envelope  ttat  . 
dies  as  an  inverse  power  law,  in  contrast  to  the  Thomas^ 
Fermi  result  that  predicts  an  exponential  decay  rate. 
ever,  at  finite  temperatures,  these  Friedel  oscillations  die  out  ■ 
with  an  exponential  factorproportional  to  A:*  r.  Thus,  in  r^ 
systems  at  finite  temperature,  the  signifidmce  of  these  o^- 
lations  is  negligible.  The  figures  compare  the  self-conas^t_ 
charge  density  with  the  Thomas-Fenni  charge  den^.  In 
the  My  depleted  region  which  occurs  in  Fig.  5(b),  these  " 
two  charge  densities  resemble  one  another.  However,  injhe 
‘ ;  region  beyond  thia,  the  difference  is  significant  This  is  . .  • 

. '  ly  due  to  the  difference  in  the  long-range  screening  prop^^  > 
ties  predicted  by  the  two  theories.  There  is  also  no  apparent 
dipole  layer  on  this  side  of  the  barrier  unlike  in  the  accmu- 
lation  layers. 

-  In  Figs.  6(a)  and  6(b)  the  self-consistent  potential  com-  • 
sponding  to  the  wave  functions  of  Fig.  4,  and  the  charge 
density  of  Fig.  5  is  shown.  The  self-consistent  potential  and 
the  Thomas-Fenni  potential  for  the  depletion  layers  can  ^ 
significantly  different  firom  each  other.  The  agreement  . 


Fig.  6.  Gimsponding  poteudal  for  Figs.  4  and  S  above.  Note  >lat  in  deple- 
don  layen,  the  selfcoosistat  and  Thomas-Fermi  potentials  can  look  very 
different  from  each  other,  as  illustrated  in  case  (a).  The  applied  potential 
drop  is  20  mV  in  (a)  and  40  mV  in  (b). 

between  these  two  curves  is  worse  in  depletion  layers  than  in 
accumulation  layers.  The  reason  for  this  difference  is  the 
assumption  that  the  charge  gets  depleted  in  an  exponential 
manner  with  a  Thomas-Fermi  decay  length  is  not  a  valid 
one.  In  fact,  the  screening  charge  decays  in  the  Hartree  case 
with  an  envelope  fhnetion  of  an  inverse  pmwer  law  as  dis¬ 
cussed  above.  We  also  discovered  that  the  Thomas-Fermi 
solution  lies  above  the  self-consistent  solution  in  highly 
doped  cases,  and  below  the  self-consistent  solution  in  lightly 
doped  cases. 

jV.  CONCLUSION 

In  conclimon,  we  see  that  the  Thomas-Fenni  theory  is 
quite  good  at  describing  the  potential  in  accumnlation  layers 
although  it  does  a  poor  job  at  describing  the  change  dosities. 
(?ar  results  indicate  Aat  by  assuming  the  actual  potential  to 
be  given  by  the  Thomas-Fenni  form,  and  solving  fin'  its  ei¬ 
genstates  and  eigenvalues  one  can  obtain  answers  that  are 
very  close  to  the  s^-consistent  ones.  Thus,  calculating  the 
charge  density  from  the  wave  functions  of  the  Thomas-Fer- 
mi  potential  is  quite  adequate  for  describing'the 'current 
through  the  device,  and  calculating 'the  po^om the 
bound  levels.  Hence;  one  does  not  have  to  go  throng  a  full¬ 
blown  numerical  seff-consistent  calculation  to  obtain  rea¬ 
sonable  results  for  the  accumulation  layers.  On  the  other 
hand,  the  potential  and  charge  density  of  the  depletion  layers 
can  be  quite  different  from  the  self-consistent  one,  and  hu  to 
be  trrated  more  <ardhlly.  The  results  of  this  study  also  indi¬ 
cates  that  there  are  no  other  electron  levels  in  the  well  of  the 
accumulation  region,  except  for  a  single  loosely  bound  level. 
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Effects  of  barrier  phonons  on  the  tunneling  current  in  a  doable'barrier  structure 

G.  Y.WuandT.  C.  McGiU 

T.  J.  Watsoiu  Sr.  Laboratory  of  Applied  Physics,  California  Institute  ofTechnolo^,  Pasadena,  California  91125 
(Received  27  August  1987;  revised  manuscript  received  20  May  1988) 

The  efects  of  AlAs  acoustical  and  optical  phonons  on  the  tunneling  current  in  an  ideal  GaAs- 
AlAs-GaAs-AlAs-GaAs  structure  are  discussed.  The  transfer  Hamiltonian  method  was  extended  to 
inelastic  tunneling  current  in  a  double-barrier  structure.  It  is  found  that  the  current  off  resonance 
in  the  Tsu-Esaki  model  could  be  enhanced  by  orders  of  magnitude  by  inelastic  tunneling  due  to  the 
coupling  of  electrons  to  barrier  optical  phonons.  The  contribution  to  the  current  due  to  the 
deformation-potential  coupling  of  electrons  to  barrier  acoustical  phonons  is  found  to  be  much  less 
important. 


Semiconductor  heterostructures  involving  GaAs  and 
AlAs  have  been  the  subject  of  both  theoretical  and  exper¬ 
imental  studies.’"*  In  GaAs-AlAs-GaAs-AlAs-GaAs 
structures,  AlAs  layers  act  as  energy  barriers  while  the 
middle  GaAs  layer  acB  as  a  well.  Of  particular  interest 
are  GaAs-AlAs-GaAs-.AlAs-GaAs  double-barrier  struc¬ 
tures  which  are  specifically  doped  as  n-i-i-i-n.  Observa¬ 
tions  of  resonant  tunneling  of  electrons  through  double- 
barrier  structures  have  been  reported.’  Tsu  and  Esaki^ 
have  theoretically  treated  the  general  case  of  electronic 
tunneling  through  multiple  barriers.  In  the  double¬ 
barrier  case,  resonant  tunneling  was  shown  by  them  to 
give  rise  to  a  current  maximum  iJpcak^  ^be  voltage  bias 
( Vg )  where  the  Fermi  sea  of  an  electrode  is  aligned  in  en¬ 
ergy  with  one  of  the  quasibound  states  in  the  GaAs  welL 
The  current  drops  very  rapidly  at  other  voltages.  Subse¬ 
quent  works  have  refined  the  theory  for  calculating  the 
barrier  transmission  (D  and  current- voltage  {/-K,) 
characteristics.*"*  Typically,  the  peak  current  they  c^- 
culated  agrees  with  the  measured  value.  However,  the 
measured  valley  current  is  greater  than  the  theoretical 
value  by  at  least  1  order  of  magnitude.  Recently,  Frens- 
ley  has  employed  a  transport  theory  approach  which  in¬ 
cludes  losses.^  However,  all  of  these  theories  only  consid¬ 
er  elastic  tuimeling.  On  the  other  hand,  Tsu  and  Dohler* 
have  considered  inelastic  tunneling  assisted  by  the  weQ 
phonons.  In  a  recent  study  of  GaAs-AlAs-GaAs  struc¬ 
tures,  Collins,  Lambe,  and  McGill’  reported  the  observa¬ 
tion  of  inelastic  tunneHng  of  electrons.  In  such  tunneling, 
the  excitation  of  AlAs  phonons  could  take  away  finite 
amount  of  energy  and  momentum  from  the  tunneling 
electron.  Hence,  the  total  energy  and  the  tran^erse 
momentum  of  the  tunneling  electron  are  no  longer  con¬ 
served.  As  we  shall  show,  they  could  enhance  the  current 
off  resonance  by  orders  of  magnitude  in  the  Tsu-Esaki 
model.  This  suggests  the  importance  of  the  inclusion  of 
inelastic  tunneling  assisted  by  barrier  phonons  in  calcula¬ 
tions  such  as  the  self-consistent  analysis  of  Onishi  et  al.^ 
Here,  we  present  the  first  theoretical  study  of  the 
effects  on  the  tuimeling  current  in  a  double-barrier  struc¬ 
ture  due  to  the  electron  interaction  with  barrier  phonons. 
The  transfer  Hamiltonian  method  proposed  by  Bardeen’® 
originally  for  the  treatment  of  the  tuimeling  of  electrons 
through  a  single-barrier  structure  b  extended  to  calculate 

40 


phonon-induced  inelastic  tunneling  in  a  double-barrier 
structure.  For  illustration,  we  shall  consider  the  zero 
temperature  case  where  only  phonon  emission  b  possible. 
In  this  case  with  proper  approximations,  we  obtain 
analytical  expressions  which  shed  some  light  on  the  im¬ 
portant  effects  of  inelastic  tunneling.  These  effects 
change  I-V characteristics  for  a  double-barrier  structure. 

We  treat  a  standard  double  barrier  structure  consbting 
of  a  layer  of  GaAs  followed  by  a  barrier  of  Ga,_,Al,As 
followed  by  a  well  of  GaAs  followed  by  a  second  barrier 
of  Ga,_,Al,As  followed  by  the  contact  layer  of  GaAs. 
The  interfaces  between  the  layers  are  at  Xj,  x^,  Xj,  and 
X4,  respectively.  The  barrier  height  (the  conduction-band 
offset)  is  taken  to  be  55%  of  the  band-gap 

difference.”  The  cross-section  area  b  i4.  Each  electrode 
has  the  thickness  L.  The  effective  mass  of  the  electron  b 
taken  to  be  m  *  whether  in  GaAs  or  AlAs.  The  left  bar¬ 
rier  and  the  right  barrier  are  taken  to  be  equally  thick: 
d|  =d2—d,  u;  b  the  width  of  the  well.  The  voltage  bias 
Vg  b  such  that  the  lowest  quasibound  level  E^  b  lower  in 
energy  than  the  conduction-band  edge  of  the  left  elec¬ 
trode. 

We  use  the  transfer  Hamiltonian  method  proposed  by 
Bardeen’®  and  extend  the  method  to  treat  the  inelastic 
tunneling  of  an  dectron  through  the  excitation  of  an 
AlAs  phonon.  We  choose  the  states  and  so  that 
is  matched  to  the  correct  solution  for  x  bat  decays 
in  the  region  x  <X|  instead  of  satisfying  the  wave  equa¬ 
tion,  and,  similar^,  tpi  continues  to  decay  for  x^X2. 
Then  b  a  correct  solution  for  the  Hamiltonian  H  for 
X  ^X]  and  is  correct  for  x  ^X2.  With  WKB  approxi¬ 
mation,  we  have 


^fr,(x)= 


^,(x)= 


2 

AL 

1 


in 

e****sin{k,x  -f-y/), 

Oc.x. 

e  •  le  '  ‘ 


X£Xi 


x^x, 

where  k,  is  the  wave  vector  of  the  electron  in  the  left 
electrode,  ix;  b  the  imaginary  wave  vector  in  the  left  bar¬ 
rier,  and  Xi(x)  b  defined  in  the  equation.  Similarly, 


9969 


®  1989  The  American  Physical  Society 


9970 


BRIEF  REPORTS 


40 


1/2 


^,(x)= 

2 

AL 

e 

1 

^“'”sin(k,x -t-y,),  x>x^ 

^,(x)= 

L 

2 

1/2 

hL 

1/2  it, jt,  ~itAx—x,) 

e  *  *e 

.41, 

{ [sm{k^w)Ye  t4cos(k„ui)]^c^*‘*^  j 

1/2 

*- 

= 

1 

.41, 

e 

X<X2 

r 


(2) 


where  is  the  wave  vector  in  the  right  electrode,  iV,  is 
the  imaginary  wave  vector  in  the  left  barrier,  k„  is  the 
wave  vector  in  the  well,  is  the  imaginary  wave  vector 
in  the  right  barrier,  and  X,(x)  is  defined  in  the  equation. 
The  electron-phonon  coupling  in  the  left  AlAs  barrier 
gives  arise  to  the  interaction  Hamiltonian 


A. 


AIM 


(3) 


for  the  deformation-potential  (DP)  coupling  due  to  the 
longitudinal-acoustic  (LA)  phonon  mode  qX,  where  q  is 
the  phonon  wave  vector  and  A.  specifies  the  polariza¬ 
tion.’^  Here,  Dj-  is  the  deformation  potential  for  F- 
valley  electrons,  p  is  the  density,  o,  is  the  sound  velocity 
in  AlAs,  and  P^aias  “  volume  of  the  left  barrier.  On 
the  other  hand. 


ff 


a  _ 
e-ph 


Airi 

9 


fie> 

1 

'  «0 

(4) 


for  the  polar-optical  (PO)  coupling  due  to  the- 
longitudinal-optical  (LO)  phonon  mode  qA.  ejj  and 
are  the  dielectric  constants  at  zero  and  optical  freqnen- 
_ —J 


cies,  respectively.  Other  kinds  of  electron-phonon  cou¬ 
pling,  piezoelectric  coupling  for  example,  are  much  weak¬ 
er  and  much  less  important.’^  They  are  not  considered 
here.  The  transition  rate  can  be  calculated  easily  by  the 
application  of  Fermi’s  golden  rule.  The  temperature  is 
taken  to  be  zero  to  simpliiy  our  analysis;  hence,  only  pho¬ 
non  emission  needs  to  be  considered.  The  inelastic  tun¬ 
neling  current  is 

"  q,A.  k  k./ 

X[l— /(e*-)]5(et-f-cF'„— Ei— Aa)  , 

(5) 

where  fio)  is  the  phonon  energy,  and  Ej^'  are  the  kinetic 
energies  of  the  electrons  in  the  left  and  right  electrode, 
respectively,  and  Mis  the  matrix  element.  Here,  k  and  k' 
are  wave  vectors  of  the  electron  at  left  and  right  elec¬ 
trodes,  respectively,  and  /(e)  is  the  Fermi  distribution 
function.  The  general  expression  for  7^^  at  a  finite  tem¬ 
perature  may  be  derived  with  the  use  of  many-body 
theory.  Bennett  et  al.  have  given  the  derivation  of  the 
formula  for  the  single-barrier  case.  We  do  not  attempt  to 
treat  finite  temperature  case,  however.  Replacing  2  with 
f,  we  have 


'ioel 


..TX8(k||-ky-q|i)(2ir)*6(Etq-eFa-et.-fitu) . 


(6) 


Here,  we  have  '''■  ■ 

U^=47re[n(a(l/eo-l/€j/SirY^/q  ... 
for  PO  coupling,  and  /J: 

U,=Dri1iq/2pv,)^^  '«) 

for  DP  coupling.  Tis  the  overlap  integral 

Suppose  now  the  dectron  is  scattered  into  the  qoari- 
bound  level,  then  m  the  left  barrier  would  be  enhanced 
by  the  exponential  factor  This  can  be  verified  by 

setting  cos(k„ti;)  to  zero  in  Eq.  (2)  when  near  resonance. 
This  would  in  turn  increase  T  by  the  same  factor.  There¬ 
fore,  the  contribution  to  the  inelastic  current  due  to  mch 


a  process  would  be  dominant.  It  turns  out  that  for  PO 
coupling,  we  have 

(2ir)3lpK5  *  *  2 

_  fc* 

V  go  ^ 

For  DP  coupling,  we  have 

2Aem*^d\  vk^  D\ft 

(lv?fi?kl‘  2  2ir*  2pv,  , 

X{dr'(‘^r*+9o)‘'^+9oS«di~'[(go<^|)"'])  • 

(11) 
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Here,  £4  is  the  zero-bias  barrier  height,  Kq 
=(2m*E^/fi^)'^,  Vi=2E^/K^i,  and  satisfies 
^g^/2m*=eP;/2-&)-£,. 

We  now  compare  contributions  to  the  current  from  the 
elastic  and  the  inelastic  process.  We  consider  only  the  in¬ 
elastic  process  induced  by  PO  coupling  or  DP  coupling. 
Inelastic  processes  induced  by  other  couplings  will  fur¬ 
ther  increase  the  inelastic  current.  Therefore,  the  inelas¬ 
tic  current  we  consider  here  is  a  lower  bound  of  the  actu¬ 
al  value.  However,  even  without  a  complete  knowledge 
of  the  inelastic  current,  we  can  still  demonstrate  the 
dramatic  change  in  order  of  magnitude  of  the  current 
which  is  due  to  the  inelastic  processes. 

In  Fig.  1,  we  plot  the  lower  limit  of  the  ratio 

of  inelastic  current  to  elastic  current,  versxis  barrier 
thickness  for  biased  symmetric  structures  with  u)=50  A, 
£y =50  meV,  and  V^=0A  V.  Both  LO  and  LA  phonon- 
induced  effects  are  shown.  For  thin  barrier  cases  where 
d=20  A,  polar  coupling  induces  an  inelastic  current 
comparable  to  the  elastic  current.  As  can  be  seen, 
deformation-potential  coupling  has  much  smaller  effects 
than  polar  coupling.  The  reason  is  that  the  PO  coupling 
is  much  stronger  than  the  DP  coupling.  We  may  neglect 
the  contribution  &om  DP  coupling  in  comparison  to  that 


times  a  s  large  as  the  elastic  one. 

In  Fig.  2,  we  plot  the  upper  limit  of  the  ratio  of  the  res¬ 
onant  current,  which  occurs  at  K,  =0.3  V,  to  the  current 
at  Vg=0A  V,  wluch  occurs  off  resonance,  versus  barrier 
thickness.  The  dotted  line  is  obtained  with  Tsu-Esaki 
model,’^  which  includes  only  elastic  process.  In  that 
case,  the  ratio  shown  here  is  equal  to  the  peak-to-vaUey 
ratio  of  the  J-V  carve.  The  solid  line  is  obtained  with  in¬ 
elastic  tunneling  included.  The  peak  current  is  mostly 
due  to  the  elastic  tunneling,  since,  at  resonance,  the  elas¬ 
tic  current  is  much  larger  than  the  inelastic  current.  The 
current  off  resonance  is  largely  due  to  the  inelastic  tun¬ 
neling,  since  the  inelastic  current  is  much  larger  than  the 
elastic  current,  as  shown  in  Fig.  2.  The  ratios  shown  by 
the  two  curves  are  comparable  for  thin  barrier  cases. 
However,  as  the  barrier  l^omes  thicker,  the  solid  curve 
only  varies  slowly,  since  both  Jj^y  and  have  same 

exponential  dependence.  For  barrier  thickness  equal  to 
40  A,  magnitudes  of  the  two  ratios  differ  by  7  in  the  natu¬ 
ral  log  scales  The  theory  with  the  inelastic  process  in¬ 
cluded  predicts  a  much  smaller  value  of  ratio  for  thick 
barrier  cases. 


from  PO  coupling.  Effects  of  both  coupling  increase  as 
the  barrier  becomes  wider,  due  to  the  exponential 
enhancement  in  the  right  wave  function  rb^-  For  the 
thick  barrier  case  where  d=40  A,  polar  coupling  even 
gives  rise  to  an  inelastic  current  which  is  a  thousand 
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FIG.  3.  The  upper  bound  of  Jfak/J-^nOet-y  of 

peak  current  at  K,=0J  V  (at  resonance,  the  peak  in  the 
current-voltage  case)  to  the  current  at  V,  =0.4  V  (off  resonance, 
the  valley  in  the  case  of  elastic  transport  but  not  necessarily  in 
the  case  with  inelastic  processes)  vs  Fermi  energy^  for  the  sym¬ 
metric  structure  with  a  barrier  thickness  of  40  A  and  a  well 
width  of  50  A. 


In  Fig.  3,  the  upper  bound  of  the  ratio  is  shown  versio 
Fermi  ener^  for  the  symmetric  structure  with  J=40  A 
and  in =50  A.  The  solid  curve  includes  contribution  from 
the  inelastic  process  and  the  dashed  curve  includes  only 
contribution  from  the  elastic  process.  The  curves  are 


shown  for  the  region  from  Ef  =  5  meV  to  Ey=50  meV 
corresponding  to  dopant  density  from  10‘Vcm^  to 
10'*/cm^  currently  used  in  the  tunneling  experiments. 
We  see  that  both  curves  behave  similarly  as  the  Fermi  en¬ 
ergy  changes.  The  difference  between  them  is  maintained 
through  the  Fermi  energy  range  of  interest.  This  shows 
the  importance  of  the  inelastic  process  for  both  low  and 
high  doping  cases. 

In  summary,  the  inclusion  of  the  inelastic  tunneling  in 
the  theory  for  double-barrier  structures  is  very  impor¬ 
tant.  We  have  studied  specifically  the  inelastic  tunneling 
induced  by  phonons.  Two  types  of  coupling  have  been 
considered.  The  PO  coupling  is  much  stronger  tlian  the 
deformation-potential  coupling,  and  hence  the  current  in¬ 
duced  by  PO  coupling  is  much  larger  than  that  by  DP 
coupling.  Because  of  similar  reasons,  we  expect  the  effect 
of  piezoelectric  (P£)  coupling  also  to  be  negibile  in  com¬ 
parison  to  that  of  PO  coupling. 

The  electron-phonon  interaction  induces  a  channel 
through  which  the  electron  can  tunnel  much  more  readi¬ 
ly  than  through  the  elastic  channel.  With  the  excitation 
of  a  barrier  phonon,  the  electron  can  utilize  the  quasi¬ 
bound  level  to  tunnel  through  the  barrier.  The  inelastic 
resonant  tunneling  enhances  the  current  off  resonance  by 
orders  of  magnitude.  The  effects  are  best  reflected  in  the 
big  difference  between  the  magnitudes  of  current  ratios 
which  are  predicted,  with  and  without  inclusion  of 
electron-phonon  coupling,  in  the  simple  model  of  Tsu 
and  Esaki.^  We  expect  the  inclusion  of  inelastic  processes 
to  be  important  in  any  refined  theory  such  as  the  self- 
consistent  analysis  of  Ohnishi^  or  Frensley's  transport  ap¬ 
proach.^  Other  mechanisms  such  as  impurity  scattering 
could  also  be  critical  to  current  transport  in  a  double¬ 
barrier  structure. 
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GaAs-ZnSe  and  InAs-ZnTe  double-barrier  resonant  tunneling  devices  are  proposed  and 
analyzed  theoretically.  These  structures  would  provide  a  novel  way  to  study  mixed  III-V/II-VI 
materials  systems  and  growth  techniques,  and  may  also  yield  devices  that  are  faster  and  have 
greater  pealc-to-valley  current  ratios  than  GaAs-AlGaAs  tunnel  structures.  In  the  GaAs-ZnSe 
system,  the  increase  in  speed  is  due  to  the  smaller  dielectric  constant  in  ZnSe,  compared  to  that 
in  AlGaAs.  In  the  InAs-ZnTe  system,  the  increase  in  speed  is  due  to  the  small  dielectric 
constant  in  ZnTe  and  the  high  mobility  in  InAs,  and  the  larger  peak-to-valley  ratio  is  due  to 
the  large  conduction-band  offset  that  is  expected.  We  calculate  current-voltage  curves  for  these 
devices  and  show  that  their  negative  differential  resistance  current-voltage  characteristics 
should  be  comparable  or  superior  to  those  of  QaAs- AlGaAs  devices.  The  effects  of  band 
bending  and  the  general  method  used  to  calculate  current  as  a  function  of  voltage  are  briefly 


discussed. 

Over  the  past  several  years,  much  of  the  study  of  device 
physics  and  technology  has  focused  on  III-V  materials  sys¬ 
tems.  However,  many  of  the  device  concepts  developed  and 
demonstrated  in  III-V  systems  can  be  applied  to  other  mate¬ 
rials  as  well.  One  device  that  has  been  the  subject  of  consider¬ 
able  interest  recendy,  both  for  studying  very  thin  semicon¬ 
ductor  heterostructures  and  as  a  possible  high-speed  device, 
is  the  double-barrier  resonant  tunneling  structure.  In  this 
letter,  we  discuss  theoretically  some  properties  of  GaAs- 
7nSf  and  InAs-ZnTe  double-barrier  resonant  tunneling 
structures;  this  is  the  Srst  time  mixed  III-V/II-VT  tunneling 
structures  have  been  proposed.  ZnSe  in  particular  is  a  poten¬ 
tially  important  11- VI  material  because  its  lattice  parameter 
is  very  close  to  that  of  GaAs  ( ~0.27%  mismatch),  and  be¬ 
cause  its  band  gap  is  very  wide  (2.67  eV  at  300  K),  which 
makes  visible  light  emission  possible.  Furthermore,  the  suc¬ 
cessful  growth  of  ZnSe/GaAs  double  heterostructures' 
would  seem  to  indicate  that  GaAs-ZnSe  double-bamer 
structures  can  be  grown  by  the  same  techmques.  GaAs-ZnSe 
and  InAs-ZnTe  tunnel  structures  would  in  fact  provide  a 
novel  way  to  study  mixed  III-V/II-VI  materials  systems  and 
growth  techniques. 

For  actual  device  applicadons,  one  advantage  that  a 
GaAs-ZnSe  structure  may  have  over  a  GaAs-AlGaAs  de¬ 
vice  is  speed.  The  time  for  an  electron  to  tunnel  across  the 
device  will  typically  be  on  the  order  of  10~'^— 10“'*  s,  and 
therefore  will  not  be  the  limiting  factor  in  the  speed  of  the 
device.^  A  parameter  that  may  be  more  relevant  to  the  device 
speed  is  the  RC  time  constant  for  charging  the  capacitor 
formed  by  the  barriers.  This  parameter  can  be  estimated  as 
r  =  ff/ A£,  where  A£  is  the  width  of  the  transmission  reso¬ 
nance  of  the  barrier-well  structure;  typically,  ris  on  the  ot- 
der  of  1  ps.  Since  the  capacitance  is  proportional  to  the  di¬ 
electric  constant  in  the  barrier,  and  since  21nSe  has  a  lower 
dielectric  constant  than  AlGaAs  (9.1  vs  12.0),  the  capaci¬ 
tance,  and  therefore  the  RC  time  constant,  should  be  re¬ 
duced  by  replacing  AlGaAs  with  ZnSe. 

InAs-ZnTe  tunnel  structures  should  exhibit  an  even 
greater  advantage  in  speed.  The  dielectric  constant  for  ZnTe 
is  10.4,  which  again  is  smaller  than  that  in  AlGaAs.  More 


important,  InAs  has  a  much  higher  mobility  than  GaAs 
[33  000  vs  8500  cmV(V  s)  at  300  K].  The  resistivity  in 
InAs  will  therefore  be  much  lower  than  in  GaAs,  leading  to  a 
large  reduction  in  the  RC  time  constant.  The  low  resistivity 
in  InAs  will  also  permit  a  greater  fraction  of  the  voltage  drop 
across  the  device  to  occur  in  the  barrier-well  region,  so  that  a 
smaller  bias  voltage  will  be  required  to  achieve  the  Tnarimiiin 
tunneling  current 

The  material  parameter  that  is  most  important  in  deter¬ 
mining  the  electrical  properties  of  the  double  barrier  is  the 
conduedon-band  offset  The  valence-band  oSset  for  the 
( 1 10)21nSe-GaAs  heterojunction  has  been  measured  by  x- 
ray  photoelectron  spectroscopy*  and  found  to  be,  depending 
on  growth  condinons,  between  0.96  and  1.10  eV  at  3(X)  K. 
The  value  of  0.%  eV  was  obtained  for  a  heterojunction  in 
which  ZnSe  was  deposited  on  GaAs  at  ~300  *C,  while  the 
value  of  1 . 10  eV  was  for  a  sample  in  which  ZnSe  was  deposit¬ 
ed  on  GaAs  at  23  ’C,  and  annealed  at  ~  300  *C  The  corre¬ 
sponding  conduction-band  offset  is  between  0.15  and  0.29 
eV,  which  is  comparable  to  the  direct-band  ofiet  in  the 
GaAs-Al^Ga)  _,As  system  (0.25  eV  for  x  =  0J3).  As  has 
been  seen  in  the  GaAs- AlGaAs  system,  band  o&et  values  in 
this  range  make  the  fabrication  of  high  quality  high  electron 
mobility  transistor  structures  difficult,  but  are  sufficient  for 
achieving  negative  differential  resistance  in  tunnel  struc¬ 
tures. 

In  the  InAs-ZiTe  system,  the  conduction-band  ofiet  is 
probably  very  large,  which  is  another  advantage  of  this  sys¬ 
tem  over  CjaAs-AlGaAs.  The  value  of  the  band  offset  has 
not  been  measured  directly,  but  we  can  obtain  some  idea  of 
its  value  if  we  assume  that  band  offsets  are  transitive.  For 
CaAs-2^nSe,  the  valence-band  of^  is  ==  1.1  eV;  for 
InAs-GaAs,  =  0.17  eV  (Ref.  4),  and  for  ZnTe-ZnSe, 
LEy  =  0.98  eV.*  Using  these  values,  we  obtain  =;0.30 
eV  or  EEc  ss  1.6  eV  for  InAs-ZnTe.  The  large  conduction- 
band  ofi^t  wDl  reduce  thermionic  currents  through  higher 
resonances  and  over  the  barrier,  facilitating  device  operation 
at  room  temperature. 

We  have  calculated  current-voltage  characteristics  for 
these  structures,  and  for  a  comparable  GaAs-AlGaAs  de- 
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vice  as  well.  For  the  GaAs-ZnSe  system,  we  have  assumed 
that  the  double-barrier  structure  consists  of  a  50  A  GaAs 
well  sandwiched  between  two  50  A  ZnSe  barriers.  The  well 
and  barrier  are  assumed  to  be  undoped,  and  the  contact  lay¬ 
ers  are  assumed  to  be  degenerately  doped  (10'*  cm"^) 
GaAs.  For  the  InAs-ZnTe  system,  the  structure  consists  of 
an  undoped  50  A  InAs  well  sandwiched  between  two  un¬ 
doped  20  A  ZnTe  barriers.  The  contact  layers  are  degener¬ 
ately  doped  ( 10'*  cm~^)  InAs.  Our  calculations  assume  the 
following  material  parameters  (for  T—  300  K):  for  GaAs, 
m*  =  0.066/n,,  £'^  =  1.42  eV,  £=12.9^0;  for  ZnSe,  m* 
=  0.16m,,  £,=167  eV,  £  =  9.1fo;  for  InAs,  m* 
=  0.028m,,  £,=0.36  eV,  £=14.6eo;  for  2^Te,  m* 
=  0.11m,,  £,  =126  eV,  e=  10.4£o.  These  parameters 
change  slightly  at  r=  77  K,  but  the  changes  are  not  suffi¬ 
cient  to  prevent  the  occurrence  of  negative  differential  resis¬ 
tance. 

Energy-band  diagrams  are  schematically  illustrated  in 
Fig.  1  for  zero  bias  and  for  a  voltage  F  applied  across  the 
contact  layers.  Band  bending  effects,  as  shown  in  Fig.  1  for 
nonzero  bias  voltage,  have  been  included  by  solving  for  the 
conduction-band-edge  energy,  which  acts  as  an  electrostatic 
potential  for  electrons,  using  Poisson’s  equation.  When  solv¬ 
ing  for  the  conduction-band-edge  energy,  we  do  not  include 
the  effects  of  charge  buildup  in  the  quantum  well  during 
device  operation,  i.e.,  we  do  not  solve  for  the  electron  wave 
function  and  the  position  of  the  conduction-band  edge  self- 
consistently .  We  assume  that  the  barrier  and  well  regions  are 
fully  depleted  of  carriers,  and  therefore  obtain  linear  voltage 
drops  across  these  layers.  This  assumption  is  included  only 
as  a  matter  of  convenience  and  is  not  necessary  for  achieving 
negative  differential  resistance.  We  find  that  including  band 
bending  effects  increases  the  peak  current,  since  the  bamer 
is  effectively  lowered  by  bending  in  the  accumulation  layer. 
However,  peak-to-valley  ratios  are  not  significantly  affected 
by  band  bending. 


1 

1 

1 

■>x 


FIG.  1.  Top:  Condoctioii-baiid  edge  of  the  double-barrier  jtructure  as  a 
function  of  position,  with  no  applied  voltage.  Bottom:  Conductioo-band 
edge  when  a  voltage  Kis  applied  to  the  device.  Note  dte  band  bending  in  the 
contact  layers. 


The  transmisson  coefficient  is  calculated  as  a  function 
of  energy  using  a  transfer  matrix  technique*  in  which  the 
device  is  divided  into  small  intervals,  and  the  potential  then 
assumed  to  be  constant  in  each  interval.  For  the  purposes  of 
calculating  the  transmission  coefficient,  we  consider  the  de¬ 
vice  to  include  the  well,  the  barriers,  and  electrode  regions  in 
which  band  bending  effects  are  not  negligible.  We  enforce 
the  usual  botmdary  conditions  between  intervals,  Le,  that  ^ 
and  (dij>/dx)/m*  be  continuous,  where  ^  is  the  envelope 
function  of  the  electron  wave  function  in  the  solid.  We  al^ 
assume  that  the  transverse  wave  vector  k|,  is  conserved  in  the 
tunneling  process,  and  that  the  total  energy  of  the  electron  is 
conserved.  These  two  assumptions,  combined  with  the  fact 
that  the  effective  masses  in  the  barriers  and  the  well  are  not 
equal,  lead  to  a  dependence  of  the  transmission  coefficient  on 
kg ,  but  we  found  this  effect  to  be  small 

For  electron  energies  above  the  conduction-band  edge, 
we  employ  the  effective  mass  approximation  £  =  1?k  V2m*. 
For  electron  energies  below  the  conduction-band  edge,  the 
wave  vector  is  imaginary;  its  value  is  calculated  using  Kane’s 
eight-band  model^  and  the  k*p  theory  outlined  by  Smith  and 
Mailhiot.*  To  obtain  the  transmission  coefficient,  we  assume 
plane  wave  solutions  for  outside  the  device.  Forx->  —  oo , 
we  have  incident  and  reflected  waves,  and  fora:—  -{-  oo ,  we 
have  a  transmitted  wave: 

x- —  oo, 

X--I-O0. 

The  transmission  coefficient  b  defined  to  be  the  ratio  of 
transmitted  to  incident  current,  T  =  Ar-lf  ( VAtf 

The  current  is  then  calculated  by  including  the  appro¬ 
priate  Fermi  factors  and  electron  velocities,  and  integrating 
over  the  incident  dectron  energy  distribution: 


^  =  ^(J  T(E,)AE){\-f(E  +  eV)]^-^d% 
-^T(E,)J(E+eV)\\-f(E)\^-l^d^k^  (1) 


em*kgT  f* 

Jo 


dEJ(E,) 


Xln 


l+ap[(Eg-E,)/kgT]  \ 
H-exp[(£f-£, -eV)/kgT]) 


(2) 


where  m*  is  the  effective  mass  in  the  electrodes,  is  the 
component  of  the  wave  vector  normal  to  the  bazzuis,  and  £j. 
=  fPkf/2m*.  InEq.  (1),  the  first  integral  conesponds  to 
electron  current  in  the  -|- x  direction,  and  the  second  inte¬ 
gral  to  current  in  the  —  X  direction.  An  integration  over  the 
'  transverse  wave  vector  k|  has  been  carried  out  to  obtain  (2). 
FromEq.  (1).  we  see  that  the  integrab  over  the  longitudinal 
wave  vector  both  contain  the  ffictor  idE/dk,^)dk^ 
=  d£^,  Le.,  the  effects  of  the  (Uffering  velocities  on  each  side 
of  the  device  are  canceled  by  the  effects  of  the  Afferent  one- 
dimensional  densities  of  states  at  a  given  energy.!^  cancel- 
'  lation  occurs  because  we  must  calculate  both  the  density  of 
states  and  the  electron  group  velocity  s  (dE/dk^  )/k  for 
the  side  of  the  device  on  which  the  current  is  incident,  be¬ 
cause  of  the  fact  that  we  have  defined  T(Ej^ )  as  the  ratio  of 
the  transmitted  to  incident  current  As  a  result,  the  expres- 
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sion  we  obtain  for  J  is  identical  to  that  obtained  by  Tsu  and 
Esaki,’ except  for  the  fact  that  we  have  defined  )  tobea 
transmission  probability  for  the  electron  current,  rather 
than  for  probability  densities. 

Figure  2  shows  current-voltage  curves  for  the  GaAs/ 
ZnSe  tunnel  structure  at  T  =  77  and  300  K.  The  peak  cur¬ 
rent  is  approximately  the  same  at  both  temperatures,  indi¬ 
cating  that  the  resonant  current  is  much  larger  than  the  ther¬ 
mionic  currents  over  the  barriers.  However,  the 
peak-to-valley  ratios  vary  considerably  with  temperature, 
since  the  barriers  formed  by  ZnSe  are  low  enough  to  allow 
substantial  thermionic  cunents  at  room  temperature.  At  77 
K  thermionic  currents  are  negligible,  and  we  have  calculated 
a  peak-to-valley  ratio  of  ~  1000;  at  300  K,  thermionic  cur¬ 
rents  become  important,  and  a  peak-to-valley  ratio  of  23  is 
obtained. 

Figure  3  shows  current-voltage  curves  for  InAs-ZnTe 
and  GaAs-AlGaAs  tunnel  structures.  The  dashed  curve  in 
Fig.  3  is  for  a  GaAs-AlGaAs  tunnel  structure  consisting  of 
undoped  40  A  AIq^jj  Gsq  47  As  barriers  and  an  undoped  50  A 
GaAs  well,  wirh  degenerately  doped  (10‘*  cm~^)  GaAs 
contact  layers.  The  peak  current  for  this  structure  is  some¬ 
what  higher  than  that  for  the  InAs-ZnTe  structure,  but  the 
peak-to-valley  ratio  is  much  lower  (92  vs  —  17(X3),  indicat¬ 
ing  that  InAs-ZnTe  devices  would  Lave  current-voltage 
characteristics  superior  to  those  of  GaAs-AlGaAs  struc¬ 
tures.  Furthermore,  the  large  peak-to-valley  ratio  obtained 
for  the  InAs-ZnTe  structure  would  vary  little  with  tempera¬ 
ture,  because  of  the  height  of  the  ZnTe  barrier.  Thus  InAs- 
ZnTe  tunnel  structures  would  appear  to  be  good  candidates 
for  devices  capable  of  operating  at  room  temperature. 

In  summary,  we  have  studied  theoreticaUy  GaAs-ZnSe 
and  InAs-ZnTe  double-barrier  resonant  tunneling  struc¬ 
tures.  These  structures  are  the  first  mixed  III-V/II-VI  tun¬ 
neling  devices  that  have  been  proposed  and  should  provide 
novel  ways  to  e.xamine  mixed  III-V/II-VI  materials  systems 


FIG.  2.  Current-voltage  characteristics  for  the  GaAs-ZnSe  tuniid  structure 
calculated  at  T  =  77  uid  300  K.  The  peak  currents  are  nearly  equal  for  the 
two  temperatures,  but  the  peak-to-valley  ratio  is  much  smaller  at  300  K  (23 
«  300  K  vs  - 1000  at  77  K). 
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FIG.  3.  Dashed  line:  Current-voltage  characteristic  for  a 
f*®As-AJQjj  Ga^t,  As  tunnel  structure^  Solid  line:  Current-voltage  charac¬ 
teristic  for  an  InAs-ZaTe  tunnel  structure.  Both  curves  are  calculated  at 
r=77K. 

and  growth  techniques.  Using  the  measured  values  for  the 
band  offset,  we  have  shown  that  the  GaAs-ZnSe  double-bar¬ 
rier  structure  should  produce  negative  differential  resis¬ 
tance.  We  also  note  that  the  smaller  dielectric  constant  m 
ZnSe,  compared  to  that  in  AlGaAs,  may  lead  to  somewhat 
faster  device  operation.  For  the  InAs-ZnTe  structure,  we 
calculate  a  conduction-band  offset  of  approximately  1.6  eV 
by  assuming  that  band  offsets  are  transitive  and  nqng  mea¬ 
sured  valence-band  offsets  for  GaAs-ZnSe,  ZaSe-ZnTe,  and 
InAs-GaAs.  We  obtain  very  large  peak-to-valley  ratios  for 
this  device,  and  the  size  of  the  band  offset  indicates  that  this 
device  should  work  well  at  room  temperature.  In  addition, 
the  small  dielectric  constant  in  ZnTe  and  the  high  mobility 
in  InAs  should  lead  to  a  substantial  advantage  in  speed  over 
GaAs-AlGaAs  devices,  and  the  high  mobility  in  InAs 
should  decrease  the  bias  voltage  required  to  achieve  the  max¬ 
imum  tunneling  current. 
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X-ray  photoelectron  spectroscopy  is  used  to  measure  the  valence-band  offset  in  situ  for  GaAs/ 
AlAs  ( 100)  heterojunctions  grown  by  molecular  beam  epitaxy.  Ga  3d  and  A1  2p  core  level  to 
valence-band  edge  binding  energy  differences  are  measured  in  GaAs  (100)  and  AlAs  (100) 
samples,  respectively,  and  the  A1 2/7  to  Ga  3d  core  level  binding  energy  difference  is  measured  in 
GaAs- AlAs  ( 100)  and  AlAs-GaAs  ( 100)  heterojunctions.  Measurements  of  the  A1 2p  to  Ga  3d 
core  level  energy  separations  indicate  that  the  band  offset  for  GaAs/AlAs  (100)  is  commutative; 
the  value  we  obtain  is  A£„  =0.46  ±0.07  eV.  Our  observation  of  commutativity  is  believed  to  be  a 
consequence  of  the  high  quality  of  our  GaAs/AlAs  (100)  heterojunctions,  and  of  the  inherent 
commutativity  of  the  GaAs/AlAs  ( 100)  band  offset. 


I.  INTRODUCTION 

Band  offsets  at  semiconductor  heterojunctions  are  of  great 
importance  in  solid-state  physics,  both  from  a  fundamental 
point  of  view  and  for  device  applications.  The  GaAs/AlAs 
band  offset  in  particalar  has  been  a  subject  of  e.xtensive 
study,  due  to  the  technological  importance  of  the  GaAs/ 
Al^Ga,  _j,As  heterojunction  system.  An  important  but  un¬ 
resolved  experimental  issue  in  this  system  is  that  of  commu¬ 
tativity  (i.e.,  independence  of  growth  sequence)  of  ±e  band 
offset.  Commutativity  is  consistent  with  current  theoretical 
models  and  some  experimental  results'-^;  however,  a  depen¬ 
dence  of  the  band  offkt  on  growth  sequence  has  also  been 
reported.’*  Published  experimental  values'"*^  for  the 
GaAs/AlAs  valence-band  offset  range  from  =  0.36  cV 
to  =  0.56  eV. 

In  this  paper  we  present  the  results  of  measurements  of  the 
band  offset  in  GaAs-.ALAs  (100)  and  AlAs-GaAs  (100) 
heterojunctions  using  x-ray  photoelectron  spectroscopy 
(XPS);  we  adopt  the  convention  of  listing  the  top  layer  in 
the  heterojunction  first,  e.g.,  GaAs-AlAs  indicates  a  hetero¬ 
junction  with  GaAs  grown  on  top  of  AlAs;  GaAs/AlAs  re¬ 
fers  collectively  to  both  growth  sequences.  We  describe  sam¬ 
ple  preparation  and  experimental  techniques  in  Sec.  IL  In 
Sec.  Ill  we  present  and  discuss  the  results  of  our  experi¬ 
ments.  Sec.  rv  concludes  the  paper. 

II.  EXPERIMENT 

The  samples  studied  in  this  experiment  were  grown  by 
molecular  beam  epitaxy  (MBE)  in  a  PHI  430  growth 
chamber  on  Si-doped  a*'-GaAs  { 100)  substrates.  The  sub¬ 
strate  preparation  consisted  of  a  degrease  followed  by  a  90-s 
etch  in  5:1:1  H-S04JIi02'H20.  The  native  oxide  was  de¬ 
sorbed  in  the  growth  chamber  by  heating  the  substrate  to 
~610’C  in  an  As  flux.  All  of  the  samples  were  grown  at 
600  *C  and  were  lightly  doped  n-type  with  .Si 
(/ic=lXl0'*cm~’)  to  avoid  both  sample  charging  effects 
and  excessive  band  bending  due  to  Fermi-level  pinning  at  the 
surface  of  the  specimen. 

Epitarial  GaAs  layers  >  1  pm  in  thickness  were  used  to 
make  the  bulk  GaAs  XPS  measurements.  Electron  mobili¬ 
ties  in  these  films  were  comparable  to  those  typically  report¬ 


ed  for  high-quality  MBE  growth,’  and  the  samples  exhibited 
good  photoluminescence  at  low  temperatures.  Two  bulk 
AlAs  samples  were  grown  for  this  study;  each  consisted  of  a 
2500-A-thick  epitaxial  AlAs  layer,  with  one  grown  on  top  of 
a  1.5-/im  GaAs  buffer  layer  and  the  other  on  top  of  a  1500-A 
M).3  G^o.?  As  buffer. 

Two  types  of  heterojunction  samples  were  grown,  allow¬ 
ing  the  commutativity  of  the  band  offset  to  be  checked.  One 
type  consisted  of  25  A  of  GaAs  on  top  of  either  500, 200,  or 
100  A  of  AlAs,  while  the  other  consisted  of  25  A  of  AlAs  on 
top  of  1(X)0  A  of  GaAs;  growth  interruption  was  not  used  in 
fabricating  the  heterojunction  samples.  In  both  cases,  the 
interfaces  were  beSeved  to  be  fairly  abrupt  based  upon  re¬ 
sults  from  quantum  weU  photoluminescence  (the  quantum 
well  luminescence  peak  width  was  3.8  meV  full  width  at  half¬ 
maximum  (FWHM)  at  9  /rW/cm^  incident  power)  and 
double-barrier  resonant  tunneling  experiments  on  samples 
with  similarly  grown  interfaces. 

XPS  measurements  were  obtained  using  a  Perkin-Elmer 
Model  5100  analysis  system;  samples  were  probed  using  a 
monochromatized  A1  Ka  x-ray  source  (Av=  1486.6  eV), 
and  electron  energies  were  measured  by  a  spherical  capaci¬ 
tor  analyzer.  Au  4f  core  level  peaks  had  full  widths  at  half¬ 
maximum  of  ~(X73  eV.  The  base  pressure  in  the  XPS 
chamber  is  typically  ~5Xl0“‘"  Torr,  and  the  analysis 
chamber  is  connected  to  the  MBE  chamber  via  a  UHV  trans¬ 
fer  tube,  allowing  samples  to  be  grown  and  analyzed  without 
being  exposed  to  atmosphere.  This  arrangement  eliminates 
experimental  uncertainties  associated  with  surface  passiva¬ 
tion,  sputter  cleaning,  and  annealing  of  samples  that  have 
been  transferred  through  atmosphere.  Measurements  were 
obtained  from  three  pure  GaAs  (100)  and  two  pure  AlAs 
( 100)  samples,  and  fiom  three  GaAs-AlAs  ( 100)  and  two 
AIAs-GaAs  ( 100)  heterojunctions.  Figure  1  shows  a  sche¬ 
matic  energy-band  diagram  for  the  GaAs/AlAs  interface. 
The  core  levels  ofinterest  in  the  band  offset  measurement  are 
the  A1  2/7  level  and  the  Ga  3d  level  in  AlAs  and  GaAs,  re¬ 
spectively.  From  Figure  1,  it  can  be  seen  that  the  valence- 
band  offset  is  given  by 

Aff,  =  (£:i{i:,  -  FgliS)  +  (f  8:5?  - 

(1) 
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Fig.  1.  Schematic  energy  band  diagram  for  the  GaAs/ AlAs  interface. 


The  core  level  to  valence-band  edge  binding  energy  differ¬ 
ences  are  obtained  from  measurements  on  pure  GaAs  and 
AlAs  samples,  and  the  A1  Ip  to  Ga  3d  core  level  binding 
energy  differences  are  obtained  from  the  heterojnnction 
samples. 

For  each  sample,  a  binding  energy  window  that  included 
aU  energy  levels  of  interest  was  scanned  repeatedly  for  ~  15- 
20  h.  Due  to  the  low  count  rate  for  electrons  near  the  val¬ 
ence-band  edge,  energies  in  that  region  were  sampled  for  a 
longer  time  than  energies  in  other  parts  of  the  spectrum. 
Spectra  from  pure  GaAs  and  pure  AlAs  and  from  GaAs- 


Mi 


Fio.  2.  Representative  bimSiig energy  XPS  spectra  for  (a)  GaAs  (tOO)  and 
(b)  AlAs  ( 100)  samples.  Cor  level  to  valence-band  edge  binding  energy 
differences  are  obtained  fiom  these  spectra.  The  discontinuities  in  the  tpee- 
tra  at  10  eV  binding  energy  aRdue  to  the  longer  sampling  times  used  daring 
data  acquisition  near  the  valence-band  edge;  the  valence4>and  spectra  aR 
also  shown  on  enlarged  intensity  scales,  as  indicated  in  the  8gure. 


Fig.  3.  Representative  binding  energy  X?S  spectra  for  (a)  CaAs-AlAs 
(100)  and  (b)  AlAs-GaAs  (100)  heterojunction  samples.  (2bR  level  to 
core  level  binding  energy  differences  aR  obtained  from  these  spectra. 


AlAs  and  AlAs-GaAs  heterojunctions  are  shown  in  Figs.  2 
and  3,  respectively.  In  Fig.  2,  the  discontinuities  in  the  spec¬ 
tra  at  10-eV  binding  energy  are  due  to  the  longer  sampling 
time  near  the  valence-band  edge;  the  valence-band  spectra 
are  also  shown  on  enlarged  scales,  as  indicated  in  the  figure. 

To  obtkin  accurate  core  level  to  core  level  and  core  level  to 
valence-band  edge  binding  energy  differences,  it  is  necessary 
to  determine  peak  positions  and  valence-band  edge  positions 
to  a  high  degree  of  precision.  To  find  core  level  p«k  posi¬ 
tions,  we  subtracted  from  each  core  level  peak  a  background 
function  (arising  from  inelastic  scattering  in  the  sample) 
proportional  to  the  integrated  photoelectron  intensity,  and 
defined  the  peak  energy  to  be  the  midpoint  of  the  two  ener¬ 
gies  at  which  the  intensity  was  one-half  the  maTimnm  inten¬ 
sity.  The  uncertainly  in  this  determination  of  the  peak  posi¬ 
tion  was  taken  to  be  ±  0.02  eV,  and  typically  core  level 
energy  separations  were  reproducible  to  better  than  0.01 
eV  for  a  given  sanqde: 

T o  locate  the  position  of  the  valence-band  edge  in  the  XPS 
spectrum,  we  employed  the  precision  analysis  technique  de¬ 
veloped  by  Kraut  etid.  In  this  approach,  the  XPS  spectrum 
near  the  valence-band  edge  is  modeled  as  a  convotution  of  a 
theoretical  valenceband  density  of  states  with  an  experi¬ 
mentally  determine  XPS  instrumental  resolution  function. 
This  model  function  is  then  fitted  to  the  experimental  data  to 
give  the  position  of  the  valence-band  edge.  Due  to  the  simi¬ 
larity  in  the  valence-band  spectra  for  GaAs  and  AlAs,^  we 
were  able  to  use  the  GaAs  nonlocal  pseudopotential  valence- 
band  density  of  states' '  in  analymg  data  from  both  materi¬ 
als.  For  a  given  sample,  core  level  to  valence-band  edge  bind¬ 
ing  energy  differences  were  reproducible  to  better  than 
±0.04eV. 

When  determining  the  position  of  a  valence-band  edge 
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xising  this  me±od,  complications  can  arise  from  the  presence 
of  surface  states  in  the  same  energy  region  as  the  valence 
band.’®  These  surface  sutes  produce  contributions  to  the 
XPS  spectrum  that  are  unaccounted  for  in  the  theoretical 
valence-band  density  of  states,  and  can  affect  the  calculated 
position  of  the  valence-band  edge  by  ~0.1  eV  or  more.  No 
such  effects  were  observed  in  our  AlAs  ( 100)  samples,  but  a 
Pj-Iike  surface  state  was  detected  in  the  GaAs  ( 100)  sam¬ 
ples,  in  agreement  with  previous  studies  of  the  GaAs  ( 100)- 
(2x4)  reconstructed  surface.'^  Thus,  the  position  of  the 
GaAs  valence-band  edge  was  determined  by  analyzing  the 
XPS  valence-band  spectrum  only  in  an  energy  region  in 
which  no  smface  states  were  observed. 

ill.  RESULTS  AND  DISCUSSION 

The  measured  core  level  to  valence-band  edge  binding  en¬ 
ergy  differences  are  shown  in  Table  I.  For  GaAs  ( 100),  the 
result  in  =  18.73  ±  0.03  eV,  and  for  AlAs 

(100)  we  obtain  (Zj^^  -  =  72.71  ±  0.04  eV. 

These  values  are  in  good  agreement  with  those  obtained  by 
Waldrop,  Grant,  and  Kraut."*  Table  II  lists  the  core  level 
binding  energy  differences  for  the  AlAs-GaAs  (1(X))  and 
GaAs-AlAs  ( 100)  heterojunctions.  For  the  GaAs-AlAs  in¬ 
terface,  we  find  {E  )  =  54.45  +  0.02  eV.  Using 

Eq.  ( 1 ) ,  we  see  that  this  value  corresponds  to  a  valence-band 
offset  of  0.47  ±  0.07  eV.  The  core  level  binding  energy  differ¬ 
ence  for  the  AlAs-Ga.\s  (100)  interface  was  found  to  be 
=  54.43  ±0.02  eV,  corresponding  to  a 
valence-band  offset  of  0.45  ±  0.07  eV. 

From  these  measurements  we  see  that  the  valence-band 
offsets  for  the  two  different  growth  sequences  are  within  0.02 
eV  of  each  other,  indicating  that  the  GaAs/AlAs  (100) 
band  offset  is  commutative.  The  observation  of  commutati¬ 
vity  for  GaAs/ AlAs  ( 100)  is  in  agreement  with  the  result  of 
Katnani  and  Bauer,'  but  disagrees  with  that  of  Waldrop  et 
al.*  who  reported  a  dependence  of  the  band  offset  on  both 
growth  sequence  and  crystal  orientation.  The  fact  that  XPS 
band  offset  measurements  are  not  consistently  commutative 
or  noncommutative  would  seem  to  indicate  that  growth  con¬ 
ditions  and  interface  quality  can  have  a  significant  effect  on 
the  value  of  the  band  ofikt.  If  this  were  the  case,  one  would 
expect  that  nearly  ideal  interfaces  would  exhibit  commutati¬ 
vity,  while  nonideal  interffices  would  not. 

The  value  obtained  for  the  band  offset,  0.46  eV,  is  well 
within  the  range  of  experimental  values  that  have  been  re¬ 
ported.  XPS  measurements  by  Waldrop  et  al*  yielded  a  va¬ 
lence-band  offset  of0.46eV  for  the  AlAs-GaAs  ( 100)  heter- 


Table  I.  Measured  core  level  to  valence-band  edge  binding  energy  (fiSer- 
ences  for  GaAs  ( 100)  and  AlAs  ( 100). 


Semiconductor 

Sample 

£22J-£?***(eV)  £:iJ;-5:‘**(eV) 

GaAs  (100) 

III024 

18.73 

GaAs  (100) 

111031 

18.72 

GaAs  (100) 

111054 

18.73 

AlAs (100) 

III027 

72.70 

AlAs (100) 

111038 

72.72 

Table  II.  Measured  A12pioCa3d  core  level  binding  energy  diderences  for 
GaAs-AlAs  (100)  and  AlAs-GaAs  (100)  heterojunctions. 


Heterojunction 

Sample 

£i!j;--^Si:(eV) 

GaAs-AlAs  ( 100) 

III036 

54.44 

GaAs-AlAs  (100) 

III083 

34.46 

GaAs-AlAs  (100) 

111087 

54.45 

AlAs-GaAs  (100) 

III045 

34.43 

AlAs-GaAs  (100) 

III050 

54.42 

ojunction,  in  agreement  with  our  result,  but  gave  a  value  of 
0.36  eV  for  the  GaAs-AlAs  (100)  band  offset.  The  latter 
result  suggests  that,  for  the  GaAs-AlAs  ( 100)  system,  dif¬ 
ferences  in  growth  conditions  between  our  samples  and 
those  of  Waldrop  et  oL  significantly  influenced  the  quality  of 
the  interface. 

Our  measurement  of  the  valence-band  offset  also  gives 
results  comparable  to  those  obtained  in  electrical  measure¬ 
ments.  Using  the  charge  transfer  technique,  Wang  and 
Stem’  obtained  a  valence-band  offset  of  0.45  ±  0.05  eV,  in 
excellent  agreement  with  our  results.  Batey  and  Wright* 
measured  a  valence-band  offset  of  0.55  eV  by  measuring  the 
activation  energy  for  thermionic  emission  of  holes  over  an 
AlAs  barrier;  this  value  is  somewhat  higher  than  ours,  but 
still  approximately  within  the  range  of  experimental  error. 

IV.  CONCLUSIONS 

We  have  measured  the  valence-band  offset  at  the  GaAs/ 
AlAs  ( 100)  interface  using  XPS.  Our  measurements  indi¬ 
cate  that  the  band  o£et  is  commutative,  and  the  value  we 
obtain  is  AF„  =  0.46  ±  0.07  eV.  This  value  is  in  good  agree¬ 
ment  with  other  XPS  band  offset  measurements,  as  well  as 
with  electrical  determinations  of  the  GaAs/ AlAs  (100)  ve- 
lence-band  offset.  Our  observation  of  commutativity  is  in 
agreement  with  some,  but  not  all,  published  results  for  the 
GaAs/ AlAs  ( 1(X))  interface.  This  indicates  that  we  have 
been  able  to  grow  high-quality  interfaces  between  GaAs  and 
AlAs,  since  we  expect  commutativity  in  ideal  interfaces,  and 
suggests  that  commutativity  of  the  band  offset  may  safely  be 
assumed  in  analyzing  and  designing  quantum  structnres  and 
devices. 
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Finite-size  effects  in  two-dimensional  continuum  percolation 
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We  have  investigated  the  finite-size  effects  of  anisotropic  contiauiun  percolation  in  two  dimen¬ 
sions.  The  elements  that  percolate  are  widthless  sticks.  We  have  developed  a  simple  theory  to  ex¬ 
plain  the  dependence  of  longitudinal  and  transverse  critical  leng;ths  on  anisotropy  and  the  finite 
number  of  sticks  in  the  sample.  By  comparing  the  theory  to  simulations,  «e  find  good  agreement. 


I.  INTRODUCTION 

In  recent  years  there  has  been  interest  in  the  electrical 
properties  of  composite  materials  consisting  of  conduct¬ 
ing  fibers  or  sticks,  embedded  in  an  insulating  polymer 
matrix.  Typically,  the  composites  are  prepared  by  mix¬ 
ing  the  fibers  and  polymer  in  a  liquid  state  and  letting  the 
result  cool  while  ^ing  poured.  This  leads  to  an  angular 
distribution  of  the  fibers  in  the  solid  state,  about  the  flow 
direction.  Experimentally,  such  composites  are  found  to 
have  a  threshold  dependence  of  the  electrical  conductivi¬ 
ty  on  the  fiber  length.'  Various  models  have  been  tried  to 
explain  such  behavior.  These  include  the  effective  medi¬ 
um  theory.^  This  has  not  been  very  successful  as  there  is 
a  large  difference  between  the  conductivities  of  the  fibers 
and  the  insulating  polymer.^ 

More  promising  are  percolation  simulations.  The  first 
was  done  by  Pike  and  Seager.*  They  considered  widthless 
(one-dimensional)  sticks  of  constant  length,  in  a  two- 
dimensional  medium,  isotropic  on  a  macroscopic  scale. 
Their  work  was  considerably  extended  by  Balberg  and 
Binenbaum’  to  the  case  of  anisotropic  systems  where,  as 
mentioned  above,  the  fibers  have  an  average  preferred 
orientation.  They  also  considered  various  distributions  of 
stick  length.  From  their  simulations,  the  critical  stick 
lengths  for  percolation  along  the  average  stick  orienta¬ 
tion  and  transverse  to  this  were  found  as  functions  of  an¬ 
isotropy  (defined  in  Sec.  II)  and  the  number  of  sticks  in 
the  sample.  However,  due  to  computing  constraints,  they 
only  considered  a  few  random  configurations  of  sticks. 
They  showed  that  in  the  limit  of  an  infinite  ensemble  of 
sticks,  the  longitudinal  and  transverse  percolation  thresh¬ 
olds  converge  to  a  common  function.  In  this  paper  we 
wish  to  extend  their  work  by  considering  more  thorough¬ 
ly  the  finite-size  effects  on  the  percolation  thresholds  of 
ensembles  with  small  numbers  of  sticks.  In  Ref.  5.  and 
other  papers*’ on  continuum  percolation  this  appears  to 
have  been  a  neglected  point  of  interest.  We  will  derive  a 
simple  model  for  the  dependences  of  the  thresholds  on 
anisotropy  and  numbers  of  sticks  and  compare  this  to  the 
results  of  extensive  simulations.  *  -  ^2: 

Our  paper  is  arranged  as  follows.  Section  n  describes 
the  simulations.  Section  III  contains  the  model  of  the 
percolation  thresholds.  Section  IV  compares  the  model 
with  the  simulations.  FinaUy,  we  give  a  conclusion  in 
Sec.  V. 


n.  SIMULATIONS 

In  this  section  we  outline  a  method  of  performing 
simulations  to  find  the  critical  lengths.  (This  is  based  on 
the  procedure  by  Balberg  and  Binenbaum.^)  Hrst,  we 
define  the  basic  terms.  Then,  we  give  the  procedure  for 
obtaining  the  critical  lengths. 

To  start,  let  us  define  some  quantities.  Consider  a  set 
of  M  widthless  sticks.  We  place  the  centers  of  the  sticks 
uniformly  in  the  unit  ^uare,  [0,1]  X  [0,1].  The  sticks 
have  some  angular  probability  distribution  /(0)  about 
the  y  axis,  where  9  is  the  angle  between  a  stick  and  the  y 
direction.  The  angular  distributions  we  consider  all  have 
the  feature  that 

/(9)=/(-9),  9e(0,ir/2)  (1) 

and  that  the  average  angle  is  £[9]=0.  The  sticks’ 
lengths  are  given  by  another  probability  distribution, 
g(L).  We  assume  $  and  L  to  be  independent.  Given  a 
random  configuration  of  fif  sticks  with  lengths  and  angles 
((£.{,  9; )] ,  the  following  quantities  may  be  defined: 

2  A  Iws9,  |  (2a) 

,--i 

and 

=  .  (2b) 

"  i-l 

These  are  the  average  longitudinal  and  transverse  stick 
components  with  respect  to  the  y  direction,  re^)ectiveiy. 
From  these  we  define  the  macroscopic  anisotropy  as 

P, 

For  an  isotropic  sample,  we  have  P  1.  The  larg^  P  is, 
the  more  oriented  the  sticks  are  along  the  7  direction.  A 
given  sample  of  sticks  is  considered  to  be  percolating 
along  the  y  direction  if  a  continuous  path  can  be  traced 
between  intersecting  sticks  from  y  =0  to  7  ==  1.  SmOarly 
for  the  X  direction.  The  longitudinal  critical  length, 
is  the  lowest  average  len^  that  gives  the  onset  of  y  per¬ 
colation.  The  transverse  critical  length,  is  likewise 
defined  for  x  percolation.  .  ~  . 

We  wish  to  investigate  the  dependences  of  and 
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on  the  anisotropy,  P,  and  on  the  number  of  sticks  in  the 
system,  N.  First,  consider  systems  with  a  given  iV.  To 
vary  the  anisotropy,  we  change  the  variance  of  /(5).  For 
example,  for  a  normal  distribution  of  angles  /  is  iV(0,d^), 
where  6^  is  the  standard  deviation.  By  progressively  re¬ 
ducing  6^  from  some  initial  value,  we  can  corresponding¬ 
ly  increase  the  anisotropy.  Consider  now  that  we  have 
selected  a  variance  of  /(0).  We  start  with  some  small 
average  value  of  £,  E\L].  With  the  parameters  of/ and 
g  chosen,  20  random  configurations  of  sticks  are  generat¬ 
ed,  each  with  N  sticks.  We  look  for  any  percolation. 
Twenty  random  configurations  at  each  value  of  anisotro¬ 
py  and  E[L]  were  made  to  increase  the  statistical  relia¬ 
bility  of  finding  accurate  critical  lengths.  This  is  an  im¬ 
provement  over  Balberg  and  Binenbaum  who  presented 
critical  lengths  found  from  three  to  five  configurations 
per  value  of  anisotropy.  Note  that  for  sufficiently  small 
EIL],  few  of  the  sticks  will  overlap.  Systematically, 
£[Lj  is  incremented  until  percolation  is  foimd.  Often 
this  is  longitudinal  percolation.  This  should  be  e.xpected, 
as  it  is  easier  to  percolate  along  a  preferred  direction  than 
normal  to  it.  (Occasionally,  at  low  anisotropy.  P  £  1.4, 
we  may  first  encounter  transverse  percolation.)  By  fur¬ 
ther  incrementing  £[L],  we  eventu^y  come  upon  trans¬ 
verse  percolation.  During  the  varying  of  £[L]  we  record 
P  for  each  random  configuration.  It  is  found  that  to  good 
approximation,  P  is  independent  of  £[L].  Our  results 
and  a  comparison  with  the  theory  of  Sec.  UI  will  be 
presented  in  Sec.  FV. 


FIG.  I.  A  stick  AC  of  length  L  is  at  an  angle  d*  with  respect 
to  the  y  axis,  where  ff  is  defined  in  Eq.  (3).  Sticks  at  an  angle 
— y,  with  centers  inside  the  parallelogram  HFIG  wfll  intersect 
AC. 


m.  THEORY 


Here  we  present  a  simple  model  that  attempts  to  ex¬ 
plain  the  dependence  of  the  longitudinal  and  transverse 
critical  lengths  on  the  anisotropy  and  the  number  of 
sticks  in  the  ensemble.  As  before,  let  there  be  N  sticks  in 
the  ensemble,  with  an  angular  distribution  fid)  and  a 
stick  length  distribution g{L).  Our  first  approximation  is 
to  replace  all  the  stick  lengths  L,  with  the  average  length 
LQ~'2,iLi/N.  Next,  let  us  replace  the  angular  distribu¬ 
tion  /(0)  by  one  in  which  half  the  sticks  are  oriented  at 
y  and  the  other  half  at  —  0'  with  respect  to  the  y  axis, 
where  6'  is  defined  by 

•-  <|cos0l>  P,  --*1- 

For  definitiveness,  we  now  look  at  longitudinal  percola¬ 
tion.  That  is,  we  wish  to  find  a  set  of  overlapping  sticks 
that  goes  from  y^O  to  y  =  1.  Consider  a  stick 
angle  O',  labeled  AC  in  Fig.  I,  where  £  is  the  cateTcS 
the  stick.  The  only  sticks  which  can  intersect  AC  whh 
nonxero  probability  are  oriented  at  —ff.  Furthermore 
the  centers  of  these  sticks  must  lie  in  the  parallelognGn 
shown  in  Fig.  1.  (Balberg  et  al.*  define  this  as  the  exclud¬ 
ed  area  of  the  two  sticks.)  Suppose,  in  tracing  out  a  p^ 
sible  cluster,  that  we  started  at  y  =0  and  that  AC  is  the' 
highest  stick  in  the  cluster,  thus  far.  Let  AC  be  the}th 
stick  in  the  cluster.  To  make  progress  towards  y=l  we 
ask  for  the  probability  that  a  stick  has  a  center  in  EHFG 
and  is  oriented  at  —ff.  Remembering  that  the  sticks  are 
distributed  uniformly,  the  probability  is  given  by  the 


product  of  the  probability  of  one  stick  intersecting  AC 
with  a  higher  y  center  times  the  number  of  sticks. 
Thus  we  have 


E-l  =pp, 


Pj^^Ai^FG)  =PA  J 


where  A  is  the  area.  [Since  a  probability,  if  the 

right-hand  side  of  Eq.  (4)  is  greater  thaw  i,  we  set 
.1  Furthermore,  if  such  an  intersectioa  occurs  we 
see  that  the  average  y  coordinate  of  the  intersecting  stick 
will  be  at  a  distance  F//3  higher  than  the  y  coonfinate  of 
ACs  center.  Thns^  the  average  number  of  intersections 
for  y  percolation  is 


The  probability  of  finding  one  percolating  cluster  in  the  y 
direction  is  given  by  the  product  of  n  terms  of  the  form  of 
Eq.  (4),  whichwecanwriteas  r--' ' 

Pd««:=n^. 


where  from  Eq.  (4),  Py  is  defin^  as '  . 
jK  a/<l  .-  V. 

^^=[1,  aj>l 
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aj=P.P,  .  (8) 

Note  that  in  Eq.  (6)  we  have  assumed  that  n  is  even.  If  n 
is  odd,  Eq.  (6)  is  multiplied  by  •P «/2+l- 
Equation  (6)  gives  the  probability  of  forming  one  clus¬ 
ter.  However,  it  may  be  possible  for  several  clusters  to 
form,  especially  for  longitudinal  percolation  in  highly  an¬ 
isotropic  samples.  As  mentioned  above,  the  y  coordinate 
of  the  intersecting  stick  in  the  cluster  is  at  an  average  dis¬ 
tance  of  Pi  /3  higher.  For  the  x  coordinate  of  the  inter¬ 
secting  stick,  averaging  over  ^FG  gives  0.  This  comes 
from 

Fix  <jcj)=P(x>Xj)=L  , 

where  F  is  the  probability  function.  But,  we  see  that  the 
average  x  coordinate  of  the  interesting  stick  is  P,/3,  as¬ 
suming  that  x  >Xg,  and  it  is  — P,/3  if  x  <Xg.  Hence, 
for  longitudinal  percolation  we  have  a  random  walk  in 
the  X  direction,  of  average  step  size  P,  /3  and  equal  prob¬ 
ability  of  stepping  positive  or  negative.  From  Eq.  (3),  it 
takes  n  steps  to  percolate  in  the  y  direction.  Daring  these 
steps,  the  cluster  win  span  an  approximate  distance’  in 
the  X  direction  of 


Since  we  are  bounded  by  0  and  1  in  the  x  direction,  let 
Uq.  oo-l 

“'ll,  .,>!■ 

Then  v  is  the  hractional  transverse  distance  covered  when 
we  are  looking  at  longitudinal  percolation.  Therefore, 
the  total  probability  for  a  longitudinal  percolating  cluster 
is  given  by  dividing  Eq.  (6)  by  v  to  give 


To  find  transverse  percolation,  we  merely  interchange  P; 
and  P,  in  the  above  equations.  Given  values  for  ( |  $in6| ) 
and  ( I  cosdl ),  we  can  use  the  above  equations  to  find  the 
critical  lengths  and  the  anisotropy.  ...i  . 

.  One  further  matter  needs  to  be  addressed.  Given  a 
{ 1  cosdl ),  which  b  defined  by  a  particular  /(fi)  and  a 
choice  of  variance  of  /,  there  correspond  a  unique 
<  I  sindl ).  But  for  the  set  of  arbitrary  probabifity  distri¬ 
butions  /,  subject  to  the  constraints  of  Eq.  (1),  there  will 
be  no  unique  relation  between _( | cosdl )  .ud  {|sin0(). 
How  can  we  find  an  anisotropy  that  is  in  some  sense  in¬ 
dependent  of  a  particnlar  /?  To  answer  this,  note  that  we 
have  the  following  ineq»salities:  ;  • . .  - 


•“  ^0<<lcosB|>,  {|sinB|><l  -  --'J 

;c 

and 

l<d<2,  ' 

a2b) 

where  .  t  ■ 

,  d  =  (lcosBl>+(lsinBl)  , 

-.(12c') 

and  the  triangle  inequality  was  used  to  get  the  lower 
bound  of  Eq.  (12b).  We  do  the  following.  A  value  of  d  is 
chosen  to  satisfy  Eq.  (12b).  We  then  vary  (|sind|)  in 
steps,  starting  from  d  /2  and  decrementing  to  0  (though 
not  going  all  the  way).  Then  <  1  cosd| )  is  given  by 


<  I  sindl  >  = 


d  — <|sind|>,  d  — <|sind|)<l 
1,  d  — (lsind|)>l. 


(13) 


With  these  values,  we  can  find  the  anisotropy  and  use  Eq. 
(11)  to  obtain  the  critical  lengths.  We  choose  U,  the 
midpoint  of  Eq.  (12b),  as  the  most  reasonable  value  of  d. 
It  should  be  realized  that  the  procedure  of  ming  Fqs 
(12c)  and  (13)  is  an  approximation.  In  general,  a  given  / 
will  not  yield  { |  cosdf )  and  { |  sindj )  satisfying  the  linear 
relationship  of  Eq.  (12c)  over  a  range  of  values  of  the 
variance. 

A  program  was  written  to  find  the  lowest  lengths  that 
set  Pto,  -  in  Eq.  (11)  as  a  function  of  anisotropy,  for  both 
longitudinal  and  transverse  percolation.  These  are  the 
critical  lengths.  The  number  of  sticks  in  the  sample  was 
a  parameter  in  this  calculation.,.  > 


IV.  RESULTS 


Here  we  present  and  compare  the  results  of  Sec.  n  and 
in.  First,  let  us  consider  the  predictions  of  the  model  of 
the  preceding  section.  The  values  of  critical  lengths  from 
Eq.  (11)  are  displayed  in  Figs.  2  and  3  for  1(X)  and  500 
sticks,  respectively.  In  each  figure,  the  solid  line  is  for 
longitudinal  percolation  while  the  dashed  line  is  for 
transverse  percolation.  Following  Balberg  and  Binen- 
baum,’  we  have  normalized  the  critical  lengths  in  units  of 
the  average  interstick  separation,  r,  where 


(14) 


From  the  simnlations  we  have  found  the  critical 
lengths  Le|,  as  functions  of  anisotropy,  for  ensembles 
of  100  and  300  sticks.  These  are  displayed  in  Rgs.  2  and 
3,  respectively.  For  Fig.  2,  we  have  found  the  critical 
lengths  for  three  types  of  ensembles:  Those  with  normal 
distributions  for  /and  g;  with  a  uniform  f  and  a  normal 
g;  and  with  constant  absolute  angle  and  a  S  function  for  g 
(i.e.,  all  the  sticks  in  an  ensemble  are  the  same  length). 
This  last  pair  of  distributions  corresponds  to  the 
simplified  choices  of  distributions  made  in  Sec.  IIL  Simi¬ 
lar  remarks  hold  for  Fig.  3.  The  solid  objects  in  the 
figures  are  the  critical  longitudinal  lengths.  The  hollow 
objects  are  the  critical  transverse  lengths. 

We  see  from  Figs.  2  and  3  that  for  both  the  simulations 
and  the  theory  the  transverse  critical  lengths  fie  distinctly 
above  the  longitn£nal  critical  lengths.  Consdering  sepa¬ 
rately  the  transverse  and  longitudinal  results,  we  see  that 
the  simulations  with  constant  length  and  sbsdnte  angle 
tend  to  yield  larger  critical  lengths  than  the  ether  timula- 
tions.  This  is  eapected,’  as  a  distribution  lengths  will 
cause  the  sticks  with  lengths  greater  *bwn  the  mean 
length  to  contiibate  preferentially  to  the  percolation. 
Hence  percolaticm  will  start  sooner  than  if  all  the  lengths 
in  a  sample  are  constant.  Nonetheless,  the  simulation  re- 
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FIG.  2.  The  critical  lengths  as  a  function  of  anisotropy  for 
samples  of  100  sticks.  The  anisotropy  is  defined  in  Eq.  (2).  The 
solid  and  open  symbols  are  from  simulations  for  longitudtnai 
and  transverse  percolation,  respectively.  The  lines  are  from  Eq. 
(11). 


FIG.  3.  The  critical  lengths  as  a  function  of  anisotropy  for 
samples  of  SOO  sticks.  The  anisotropy  is  defined  in  Eq.  (2).  The 
solid  and  open  symbols  are  from  simulations  for  longitudinal 
and  transverse  perctdation,  respectively.  The  lines  are  from  Eq. 
Ul).  ,  ' 


suits  for  :  .transverse  percolation  are  clustered  close 
enough,  and  likewise  for  longitudinal  percolation,  that  we 
find  a  universal  behavior  of  the  critical  lengths  on  the  an¬ 
isotropy  that  is  largely  independent  of  the  choice  of  dis¬ 
tributions.  Comparing  the  theory  with  the  simulations, 
we  see  good  agreement  for  both  1(X)  and  5(X)  sticks.  'The 
largest  disagreement  is  for  the  longitudinal  percolation  of 
1(X)  sticks,  with  the  theory  lying  above  the  smnlations. 
But  even  here^  the  theoretical  curve  exhibits  the  s^e 
trends  as  the  simnlations  and  the  disagreement  is  <^y 
semiquantitative.  tS-  :  . 

.^(Note  that  in  comparing  Figs.  2  and  3,  the  longitm&^ 
and  transverse  results  tend  to  converge  together  as  m 'go 
from  100  sticks  to  300  sticks.  This  is  in  accordance  ^ 
Balberg  and  Binenbanm^  who  showed  by  a  topolopcal 
argument  that  in  the  fimit  of  infinite  the  two  types  of  ^ 
critical  lengths  coindde.  r 

^';Let  us  add  a  comment  about  the  midrange  choice  of,_ 
d  =  1.5  made  in  Eqs.  (12)  and  (13)  in  the  prece«fing1w^5 
tion.  This  was  done  to  give  a  definite  prescription  for  Ae 
anisotropy.  If  we  set  d  equal  to  the  lower  limit,  1,  we  g^ 
curves  that  start  at  a  critical  length  mund  S  at  notropyT 
and  remain  above  the  displayed  curves  as  we  increase  the 
anisotropy.  This  was  found  for  both  longitndinal  and 


transverse  percolation.  The  agreement  between  the  (un¬ 
displayed)  d=  1  cmves  and  the  simulations  is  not  as  good 
as  for  the  d  =  l.S  carves  shown  here  in  Figs.  2  and  3.  Al¬ 
ternatively,  choosing  d  equal  to  the  upper  Smh  of  2  gives 
curves  that  start  and  remain  slightly  lower  than  the 
displayed  curves.  Hencej,'  we  believe  bur  choice  of  d  =  1. 5 
to  be  a  reasonable  one. '  : 

"  Therefore,  we  suggest  that  our  theory  erf*  Sec.  HI  gives 
good  predictions  of  critical  length.  The  theory  also  has 
the' computational  advantage  of  being  much  &ster  to  Iron. 
.'The  curves  in  Figs.  2  and  3  took  less  than  a  minute  of 
CPU  (central-processing-iinit)  time  each  to’ find,  on  a 
^'VAX  1 1/785.  By  contrast,  the  CPU  time  for' die'smula- 
'Tions  is  measored  in  hours.  This  is  especially 'hue  for  the 
^simulations  of  5(X)  sticks.  The  CPU  time  teqiined  scales 
roughly  as  the  sqnare  of  the  number  of  sticky  as  the  most 
'  intensive  task  a  to  find  the  possible  intersections  amongst 
^fsnicki.  - 

•^1:7  :  V.  coNcusiON  ' 

We  have  investigated  the  finite-size  efifects  of  anisotrop¬ 
ic  continuum  percolation  in  two  dimensions.  The  ele¬ 
ments  that  peredate  are  widthless  sticks.'  We  have 
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developed  a  simple  theory  to  explain  the  dependence  of 
longitudinal  and  transverse  critical  lengths  on  anisotropy 
and  the  finite  number  of  sticks  in  the  saimple.  By  com¬ 
paring  the  theory  to  emulations,  we  find  good  agreement. 
It  also  affords  significant  computational  advantages  over 
performing  simulations.  We  believe  that  our  theory  is 
the  first  to  explain  the  finite-size  simulation  results  for  an¬ 
isotropic  continuum  percolation  in  two  dimensions. 


Thus,  we  suggest  that  the  theory  offers  a  useful  comple¬ 
ment  to  the  running  of  simulations. 
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The  tunneling  time  for  electrons  to  escape  from  the  lowest  quasibound  state  in  the  quantum 
wells  of  GaAs/AlAs/GaAs/AlAs/GaAs  double-barrier  hetcrostructures  with  barriers 
between  16  and  62  A  has  been  measured  at  80  K  using  photoluminescence  excitation 
correlation  spectroscopy.  The  decay  time  for  samples  with  barrier  thicknesses  from  16  A  ( =:  12 
ps)  to  34  A(  =800  ps)  depends  exponentially  on  barrier  thickness,  in  good  agreement  with 
calculations  of  electron  tunneling  time  derived  from  the  energy  width  of  the  resonance. 

Electron  and  heavy  hole  carrier  densities  are  observed  to  decay  at  the  same  rate,  indicating  a 
coupling  between  the  two  decay  processes. 


The  electrical  properties  of  the  double-barrier  (DB) 
heterostructure  have  been  of  great  interest  since  its  proposal 
by  Tsu  and  Esaki.'  The  desire  to  characterize  the  high-fre¬ 
quency  behavior  of  the  DB  stems  from  interest  in  its  use  as  an 
oscillator'  and  as  a  switching  element.^'^  However,  the  time 
associated  with  the  tunneling  of  electrons  has  been  the  sub¬ 
ject  of  more  than  20  years  of  discussion.^’  Experimental 
measurements  of  tunneling  times  have  required  the  develop¬ 
ment  of  high-speed  measurement  techniques.  Recently,  two 
groups  have  reported  experimental  studies  of  the  temporal 
response  of  DB  structures.  Whitaker  et  al.'^  have  used  elec¬ 
tro-optic  sampling  measurements  to  study  a  single  tunnel 
device.  Tsuchiya  et  c/."  used  the  photoluminescence  (PL) 
from  carriers  in  the  quasibound  states  in  the  quantum  weO  to 
study  the  decay  of  the  electron  population  in  the  quantum 
well  as  a  function  of  the  barrier  thickness.  However,  the  ex¬ 
periments  of  Ref.  1 1  were  limited  to  times  greater  than  60  ps, 
and  Guo  et  al.^  have  raised  some  serious  questions  about  the 
theoretical  calculatious  that  were  used  for  comparison  with 
experimental  results. 

In  this  letter  we  report  a  study  of  the  decay  of  photocx- 
cited  carriers  in  double-barrier  hetcrostructures  as  a  func¬ 
tion  of  the  thickness  of  the  barrier  layers.  By  employing  an 
e.xcitation  correlation  method '  ^  we  have  extended  the  results 
of  Tsuchiya  et  a/."  to  significantly  shorter  times  and  are 
comparing  our  data  with  results  of  corrected  calculations  of 
the  tunneling  times. 

Double-barrier  (DB)  structures  were  grown  on 
(100)  GaAs  substrates  by  molecular  beam  epitaxy  in  a  Per- 
kin-Elmer  430  system  at  600  *C  After  grow^  of  0 J  iita  of 
GaAs,  a  superlattice  buffer  layer  consisting  of  five  periods  of 
50  A  A1<,.j5  Gao.6j  As/500  A  GaAs  was  grown.  This  was  fol¬ 
lowed  by  growA  of  a  0.7  (im  layer  of  GaAs,  which  provided 
a  high  quality  layer  on  which  to  grow  the  DB  and  eliminated 
any  optical  effects  from  the  superlattice.  Then  a  symmetrical 
GaAs/AIAs/GaAs/AiAs/GaAs  DB  was  grown,  with  a 
well  thickness  of  58  A.  The  final  layer  was  a  300  A  GaAs  cap. 
All  layers  were  nominally  undoped  with  an  estimated  resid¬ 
ual  carbon  acceptor  concentration  of  10'^  cm  Seven  sam¬ 
ples  were  studied,  with  bulk  growth  rate  information  pre¬ 


dicting  barrier  thicknesses  of  16, 22, 28, 34, 34, 48,  and  62  A. 
High-resolution  transmission  electron  microscopy  con¬ 
firmed  the  barrier  thicknesses  of  the  16  A  sample  and  one  of 
the  34  A  samples,  within  an  uncertainty  of  two  monolayers. 
We  estimate  an  uncertainty  in  barrier  thickness  of  two  mon¬ 
olayers  for  all  of  the  samples. 

The  experimental  technique  has  been  desenbed  pre¬ 
viously.'^  A  colliding  pulse  mode-locked  ring  dye  lastf  is 
used  to  generate  a  train  of  pulses  2(X)  fs  full  width  at  half 
maximum  (FWHM),  at  a  repetition  frequency  of  120  MHz. 
The  laser  output  is  centered  at  6200  A  and  has  a  spectral 
width  of  20  A  FWHM.  The  pulse  train  is  equally  divided 
into  two  separate  beams  which  are  independently  chopped  at 
/,  =  1600  Hz  and/-  =  2100  Hz,  and  d^yed  with  respect  to 
one  another  by  time  y  (  —  5(X)<y<500  ps)  before  b^g  re¬ 
combined  and  focused  to  a  2S-jum-diam  spot  on  the  surface 
of  the  sample.  The  typical  average  power  used  was  1  mW  per 
beam  before  chopping.  The  PL  is  spectrally  resolved,  and 
detected  by  a  GaAs  photomultiplier  tube  (PMT).  After  am¬ 
plification,  the  PMT  signal  is  synchronously  detected  by  a 
lock-in  amplifier  at  ehber  the  fundamental  frequency /,  or  at 
the  sum  frequency/,,,  =/,  -b/.  All  of  the  re^ts  reported 
here  were  taken  with  the  sample  at  80  K. 

In  Fig.  1  we  present  a  schematic  diagram  of  the  pro- 


FIG.  1.  Schematic  diagnon  of  relevant  carrier  processes  involved  in  the 
double-barrier  samples  dnring  the  experiment  Shown  are  pbotoexdtation 
of  carriers  in  the  nmneiing  of  caniers  out  of  the  well,  and  recombina¬ 

tion  of  carriers  in  the  wdL 
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cesses  of  excitation,  tunneling  from  the  well,  and  the  radia¬ 
tive  recombination  of  carriers  within  the  well.  Recent  obser¬ 
vations  of  the  thermalization  times  of  electrons  between 
subbands  have  given  times  less  than  200  fs.'^  Hence  we  will 
assume  that  the  thennalization  of  electrons  and  holes  to  the 
lowest  subband  to  be  fast  compared  to  the  times  of  interest 
here.  A  model  of  the  measurement  process  has  been  given  in 
Ref.  12.  We  will  simply  sketch  the  important  results.  If  the 
electron  and  hole  populations  created  by  an  optical  pulse  are 
n(t)  and  pit),  and  the  populations  created  by  two  optical 
pulses  are  independent,  then  the  sum  frequency  signal/^  is 
proportional  to  the  cross  correlation 

Aun,{r)«  J  [fiinpit  -  r)  +  nit  -  r)pit)]dt. 

This  signal  is  due  to  the  recombination  of  electrons  created 
by  the  first  pulse  with  holes  created  by  the  second  pulse,  and 
vice  versa.  If  we  assume  that  the  electron  and  heavy  hole 
densities  decay  exponentially^  with  time  constants  r,  and 

respectively,  then  the  sum  frequency  signal  is  propor¬ 
tional  to  the  sum  of  two  exponentials 

Aum(/)  “  [exp(  -  |y|/rj  -fexp(  -  \y\/'v>,)\-  (D 
We  have  used  this  reult  as  the  basis  for  interpreting  our 
data. 

In  Fig.  2  we  present  typical  photolumincscence  spectra 
taken  at  80  K  at  the  fundamental  and  sum  chopping  frequen¬ 
cies.  The  spectrum  at  the  fundamental  frequency  consists  of 
a  single  distinguishable  feature  centered  at  7650  A-  The 
wavelength  of  the  feature  in  the  fundamental  frequency 
spectrum  is  in  approximate  agreement  with  the  calculated 
position  of  7730  A  for  transition  from  the  lowest  electron 
subband  to  the  lowest  heavy  hole  subband  in  aS8  A  quantum 
well.  The  peak  of  the  corresponding  feature  in  the  sum  fre¬ 
quency  spectrum  is  shifted  slightly  to  longer  wavelengths. 
We  do  not,  at  present,  have  an  explanation  for  this  shift 

In  Fig.  3  we  present  a  semilogarithmic  plot  of  a  typical 
scan  of  the  photolominescence  at  the  sum  chopping  frequen¬ 
cy  as  a  function  of  the  time  delay  between  the  two  pulses.  The 
scan  was  taken  at  80  S,  at  a  wavelength  of 7665  A,  the  peak 
of  the  sum  PL  spectrum  shown  in  Fg.  2.  This  delay  scan 
consists  of  a  peak  at  zero  delay  with  wings  extending  to  much 
longer  times.  The  peak  at  y  =  0  is  a  coherence  peak  due  to 


Wavelength  (A) 


fig.  2.  Typical  exciiatiaa  correlated  Itiniinesceiice  for  2S  A  barrier  junple 
at  80  K.  Shown  are  haameacence  signals  at  the  Audamenol  chopping  fre¬ 
quency  (solid  tine)  and  the  sum  frequency  (dashed  line).  Both  teas  were 
taken  with  delay  y  w  0.  The  sum  frequency  spectrum  has  been  muldptied  by 
2. 
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FIG.  3.  Semilogaritheac  plot  of  the  variation  of  the  sum  fiequency  lumines. 
cence  signal  (solid  Sae)  with  delay  y.  The  sample  is  the  same  as  in  Fig.  2. 
The  scan  was  taken  at  766S  A,  the  peak  of  the  sum  frequency  spectrum.  The 
coherence  peak  at  y  =  0  is  not  resolved  in  this  scan.  The  riMhwt  lines  are  the 
fits  discuss^  in  the  xn,  which  give  a  decay  time  of 236  ^  20  ps. 

the  Optical  intcticrence  of  the  two  incident  pulses  on  the 
sample  and  is  net  resolved  in  this  scan.  A  single  exponential 
fits  the  wings  in  the  sum  frequency  delay  scan  shown  in  Fig. 
3.  Delay  scans  for  one  sample,  taken  at  4  K,  show  evidence  of 
a  second  exponential  with  a  faster  decay.  Fits  to  the  sum 
frequency  delay  data  at  80  K  using  a  single  exponential  are 
shown  as  dashed  lines  in  Fig.  3.  The  coherence  peak  is  not 
included  in  fitting  the  delay  scans. 

In  Fig.  4  we  have  plotted  the  exponential  decay  time  at 
80  K  as  a  function  of  barrier  thickness  for  the  seven  samples 
described  above.  The  decay  time  depends  exponentially  on 
barrier  thickness  for  barriers  up  to  approximately  34  A. 
Over  this  range  of  exponential  dependence  the  decay  time 
varies  by  two  orders  of  magnitude.  For  thicker  barriers  the 
decay  time  seems  to  approach  a  value  that  is  independent  of 
the  barrier  thickness.  With  the  120  MHz  repetition  rate  of 
the  excitation  pulses,  there  is  an  upper  limit  on  the  order  of  2 
ns  to  the  decay  times  that  can  be  measured  accurately  with 
our  technique.  Consequently,  the  result  for  the  sample  with  a 
barrier  thickness  of  62  A  should  be  viewed  with  some  cau- 


Barrier  Thickness  (A) 


FIG.  4.  Measured  deay  times  as  a  function  of  harrier  tfricknas.  The  data 
points  are  the  mesmed  decay  times,  with  error  bars  on  the  thickness  based 
on  uncertainty  in  the  barrier  thickness  and  error  baum  on  the  decay  times 
from  uncertaintie  in  (he  fits  to  the  dday  scans.  The  solid  fine  is  the  dectron 
tunneling  time  cairnbteri  from  the  width  of  the  lowest  coodaction.band 
transmission  resonance,  using  a  tsn>4iaad  expressioa  for  the  wave  vector  in 
AlAs  harriers.  Thedsbed  line  shows  the  tunneling  time  for  heavy  boles  and 
was  calculated  usng  a  one.band  expression  for  the  barrier  wave  vector. 
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tion.  Our  measured  decay  times  are  about  a  factor  of  4  longer 
than  those  of  Tsuchiya  el  al. ' ' 

In  Fig.  4  we  have  also  plotted  the  time  calculated  for 
electrons  and  heavy  holes  to  tunnel  out  of  the  quantum  well. 
The  solid  line  is  the  electron  tunneling  time,  and  the  dashed 
line  is  the  heavy  hole  tunneling  time.  We  will  assume  that  the 
light  holes  thermalize  to  become  heavy  holes  on  a  time  scale 
short  compared  to  times  of  interest  here.  From  Ref.  7  the 
time  for  a  particle  to  tunnel  out  of  a  quasibound  state  in  the 
quantum  well  is  related  to  the  energy  width  of  the  corre¬ 
sponding  resonance  in  the  transmission  probability  by 
r  =  •  The  transmission  probability  is  calculated 

using  the  transfer  matris  approach  of  Kane,  modified  to 
account  for  the  different  effective  masses  of  the  particle  in 
the  quantum  well  and  in  the  barrier. 

For  electrons,  we  have  considered  only  F-point  barriers. 
It  is  appropriate  to  use  a  simple  one-band  expression  for  the 
wave  vector  in  the  well,  k  =  {2mZm,E /fr)  where is 
the  effective  mass  in  the  GaAs  well,  m,  is  the  free-electron 
mass,  and  E  is  the  energy  of  the  particle  with  respect  to  the 
GaAs  band  edge.  With  the  pure  AlAs  barriers  in  our  sam¬ 
ples,  the  lowest  quasibound  electron  state  has  an  energy  far 
from  the  band  edge  in  the.AJAs  barriers.  Hence  we  have  used 
a  two-band  model'*  to  calculate  the  wave  vector  in  the  bar¬ 
riers.  The  barrier  height  used  in  these  calculations  was  1.07 
eV,  corresponding  to  a  valence-band  offset  of  0.55  eV'*and 
an  AlAs  band  gap  of  3. 1 3  eV.  The  effective  masses  in  the  well 
and  the  barriers  were  taken  to  be  0.061  m,  and  0.15/n,,  re¬ 
spectively.  For  the  heaA^y  holes,  we  used  a  one-band  expres¬ 
sion  to  estimate  the  wave  vector  in  the  barriers.  From  the 
theoretical  curves  in  Fig.  4,  we  can  see  that  the  tunneling 
time  for  electrons  u  much  shorter  than  that  for  heavy  holes. 

Comparing  these  theoretical  estimates  of  the  tunneling 
times  with  the  decay  times  observed  experimentally,  we  note 
that  the  decay  times  agree  well  with  the  shorter  tunneling 
time  of  the  electrons.  In  contrast,  if  we  expected  Eq.  (1)  to 
explain  the  data,  then  we  should  observe  two  decay  times,  a 
short  one  near  zero  delay  and  a  longer  time  at  much  longer 
delays.  We  might  e.xpect  the  longer  decay  time  to  be  that  for 
the  heavy  holes.  There  are  a  number  of  phenomena  that 
could  contribute  to  the  observation  of  a  single  time  that  is 
close  to  the  electron  tunneling  time.  First,  with  the  widely 
differing  decay  rates  for  electrons  and  heavy  holes,  h  is  very 
likely  that  the  quantum  well  becomes  charged.  Rapidly  es¬ 
caping  electrons  could  leave  slowly  decaying  holes  behind, 
The  areal  densities  of  holes  and  electrons  produced  in  our 
experiment  are  roughly  10"  cm~^,  which  are  high  com¬ 
pared  to  the  background  doping  of  10*  cm“^.  These  areal 
densities  produced  by  the  photoexcitation  could  result-  in 
very  substantial  amounts  of  band  bending.  We  would  expect 
the  potential  across  the  harriers  due  to  charging  created  by 
the  differences  in  escape  rates  to  increase  the  tunneling  es¬ 
cape  rate  for  holes.  This  effect  would  tend  to  make  the  elec¬ 
tron  and  hole  densities  in  the  well  follow  each  other  and 
produce  a  decrease  in  the  tunneling  time  for  holes.  Second,  at 
higher  temperatures,  we  would  expect  significant  occupancy 
of  the  lowest  light  hole  subband.  The  tunneling  time  for  the 
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light  holes  is  even  shorter  than  that  for  electrons  and,  hence, 
could  offer  an  escape  channel  for  the  heavy  holes.  This  phe¬ 
nomenon  could  be  responsible  for  the  changes  in  the  ob¬ 
served  delay  scans  with  temperature.  Finally,  if  we  couple 
the  finite  spot  size  of  the  laser  with  the  variations  in  decay 
rates,  we  might  expect  density  gradients  and  electric  fields  to 
develop  in  the  plane  of  the  quantum  wells.  These  density 
gradients  and  electric  fields  could  act  to  make  the  hole  den¬ 
sity  follow  the  electron  density  by  inducing  hole  transport  in 
the  plane  of  the  quantum  well.  We  are  attempting  to  carry 
out  a  detailed  numerical  simulation  of  this  complex  problem 
which  will  be  important  in  understanding  the  results  of  any 
of  the  measurements  reported  to  date  that  are  based  on  opti¬ 
cal  excitation. 

In  summary,  the  tunneling  time  for  electrons  to  escape 
from  the  lowest  quasibound  state  in  the  quantum  wells  of 
GaAs/AlAs/GaAs/AlAs/GaAs  double-barrier  hetero¬ 
structures  has  been  measured  at  80  K  using  the  technique  of 
photoluminescence  excitation  correlation  spectroscopy.  Re¬ 
sults  are  in  good  agreement  with  calculations  of  electron 
tuimeling  time  based  on  the  energy  width  of  the  resonance  in 
the  transmission  coeScient. 
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We  have  theoretically  investigated  the  phenomenon  of  alloy  clustering  in  epitaxial  films  and,  in 
particular,  the  case  of  Al^,  Ga,  _  ,  As  grown  epitaxially  on  GaAs.  We  show  in  the  general  case  that 
the  elastic  strain  in  the  film,  caused  by  the  finite  mismatch  between  the  lattice  parameters  of  the 
film  and  the  substrate,  may  change  the  miscibility  properties  of  the  film.  A  miscibility  gap  could 
be  opened  in  the  material  composing  the  film,  evra  if  there  is  no  such  gap  in  the  bulk  material. 
Conversely,  the  strain  could  stabilize  a  solid  solution  in  the  film  in  cases  where  there  is  a 
miscibility  gap  in  the  bulk  material  of  which  the  film  is  composed.  We  estimate  that  these  effects 
are  too  small  to  account  for  phase  separation  in  the  case  of  Al,  Ga,  _  As  on  GaAs,  but  could  be 
important  in  other  systems,  such  as  In(As,  Sb). 


I.  INTRODUCTION 

Devices  consisting  of  epitaxial  layers  of  Al^^Gai.^As 
grown  on  GaAs  substrate  have  gained  prominence  in  elec¬ 
tronic  and  optoelectronic  research.'  In  many  cases,  the  per¬ 
formance  of  such  devices  depends  on  controlling  the  compo¬ 
sition  of  the  epitaxial  layer  (the  x  in  Al,  Ga,  _  ,  As).  In  recent 
years,  several  experimental  papers  have  been  published,  in¬ 
dicating  that  under  certain  growth  conditions,  the  composi¬ 
tion  of  the  epitaxial  layer  is  not  uniform.  Holonyak  et  al.^ 
observed  photoluminescence  features  that  were  compatible 
with  alloy  clustering  in  MOCVD  samples  of  Al^  Ga,  As 
grown  on  GaAs.  Petroff  et  al.^  have  observed  alternating 
layers  of  Al-  and  Ga-ricfa  Al;,Ga,  _,As  in  MBE  samples 
grown  on  GaAs(  1 10).  These  layers  were  parallel  to  the  inter¬ 
face,  about  100  A  tok,  and  the  variation  in  composition 
was  estimated  by  Petroff  ere/,  to  be  at  least  5%.  These  results 
are  somewhat  surprising  in  light  of  the  fact  that  bulk 
Al,  Ga,  _  ,  As  is  known  not  to  have  a  miscibility  gap.* 

Several  possible  mechanisms  could  lead  to  the  clustering 
in  the  epitaxial  layer.  First,  the  relatively  low  growth  tem¬ 
perature  suggests  that  the  observed  clustering  could  be  me¬ 
tastable;  that  would  not  explain  how  the  clustering  forms, 
only  why  it  is  observed.  Second,  the  observed  clustering 
could  be  stable  in  bulk  Al,  Ga,  _  ,  As  at  low  temperatures, 
and  the  thermodynamical  data  (obtained  from  high  tem¬ 
perature  experiments)  should  be  reevaluated  for  low  tem¬ 
peratures.  Third,  the  clustering  was  only  observed  in  thin 
films  grown  epitaitially  on  a  GaAs  substrate;  the  clustering 
may,  therefore,  be  induced  by  substrate  * . 

In  this  paper  we  examine  this  third  possibility,  namdy, 
that  the  substrate  induces  the  clustering.  We  investigate  &e 
effects  of  the  strain  in  the  film,  induced  by  the  substrate,  bn 
the  miscibility  properties  of  the  film.  (This  mechanism  was 
proposed  by  MaximmA  to  explain  clustering  effects  in 
GaAs,P,  _,.)  We  conclude  that  the  strain  term  conirk  in 
general,  change  these  miscibility  properties,  i.e.,  a  miscibility 
gap  could  be  either  opened  or  closdl  by  the  straiiL  To  be 
more  specific,  in  the  case  of  a  particular  system  like 
Al,Ga,  As  on  GaAs,  we  have  to  know  the  elastic  con¬ 
stants  of  the  film  over  the  entire  composition  range.  So  far. 


these  elastic  constants  have  not  been  measured  for  any  of  the 
ni-V  ternary  systems. 

This  paper  is  organized  as  follows.  In  Sec.  n,  following 
this  introduction,  we  discuss  the  effects  of  strain  on  phase 
separation  in  the  epitaxial  layer.  The  treatment  is  rather  gen¬ 
eral  in  that  section.  In  Sec.  m,  we  consider  some  m-V  sys¬ 
tems  in  which  the  strain  in  the  film  (induced  by  the  substrate) 
may  change  the  miscibility  properties  of  the  film.  In  Sec.  IV, 
we  attempt  to  estimate  the  relevant  parameters  for  the  parti¬ 
cular  case  of  Al,  Ga,  As  on  GaAs.  Our  estimates  of  the 
relevant  parameter  would  indicate  that  the  strain  term  is 
too  small  to  open  a  miscibility  gap  in  the  Al,  Ga,  _,  As  film, 
suggesting  that  the  observed  clustering  is  metastable.  In  Sec. 
y,  we  draw  conclusions  based  on  our  study. 


II.  EFFECTS  OF  STRAIN  ON  PHASE  SEPARATION 

In  this  section  we  estimate  the  excess  firee  energy  of  an 
epitaxial  film  due  to  the  induced  elastic  strain.  We  will  show 
t^t  this  excess  free  energy  may  have  both  positive  and  nega¬ 
tive  second  derivatives  with  respect  to  the  aUoy  composition 
(in  the  general  case  of  a  film  composed  of  a  binary  or  a  pseu¬ 
dobinary  alloy).  This  would  imply  that  the  strain  in  the  epi¬ 
taxial  layer  may  in  some  cases  stabilize  and  in  others  destabi¬ 
lize  clustering.  The  magnitude  of  this  excess  firee  energy  de¬ 
pends  on  the  elastic  constants  of  the  aOoy,  as  well  as  on  the 
lattice  mismatch.  If  the  lattice  mismatch  is  too  small,  the 
extra  stabilization  (or  destabilization)  of  the  dustering;  re¬ 
sulting  firom  the  dstic  strain,  could  be  negligible  compared 
to  the  enthalpy  and  entropy  of  mixing.  In  Sec.  IV,  we  will 
estimate  that  this  is  the  case  in  Al,Ga,  .^As  grown  on 
GaAs. 

When  there  is  amismatch  between  the  lattme parameters 
of  the  substrate  and  tile  film,  strain  effects  are  observed  in  the 
crystal  When  the  epitaxial  film  b  very  thin,  it  is  energetical¬ 
ly  favorable  for  the  film  to  be  strained  uniformly,  sneh  that  it 
adopts  the  lattice  parameters  of  the  substrate  parallel  to  the 
interface.  No  bonds  have  to  be  broken  in  this  case,  but  the 
elastic  strain  energy  accumulates  as  the  film  grows  thicker 
and  thicker.  At  some  critical  thirlrm^  misfit  dislocations 
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would  relieve  part  of  tliis  strain  (see,  for  example,  van  der 
Merwe  and  Ball*  and  references  therein).  In  our  treatment 
here,  we  should  only  treat  the  case  of  a  thin  film,  having  no 
misfit  dislocations. 

When  the  film  adopts  the  lattice  parameters  of  the  sub¬ 
strate,  its  elastic  strain  energy  is  given  by  Hook’s  law 

AE  =  \V.cr-,  (1) 

where  V  is  the  volume  of  the  film,  C  is  some  elastic  constant 
of  the  film,  and /is  the  strain  in  the  film.  The  strain  energy 
in  Eq.  { 1  j  could  depend  on  the  composition  of  the  alloy. 
For  discussing  phase  separation,  we  have  to  consider  only 
the  second  derivative  of  the  energy  with  respect  to  composi¬ 
tion.  Constant  and  linear  terms  in  composition  would  not 
alter  the  total  energy  of  the  system  upon  either  separation  or 
mixing.  If  the  second  derivative  of  the  strain  energy  with 
respect  to  x  is  positive,  the  strain  would  drive  the  system 
toward  mixing,  if  it  is  negative,  toward  phase  separation. 
The  strain  energy  is  only  part  of  the  total  free  energy,  which 
also  includes  enthalpy  of  formation  (the  linear  part  of  the 
enthalpy),  enthalpy  of  mixing  (the  nonlinear  part),  and  en¬ 
tropy  of  mixing  terms.  Consequently,  the  driving  toward 
phase  separation  or  mixing  due  to  strain  competes  with  the 
driving  toward  separation  or  mixing  due  to  the  other  terms. 
One  can  see,  however,  that  the  sign  of  the  curvature  of  the 
free  energy  could  sometimes  be  reversed  by  the  strain  term. 
For  example,  in  the  vicinity  of  coexistence  points,  where  the 
curvature  of  the  bulk  free  energy  is  zero,  it  is  the  elastic  strain 
term  that  could  make  the  difference. 

In  our  treatment,  we  assume  a  general  model  of  a  binary, 
or  a  pscudobinary  alloy,  say  Ai  We  further  assume 
that  this  alloy  is  grown  on  a  substrate  consisting  of  one  of  the 
two  materials  that  form  the  alloy  (say  A).  In  our  notation, 
therefore,  when  x  =  0  the  film  composition  is  identical  to 
that  of  the  substrate.  The  composition  dependence  of  the 
strain  energy  comes,  in  this  case,  from  the  composition  de¬ 
pendence  of  the  three  terms  in  Eq.  (1),  namely  V,  C,  and  f. 
Usually,  the  materials  in  question  are  nearly  lattice  matched, 
and,  therefore,  the  specific  volume  is  almost  independent  of 
composition.  We  shall  asstune,  therefore,  that  the  film  vol¬ 
ume  Fin  Eq.  (1)  is  independent  of  compositiciL  The  second 
term  in  Eq.  (1),  C,  can  be  expressed  in  terms  of  the  usual 
elastic  constants  of  the  film  and  depends  on  the  crystal  direc¬ 
tion  of  the  epitaxy.  In  the  case  of  a  cubic  film  grown  on  a 
cubic  substrate,  C  is  given  by^  -1  ■ 


2Cii  -h  2Ci2 


for  (100)/(100) 


12CJC.,  +  ^12)  fo,  (111, /(111). 

The  composition  dependence  of  C  can  be  calculated  from  the 
composition  dependence  of  C,„  C,2,  and  Cm.  Finally,  the 
composition  dependence  of  the  strain  is  linear,  according  to 
Vegard’s  law.  Therefore,  if  da  is  the  difference  between  the 
lattice  parameters  of  the  two  end  materials  which  form  the 
alloy,  then 


d£=F(^)'.yC(x)x^.  (3) 


In  Eq.  (3),  the  compositional  dependence  of  the  strain  energy 


is  lumped  in  the  term  C(x)x^.  It  should  be  emphasized  here 
that  this  term  may  have  a  second  derivative  with  respect  to  x 
of  both  positive  and  negative  sign.  To  demonstrate  this,  as¬ 
sume  that  C  (.x)  is  a  linear  function  of  the  composition 

C (x)  =  Co(l  —  x)  -t-  C,x  .  (4) 

In  that  case,  the  second  derivative  of  C(x)x^  is  also  a  linear 
function  of  x: 

[ C (x)x=]  =  2Q(  1  -  X)  +  (6C,  -  4Co)x .  .  (5) 

Therefore,  in  this  example,  C  (x)x^  always  curves  upward  for 
small  values  ofx  but  may  curve  downward  nearx  =  1  if  the 
slope  of  C  (x)  is  negative  enough  (C,  <  |  Cq). 


III.  TERNARY  lll-V  SYSTEMS 


As  we  have  mentioned  in  the  previous  section,  the  energy 
term  due  to  strain  in  the  film  has  to  compete  with  the  bulk 
terms  (enthalpy  and  entropy  of  mixing)  in  the  film  material. 
While  exact  numbers  are  not  available  in  most  cases,  we 
would  like  to  estimate  the  order  of  magnitude  of  each  of 
these  terms. 

The  entropy-of-mixing  term  dominates  the  free  energy  at 
high  temperatures.  In  many  of  the  III-V  ternary  systems. 


the  entropy  of  mixing  can  be  approximated  by  assuming  that 


the  solution  is  ideal: 


«  -  iVX:  [x logx  +  (1  - x)  log(l  - x)]  ,  (6) 


wheteN  'vs  the  number  of  formub  units.  Taking  the  second 
derivative  of  Eq.  (6)  with  respect  to  concentration 


C*  lx  1— xJ 


4kris  approximately  300  meV  at  600  'C.  This  temi  is,  there¬ 
fore,  quite  large  at  normal  growth  temperatures. 

We  would  like  to  comment  on  the  correction  of  Eq.  (6)  if 
the  solid  solution  is  not  ideal.  In  that  case,  different  arrange¬ 
ments  of  the  atoms  would  have  slightly  diffjiTnt 

urationsl  is  very  large,  we  may  assume  that  it  has  aCanssian 
energy  distribution,  with  average  n  and  standard  deviation 
<7,  which  depend  on  the  composition.  Now,  it  is  rather 
straightforward  to  use.  this  Gaussian  distribution  to  obtain 
both/r  and  S  as  functions  of  temperature: 


and  _  - '  - 

where  Su^.,  is  given  by  Eq.  (6).  It  should  be  noted  here  that 
the  expressions  (8)  and  (9)  ^ould  be  viewed  as  aqnnptotic 
expansions,  correct  only  in  the  limit  of  high  temperature. 
Expressions  such  as  Eqs.  (8)  and  (9),  as  well  as  hi^er  order 
terms,  were  given  by  Hild^rand  and  Scott.’  It  should  be 
noted  that  <r  vanishes  at  the  adpoints  of  the  composition 
range(x  =  Oorx  =  1)  and  attains  a  maximum  In  between.  At 
this  maximum  d^i^dx^cO.  Therefore,  the  conection  to 
the  free  energy  of  mixing  due  to  finite  temperatures  in  Eqs. 
(8)  and  (9)  has  a  positive  curvature,  thus  enhancing  mixing 
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Table  I.  Elastic  constants  dmes  the  volume  of  a  formula  unit  for  various 
ni-V  semiconductors  (in  units  of  eV).  C|,oo)  andC,iiii  were  defined  in  Eg. 
(2).  From  this  table  it  seems  that  a  value  of  100  eV  has  the  right  order  of 
magnitude  for  the  III-V  semiconductors. 


GaP  74.9  103.7 

InP  39.9  94.4 

GaAs  69.6  97.6 

InAs  33.2  83.5 

AlSb  60.6  87.0 

GaSb  62.0  86.3 

InSb  33.7  82-5 

rather  than  separation.  The  magnitude  of  these  correction 
terms  is  usually  smalL^ 

The  enthalpy  of  mixing  term  is  ±e  one  that  poses  the  main 
difficulty.  It  has  to  have  a  magnitude  of  more  than  a  few  teas 
of  meV  to  dominate  in  temperatures  of  a  typical  crystal 
growth.  The  enthalpy  of  mixing  in  most  ternary  (pseudobia- 
ary)  III-V  alloys  is  usually  approximated*  by  regular  solu¬ 
tion  theory: 

AIT' =  ax[\  -  X) ,  (10) 

where  a  is  an  interaction  parameter.  Experimental  estimates 
of  o  in  Eq.  ( 10)  for  several  III-V  ternary  systems  are  givm  by 
Panish  and  Ilegems.*  Various  phenomenological  models,  us¬ 
ing  several  unknown  parameters,  were  used  to  obtain  these 
values.  The  uncertainty  in  the  accuracy  of  these  parameters, 
therefore,  is  still  unknown.  Stringfellow®  has  observed  a  scal¬ 
ing  of  a  in  Eq.  (10)  like  the  square  of  the  mismatch  in  most 
III-V  alloys.  Martins  and  Zunger  showed  that  such  a  scal¬ 
ing  is  expected  on  theoretical  grounds  from  the  strain  com¬ 
ponent  of  AIT' .  They  have  further  showed  that  the  chemical 
effects,  due  to  charge  transfer  between  bonds,  do  not  scale 
with  Aa^,  but  for  random  alloys  this  contribution  is  small. 

We  shall  now  discuss  the  contribution  of  the  elastic  strain 
in  the  film.  The  excess  energy  due  to  elastic  strain  in  the  film, 
induced  by  the  substrate,  is  given  by  Eq.  (3).  To  calculate  this 
term,  we  have  to  know  the  elastic  constants  of  the  alloy  for 
the  whole  composition  range.  Such  a  complete  set  of  data  has 
not  been  measured  for  any  of  the  III-V  ternary  systems.  We 
can,  however,  estimate  the  magnitude  of  this  term. 

We  assume  that  the  second  derivative  of  J  C  (x)x^  with  re¬ 
spect  to  X  is  of  the  same  order  of  magnitude  as  C  (x).  This 
assumption  is  equivalent  to  assuming  that  there  are  no  . 
“spikes”  in  C  (x).  To  estimate  this  value,  we  present  in  Table  I 
the  relevant  elastic  constants,  C,,oo,  and  C„jt,  for  sey^ 
binary  III-V  semiconductors.  These  elastic  constants  ww  ' 
calculated  [using  Eq.  (2)]  and  later  multiplied  by  the  volume 
of  a  formula  unit  (ff’/4)  to  obtain  the  value  of  the  elastic  ' 
constants  per  formula  unit.  From  Table  I,  it  foQows  thar' 
C  pc)-a^/4  is  on  the  order  of  100  eV  per  formula  unit.  Multi- 
plying  this  value  by  the  square  of  the  strain,  we  obtain  that 
the  elastic  strain  contribution  to  the  free  energy  is  on  the 
order  of  0.1  meV  per  formula  unit  when  the  mismatch  is 
0.1%,  10  meV  per  formula  unit  when  it  is  1%,  and  up  to  • 

when  the  mismatch  is  10%.  ‘ 

Only  when  the  mismatch  between  the  fihn  and  the  sub-  .. 
strate  is  on  the  order  of  several  percent,  the  strain  term  is  , 


Self  di-fusion  of  Go  in  GaAs 


TEMPERATURE  (C) 

Fig.  1.  Characterisnc  sej^difiusion  times  of  Gs.  in  GaAs  according  to  the 
expeiunental  data  of  Goidstein  (Ref.  13,  solid  lines}  and  Palfrey  tt  al  (Ref. 
14,  dashed  lines).  Tlie  time  required  to  diffuse  over  characteristic 
of  5»  10,  and  100  A  is  shown  vs  temperature. 

appreciable  (on  the  order  of  10-100  meV).  It  may  still  be 
smaller  than  the  entropy-of-mixing  term,  but  when  the  en- 
tropy-of-mixing  term  is  nearly  balanced  by  the  enthalpy  of 
mixing,  the  strain  can  make  the  difference. 

For  a  precise  calculation,  one  has  to  know  both  the  elastic 
constants  of  the  film  material  as  a  function  of  composition 
and  the  thermodynamical  data.  At  present,  the  elastic  con¬ 
stants  of  the  alloys  have  not  been  measured,  artri  there  is  a  lot 
of  uncertainty  abom  the  thermodynamical  functions  of  these 
alloys  (see  our  discussion  of  Al,  Ga,  _  ,  As,  which  is  one  of 
the  best  characterized  alloys,  in  the  following  section).  We, 
therefore,  caimot  hope  at  this  moment  to  any  definite 
statements,  only  to  point  at  some  alloy  systems  that  could 
prove  to  be  interesting. 

We  are  interested  in  systems  that  have  a  large  mismatch 
( 1  %- 10%),  and  a  large,  positive  enthalpy  of  mixing.  Among 
the  ternary  HI— V  systems,  these  two  properties  are  correlat¬ 
ed’  since  the  chemical  contribution  to  the  enthalpy  of  mixing 
b  usually  small  for  random  alloys."  The  four  ternary  ni-V 
systems  that  include  InSb  are  of  particular  interest,  since 
they  all  satisfy  the  large  mismatch  condition.  In  adHitinn, 
>In(As,  Sb)  is  a  cantfidate  for  strained  layer  supeilattices  with 
small  band  gaps  that  may  be  used  as  IR  materials."  ' 

IV.  AI,Gax_;.  As  on  GaAs 

In  this  section  we  estimate  the  relevant  parameters  for  the 
Al;,Ga|_xAs  on  GaAs  interface.  Since  the  experimental 
data  is  incomplete,  we  have  to  estimate  several  parameters 
based  on  theoretical  considerations. 

A.  Kinetics 

■  Knowing  the  kinetics  of  the  system  could  help  clarifying 
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whether  the  clustering  in  Al^,  C&i  _ As  is  stable  or  metasta¬ 
ble.  Unfortunately,  only  partial  self-diffusion  data  are  avail¬ 
able  for  the  Al,  Ga,  _  ,  As  system,  and  they  are  mainly  high 
temperature  data.  In  Eg.  1,  we  show  characteristic  times  for 
self  diffusion  of  Ga  in  GaAs  vs  temperature.  There  are  two 
sets  of  curves  in  Fig.  1,  corresponding  to  diffusion  data  of 
Goldstein‘S  and  of  Palfrey  et  al.'*  In  each  set  there  are  three 


curves,  corresponding  to  diffusion  over  characteristic  dis¬ 
tances  of  5,  10,  and  100  A.  The  relation  between  diffusion 
distance  {d  ]  and  diffusion  time  [t )  was  calculated  according  to 

t  =  dyD,  (11) 

where  the  diffusion  coefScient  D  is  given  by 


\  —Q]  f?o=  iXlO’cm^s"' ,  2  =  5.6 eV,  (Goldstein'S) 

Z)  =  Z)oexp^-^J,  2)o  =  3.9xlO'ScmSs-*,  =  2.6 eV,  (Palfrey  etal.'*]. 


(The  discrepancy  between  the  two  sets  of  data  was  attributed 
by  Palfrey  et  al.  to  the  difference  in  As  overpressure  during 
growth.)  One  can  easily  see  in  Fig.  1,  that  nonequilibrium 
features  as  narrow  as  5  A  would  require  about  a  day  of  an¬ 
neal  at  800  C  to  dissipate,  if  we  taJre  Goldstein’s  data  as 
representing  the  actual  diffusion  rates.  This  would  imply 
that  the  features  as  large  as  those  observed  by  Petroff  et  al? 
do  not  have  enough  time  to  equilibrate.  On  the  other  hand,  if 
Palfrey’s  data  represent  the  actual  diffusion  rates,  then  at 
800  C  it  will  take  only  few  minutes  to  establish  quasiequili¬ 
brium  on  a  5  A  scale  and  several  hours  to  establish  quasiequi¬ 
librium  on  a  100  A  scale.  In  that  case  the  features  observ^ 
by  Petroff  er  al.  could  be  formed  spontaneously  from  a  uni¬ 
form  alloy  if  they  are  more  stable. 

The  large  discrepancy  between  the  two  sets  of  self-diffu¬ 
sion  data,  as  well  as  lack  of  self-diffusion  data  of 
Al^tGa,  _xAs,  does  not  enable  us  to  determine  whether 
thermal  equilibrium  was  established  based  on  annealing 
time.  We,  therefore,  have  to  determine  if  the  clustered  alloy 
is  more  stable  than  the  uniform  alloy,  based  on  energy  con¬ 
siderations. 

B.  Thermodynamics 

In  our  treatment  we  would  like  to  use  the  ideal-solution 
expression  for  the  entropy  of  mixing.  Recent  experimental” 
and  theoretical'  ‘  papers  have  questioned  the  ideality  of 
Al,  Ga,  _  ,  As  grown  at  the  600-800  “C  range.  We  may  have 
to  correct  the  entropy  of  mixing  to  account  for  nonideality, 
but  it  was  shown  in  the  previous  section  that  such  a  correc¬ 
tion,  at  least  to  first  order  in  1/T,  would  enhance  mixing 
rather  than  separation.  Assuming  that  the  Al,Ga,_,As 
solid  is  an  ideal  solution,”  the  entropy  of  mixing  was  gjven  in 
Eq.  (6).  The  entropy  contribution  to  the  second  derivative  of 
the  free  energy  is  at  least  ANkT,  which  is  300  meV  per  for¬ 
mula  unit  at  600  ’C  We  will  show  that  it  is  almost  an  ordff  of 
magnitude  larger  (in  absolute  value)  than  the  contiibutioia 
of  the  enthalpy  of  mixing,  hence,  dominating  the  mixmg 
properties  of  Al,  Ga,  As  above  600  *C. 

TTie  enthalpy  ofmixing  of  random  Al,Ga,  Ashasbea 
estimated  to  be  zero,  or  very  close  to  zero.  These  estimate^ 
were,  however,  based  on  high-temperature  measurements. 
Foster  et  a/.'*  have  investigated  the  solidus  line  of  the 

Al,  Ga,  _,  As  system  and  concluded  that  this  system  is  id^ 

[i.c.,  a  ==  0  in  Fxj.  (10)].  Stringfellow’  also  concluded  that 
a  =  0  in  the  Al,Ga,  _,As  system.  Other  figures  quoted  in 
the  literature  are  a——  500  cal/mol*  at  1000  *C, 
a  =  -f-  400  cal/mol'  at  700  *C,  and  or  =  —  3892  -1-  4r cal/ 
mol'  for  1073  K<r<1273  K.  (1  kcal/mol  =  43J  meV/ 


formula  unit).  The  experimental  figures,  therefore,  are  some¬ 
where  in  the  range 

[al  =  <25  meV/formula  unit,  (13) 

and  were  obtained  by  fitting  the  phase  diagram  of 
Al,  Ga,  _  ,  As.  To  obtain  yet  another  estimate  of  a,  we  have 
performed  ab  initio  self-consistent  calculations  on  As- 
(Ga,  A1)4H'2~  clusters.  Our  calculations  (to  be  described  be¬ 
low)  give  an  estimate  of  a  which  is  included  in  the  experi¬ 
mental  range  given  in  Eq.  (13)  in  the  separating  region 
[a  >  0),  but  extends  beyond  the  experimental  range  in  the 
mixing  region  (a  <0). 

Our  calculations  were  performed  as  follows.  Self-consis¬ 
tent  generalized-valence-bond'^  formalism  was  used  to  cal¬ 
culate  the  tot^  energy  of  AsGa4Hj2~,  AsGa2AliHf2“,  and 
AsAl4Hf2~  clusters.  The  quasichemical  interaction  was  cal¬ 
culated  according  to 

■Q  =  ■  (14) 

In  Eq.  ( 14),  .E  ajc^ai,  “  the  total  energy  of  the  AsGajAljHjj" 
cluster,  etc. 

.  All  the  clusters  in  our  calculations  consisted  of  a  central 
arsenic  atom  surrounded  by  a  shell  of  four  cations  (either 
galliuTn  or  aluminium  atoms),  which  was  in  tom  surrounded 
by  a  shell  of  12  hydrogen  atoms.  The  hydrogem  atoms  were 
used  to  terminate  the  dangling  bonds  of  the  cations.  The 
angles  between  the  bonds  formed  by  each  atom  and  its  near¬ 
est  neighbors  were  kept  tetrahedral  (109.47*)  throughout  the 
calculation.  The  cfamers  had  a  net  charge  of  3-  to  ensure 
covalent  bonding  smDar  to  that  in  the  crystal  The  alumi- 
num,  gallium,  and  arsenic  atoms  were  represented  by  a  dou¬ 
ble-^  basis  set,  containing  r-  and  p-like  wave  fhnctimis.  The 
hydrogen  atoms  were  represented  by  a  single^  s4ike  wave 
function  (a  minimal  basb  set).  The  distances  between  the 
: .  - '  hydrogen  atoms  and  their  neighbor  cation  atoms  were  kept 
at  a  prescribed  value  of  =  2.35  A,  rather  than  at  the  opd- 

;  inaldistance,andthehydrogenicwavefunctionswereinflat- 
v .  edbyafactorofrj,  =2.5.Thesetwoparametecs,>’f  andd„. 

>  were  used  to  control  die  electronegativity  of  the  hydrogenic 
A '  shell”  (that  was  used  to  mimic  the  bulk  Al,  Ga,  1,  As). 

The  value  of/2  niEq.  (14)  depends  on  the  iiufividn^  Ga¬ 
As  and  Al-As  bond  Iragths.  b  Fig.  2  we  show  the  total 
energy  of  the  AsCsjAliHu  cluster  (dashed  carves)  and /3  of 
Eq.  (14)  (continuous  curves)  as  contour  maps  in  the  plane  of 
;Ga-As  1^  Al-As  bond  lengths.  One  can  observe  in  Fig.  2 
that  a  value  of 

-0.020</2< -|-a004  eV  (15) 
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2.2  2.3  2.'t  2.5  2.6 


Fig.  2.  Total  energy  of  the  cluster  (dashed  lines},  and  the 

quasichemical  intemctiocs,  [li )  (solid  lines)  as  a  function  of  the  Ga-As  and 
Ai-As  bond  lengths.  The  interaction  parameter  i2  was  calculated  using  Eq. 
(14).  It  is  acting  toward  phase  separation  when  positive  and  toward  phase 
mixing  or  compound  formation  when  negative.  Its  value  is  rather  small  for 
all  reasonable  values  of  the  bond  lengths,  in  good  agreement  with  £({.  (13). 
Solid  triangles  represent  the  values  that  were  actually  calculated  (using  gen* 
eralized  valence  bond  medmd).  The  values  between  the  triangles  were  inter¬ 
polated. 


is  obtained  for  all  reasonable  values  of  the  Ga— As  and  AI-As 
distance.  The  /2  calculated  in  Eq.  (14)  include  all  the  bonds 
between  one  anion  and  its  neighboring  cations.  The  relation 
between  12  in  Eq.  ( 14)  and  a  in  Eq.  ( 10),  was  calculated  using 
the  quasichemical  model.''  We  assume  that  the  contribution 
to  the  excess  enthalpy  of  mixing  is  given  by  12  from  each 
arsenic  atom  that  is  bonded  to  two  gallium  and  two  alumi¬ 
num  atoms;  S2  /2  from  each  arsenic  atom  that  is  bonded  to 
one  cation  of  one  kind,  and  three  of  the  other  kind,  and  0 
'  from  each  arsenic  atom  that  is  bended  to  four  cations  of  the 
same  kind.  If  we  further  assume  that  all  the  configurations  of 
cations  are  equally  probable  [7'=  oo  in  Eq.  (8)],  then  the 
enthalpy  of  mixing  per  formula  unit  for  the  random  alloy  is 
given  by 

^  (o)  ^  (i)  f  (2) 

=  20.x(l-x)(l -Hx-x^).  (16) 

The  polynomial  ( 1 -i- X  —  x^)  varies  between  1  and  1.25  in  the 
range  0  <x  <  1.  Therefore,  by  comparing  Eqs.  (16)  to  (10),  it 
follows  that  a  ranges  between  2  and  2,5/2,  using  this  model.  . 
Therefore,  our  calculations  yield  a  value  of  a  in  the  range 
—  S0<a<  +  10  meV  per  formula  unit  This  estimate  is  in¬ 
cluded  in  the  experimental  estimate  given  in  Eq.  (13)  in  the 
separating  [a  >  0)  range,  but  extends  it  somewhat  in  the  mix¬ 
ing  (a  <0)  range. 

We  now  the  second  derivative  of  the  enthalpy  of  mix¬ 

ing  in  Eq.  ( 10)  with  respect  to  the  concentration  x  to  obtam  ai 
value  of  —  2a  per  formula  unit  which  is  more  than  —  50 
meV  per  formida  unh,  using  both  our  calculations  and  the 
experimental  data.  This  is  almost  an  order  of  magnitude 
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smaller  than  the  300  meV  contributed  by  the  entropy  of  mix¬ 
ing.  Om  calculadons,  therefore,  suggest  that  there  is  no  mis¬ 
cibility  gap  in  the  Al,  Ga,  _  j.  As. 

We  would  also  like  to  comment  here  about  the  formation 
of  an  ordered  compound  ( Al^,. j  Gao.3  As)  as  was  reported  re¬ 
cently  by  Kuan  et  al.  Such  a  compound  will  be  thermody¬ 
namically  stable  if  the  gain  in  enthalpy  of  mixing  more  rba  n 
offsets  the  loss  in  entropy  of  mixing.  According  to  the  quasi¬ 
chemical  model  described  above,  the  enthalpy  of  mixing  of 
such  an  ordered  compound  is  12  per  formula  unit,  compared 
to  1 12  of  the  random  alloy  of  the  same  composition.  Q 12  is 
the  value  of  Eq.  (16)  atx  =  J].  The  gain  in  enthalpy  of  mixing 
is,  therefore,  J 12  per  formula  unit,  according  to  this  model. 
The  entropy  of  mixing  of  the  compound  vanishes  due  to  the 
ordering,  while  the  entropy  of  mixing  of  the  alloy  equals 
—  k  log  |,  representing  a  free  energy  loss  of  about  0.7/cT, 
which  is  approximately  50  meV  per  formula  unit  at  600  *C, 
compared  to  enthalpy  gain  of  J 12  which  is  at  most  7.5  me'V 
according  to  our  calculations.  Srivastva  era/."  have  suggest¬ 
ed  that  charge  transfer  from  the  weaker  Ga— .As  bonds  to  the 
stronger  Al— As  bonds  may  result  in  a  negative  enthalpy  of 
mixing,  which  in  turn  results  in  compound  formation  at  low 
temperatures.  Our  calculation  of  12  described  above  agree 
with  the  sign  of  the  enthalpy  suggested  by  Srivastva  er  al., 
but  the  critical  temperature  required  for  such  compound  for¬ 
mation  is,  according  to  our  calculations,  much  lower  than 
600  ‘C  Our  calculations,  therefore,  suggest  that  the  ordered 
structure  observed  by  Kuan  et  al.  “  may  not  be  the  equilibri¬ 
um  state  of  Aljt  Ga,  _  As. 

Rnally,  the  elastic  constants  of  Al,  Ga,  _  ,  As  as  a  func¬ 
tion  of  X  have  not  been  measured  experimentally.  Fortunate¬ 
ly,  the  strain  in  an  AIxGa,  As  film  grown  on  GaAs  is  so 
small  that  we  can  easily  establish  the  fact  that  the  strain 
contributes  less  than  1  me'V/formula  unit  to  the  total  free 
energy.  The  difference  in  lattice  parameters  between  GaAs 
and  AlAs  is  extensively  documented  in  the  literature.  Etten- 
berg  and  Faff*’  showed  that  the  mismatch  varies  between 
0.14%  at  room  temperature,  to  0.0%  at  approximately 
950  *C  We  can,  therefore,  assume  that  above  600  *C, 

(^)’<I0-.  (171 

The  contribution  of  the  strain  for  such  a  snmn  misnmtrh  was 
estimated  in  the  previous  section  to  be  on  the  order  of  0.1 
me'V  per  formula  unit  Even  if  we  allow  an  order  of  magni¬ 
tude  error  in  the  estimate  of  the  second  derivative  of  C  (x)-x^, 
the  contribution  cf  the  strain  is  still  at  most  1  meV.  Su^  a 
small  contribution  cannot  affect  the  miscibility  properties, 
unless  the  entha^y  of  mixing  and  the  entropy-of-mizing 
terms  exactly  cancel  This  was  shown  not  to  be  the  case. 

V.  SUMMARY  AND  CONCLUSIONS 

In  this  paper  we  have  investigated  the  effect  of  substrate- 
induced  strain  on  the  miscibility  properties  of  an  qntaxial 
-film  consisting  (ff  an  alloy.  This  investigation  was  motivated 
by  the  observed  alloy  clustering  in  Al,  Ga,  _  ,  As  grown  epi¬ 
taxially  on  GaAs.  We  have  shown  that  in  some  cases  the 
miscibility  properties  of  the  film  can  be  reversed  by  the  strain 
induced  by  the  substrate.  A  miscibility  gap  could  be  opened 
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in  the  film  even  if  no  such  gap  exists  in  the  bulk  alloy,  and 
conversely,  a  miscibility  gap  in  the  bulk  alloy  could  be  closed 
by  the  strain. 

The  substrate-induced  strain  in  the  film  is  small  for  alloys 
that  are  closely  matched  to  the  substrate.  For  a  mismatch  of 
0.1%,  for  example,  we  have  estimated  the  strain  effect  to  be 
only  a  fraction  of  a  meV  per  atom  (in  the  case  of  III-V  ter¬ 
nary  alloys).  For  a  larger  mismatch,  say  1%,  the  contribu¬ 
tion  to  the  free  energy  is  on  the  order  of  10  meV/molecule. 
This  contribution  is  still  too  small  to  offset  the  entropy  of 
mixing  at  typical  epitaxial  growth  temperatures.  In  cases  in 
which  the  entropy  of  mixing  is  nearly  balanced  by  the  enth¬ 
alpy  of  mixing,  however,  the  strain  term  could  tip  the  ba¬ 
lance.  We  have  suggested  that  this  effect  may  be  important  in 
strained  layer  superlattices  based  on  In(As,  Sb).  Since 
strained-layer  superlartices  are  grown  at  relatively  low  tem¬ 
peratures,  miscibility  gaps  may  be  present  in  In(As,  Sb)  al¬ 
loys.  It  would  be  of  great  technological  interest  to  find 
whether  the  strain  can  lower  the  critical  temperature  of  the 
miscibility  gap,  or  alternatively,  enhance  directional  phase 
separation.  Currently,  lack  of  elastic  data  regarding  these 
alloys  as  well  as  uncertainty  about  the  thermodynamical 
data  prohibit  any  serious  attempt  to  calculate  these  effects. 

We  ha ve  also  investigated  the  Al,  Ga ,  _  As  on  GaAs  sys¬ 
tem.  In  that  case  we  have  concluded  that  the  strain  is  too 
small  to  affect  the  miscibility.  Since  clustering  was  observed 
in  this  system,  and  we  concluded  that  strain  effects  could  not 
have  caused  it,  we  have  examined  tne  thermodynamical  data 
in  more  detail.  We  have  performed  ab  initio  calculations  of 
the  enthalpy  of  mixing  Al;,Gaj  _,As  and  concluded  that 
d^H„/d:r  is  about  an  order  of  magnitude  smaller  tto 
<9  ‘TS„  /dx-  at  the  nonnal  growth  temperature  (r>  600  ‘Q. 
Our  calculations  are  supported  by  the  experimental  data  of 
Foster  et  which  indicate  that  the  Al^,  Gai  _  *  As  system 
behaves  as  a  nearly  ideal  soludon.  In  light  of  this  big  differ¬ 
ence  between  the  contributions  of  the  enthalpy  and  entropy 
of  mixing,  the  very  ^matl  additional  term  due  to  strain  can¬ 
not  open  a  miscibility  gap.  It  seems  likely,  therefore,  that  the 
uniform  random  alloy  is  the  most  stable  state  of  the 
Al,  Ga,  _  ,  As  system  at  temperatures  above  600  *C  We  at¬ 
tribute  the  fact  that  clustering  was  observed  in  that  system  to 


the  slow  kinetics  at  600-800  *C. 
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Growth  of  superconducting  VsSi  on  Si  by  molecular  beam  epitaxial 
techniques 
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California  Institute  of  Technology,  Pasadena,  California  91125 

(Received  26  April  1988,-  accepted  for  publication  6  June  1988) 

Superconducting  thin  films  (5(X)-1000  A  in  thickness)  of  ^  15  VjSi  have  been  grown  on 
( 1 1 1  )Si  substrates  through  molecular  beam  epitaxial  techniques.  V  and  Si  are  codeposited  in  a 
feedback-stabilized,  3:1  flux  ratio  from  dual  e-beam  evaporators  onto  clean,  heated  Si 
substrates  in  an  uitrahigh  vacuum  chamber.  For  the  first  dme,  the  superconducting  A  15  VjSi 
phase  is  metastably  fabricated  directly  on  a  crystalline  Si  substrate  without  the  formation  of 
thermodynamically  favored  Si-rich  Vj,Si,  phases.  Our  films  exhibit  superconductivity  only 
within  a  narrow  range  of  intermediate  growth  temperatures  (centered  near  4(X)  *C),.  with  an 
optimum  of  12.5  K.  X-ray  diffraction  shows  the  superconducting  films  to  be  polycrystalline 
A  15  VjSi.  For  slightly  higher  growth  temperatures  ( s530  *C),  Auger  profile,  x-ray,  and 
transmission  electron  microscopy  measurements  indicate  the  appearance  of  an  intermediate 
layer  of  roughly  1:1  V:Si  composition.  Our  observations  demonstrate  that  a  growth 
temperamre  near  400  *C  results  in  successful  nucleation  of  the  superconducting  .,415  phase 
while  minimizing  solid  phase  reaction  with  the  substrate  during  VjSi  growth. 


Materials  which  crystallize  in  the  15  structure  have 
long  been  of  interest  for  their  superconducting  properties. 
Recently,  one  of  these  .4  15  materials,  VjSi,  has  been  shown 
to  exhibit  a  close  geometrical  lattice  match  to  ( 1 1 1  )SL'  The 
lattice  match  between  VjSi  and  Si  is  shown  in  Fig.  1.  There 
we  see  that  the  ( 1 X 1)  surface  unit  cell  of  ( 1 1 1 )  VjSi  coin¬ 
cides  in  size  almost  exactly  (  <0.4%)  with  a  (V3X'\^)R30r 
surface  cell  of  ( 1 1 1  )SL-  The  nature  of  the  chemical  bonds  at 
such  an  interface  appears  favorable  to  epitaxy,  although  a 
significant  amount  of  reconstruction  may  be  necessary.  This 
lattice  match,  along  with  the  recent  advances  in  Si  molecular 
beam  epitaxy  (MBE)  which  have  made  possible  the  &brica- 
tion  of  epitaxial,  single-crystal  metal^  ( NiSij  and  CoSij)  and 
insulating^  (CaF2)  61““  on  Si,  invites  the  possibility  of 
growth  of  multilayered  device  structures  based  on  various 
combinations  of  epitaxial  superconductors,  semiconductors, 
normal  metals,  and  insolators,  on  Si. 

Unfortunately,  it  has  previously  been  impossible  to  sta¬ 
bilize  growth  of  VjSi  films  directly  on  crystalline  SL  Studies 
have  shown  that  the  thermodynamics^  and  ( in  some  cases) 
kinetics^  favor  the  formation  of  the  most  Si-iadi  va¬ 

nadium  silicide  phase,  VSi^,  on  a  crystalline  Si  substrate.  In 
particular,  the  V-Si  phase  diagram  contains  several  V^Si, 
compounds,  none  of  which  exhibits  superconductivity  ex¬ 
cept  VjSi,  the  most  V-rich  (and  in  the  presence  of  excess  Si, 
the  least  stable)  phase.*  This  previous  work  was  based  on 
solid  phase  reaction  (SPR)  of  thin  V  films  with  Si  rolv 
strates*'*  and  also  on  SPR  of  sequentiidly  deposited  Si-V  dif¬ 
fusion  couples  onto  Si02  or  AI2O3  substrates.*"^ 

In  this  letter,  we  report  the  first  successful  fabrication  of 
superconducting  VjSi  films  directly  on  a  crystalline  Si  sul^ 
strate.  Our  VjSi  films  were  fabricated  by  molecular  beam 
codeposition  of  V  and  Si  ( 3: 1  flux  ratio)  directly  onto  he^ed 
substrates  in  an  nltiahigh  vacuum  (UHV)  chamber. 
Growth  temperatures  significantly  lower  than  those  re- 

*’  Present  address:  Hugbes  Kesesrch  Labontorics,  Malibu.  CA  90263. 


quired  for  SPR  silicide  formation***  were  used,  thereby 
avoiding  the  problem  of  reaction  with  the  substrate  while 
retaining  sufficient  surface  mobilities  to  allow  formation  of 
the  .,4  15  structure:  For  an  appropriate  choice  of  growth  pa¬ 
rameters,  our  VjSi  films  consistently  exhibit  superconduc¬ 
tivity  with  7',,’sashighas  12.5  K  with  sharp  transitions  (Ar, 
sO.l  K).  While  all  the  films  studied  here  proved  to  be  poly¬ 
crystalline,  we  believe  that  this  work  takes  the  first  step  to¬ 
ward  achieving  an  epitaxial  superconductor  on  SL 

The  thin  films  studied  here  were  grown  in  a  custom- 
built  UHV  system  specifically  designed  for  silicide  growth. 
Silicon  and  a  transition  metal  (in  this  case,  V)  are  codeposit¬ 
ed  at  a  predetermined  flux  ratio  from  dual  e-beam  sources 
onto  heated  substrates!  A  tungsten  heater  filament  is  used  to 
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structures.  For  claiity.  only  OK  atoaic  hjer  of  the  Si  and  dine  of  die  VyS 
an  shown.  The  ftmdaiiiaitai  mtee  unit  eeO  ofVjSi  is  almoeteuetljr  the 
same  size  as  the  JCTsoiftoeoeO  on  ( 1 1 1  )SL  However,  snoeoaiy 

1/3  of  the  Si  atoms  fie  under  theSatoms  in  (comexsoftheibonAus), 

then  an  three  translatianally  eqnivaleat  lattice  matches  and.  for  eadi  of 
these,  a  pair  of  180*  ratadon  twins,  maldni  a  total  of  six  possMe  epfoudal 
matches  of  V,Si  to  a  pnn  ( 1 11  )Si  surface 
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radiatively  heat  the  substrate.  Substrate  temperatures  were 
calibrated  to  the  heater  power  level  through  observation  of 
Au/Si  and  Al/Si  eutectics,  giving  reasonably  accurate  tem¬ 
perature-power  level  correlation  in  the  range  363-5T7  *C 

Nominally,  50O-10CX)  A  films  consisting  of  a  3:1  ratio  of 
V:Si  were  codeposited  on  (lll)Si  and  (ll02)AljOj  sub¬ 
strates  at  a  rate  of  2.6  A/s.  The  fluxes  from  the  dual  e-bcam 
sources  were  closed-loop  controlled  using  two  quartz  crystal 
sensors  and  monitors  (Inficon  model  XTC)  working  in  par¬ 
allel.  A  3:1  V:Si  flux  ratio  was  calculated  to  result  from  a 
2.07: 1  V:Si  deposition  rate  ratio.  Our  growth  rate  for  VjSi  of 
2.6  A/s  was  calculated  assuming  that  the  V  and  Si  deposition 
rates  were  2. 1  and  1.0  A/s,  respectively.  Prior  to  deposition. 

Si  and  AUOj  substrates  were  degreased  in  trichloroethylene, 
acetone,  methanol,  and  de-ionized  water.  Si  substrates  were 
dipped  in  a  50%  by  volume  solution  of  HF  and  de-ionized 
water  in  order  to  remove  the  oxide  layer.  The  substrates  were 
outgassed  at  approximately  750  “C,  and  Si  substrates  were 
also  exposed  to  a  0.5  A/s  Si  flux  in  order  to  assist  in  the 
removal  of  the  oxide"'*  and  to  ensure  an  atomically  clean  Si 
starting  surface.  Depositions  were  made  in  a  low-hydrocar¬ 
bon  UHV  environment  with  a  base  pressure  of  8xl0~" 
Ton  and  with  typical  pressures  during  growth  of  approxi¬ 
mately  5x  lO”’  Ton  (dominated  by  Hj).  Following  each 
codeposition  sequence,  the  substrate  was  immediately  al¬ 
lowed  to  cool  in  order  to  minimize  the  amount  of  time  that 
the  films  were  held  at  elevated  temperature. 

Our  films  were  codeposited  on  Si  substrates  at  tempera¬ 
tures  between  room  temperature  and  approximately  6(X)*C. 
For  comparison,  similar  codepositions  wre  also  carried  out 
on  “inert”  AKOj  substrates.  The  films  grown  at  room  tem¬ 
perature  appear  minor  smooth,  yet  show  no  evidence  of  su¬ 
perconductivity  on  either  substrate.  Depositions  made  on 
AI1O3  at  higher  temperatures  also  appear  minor  smooth 
and  show  increasing  approaching  bulk  values  for  VjSi 
( =:  17  K).  Films  grown  on  Si  at  temperatures  in  excess  of 
550  *C  have  a  cloudy  appearance  and  show  no  evidence  of 
superconductivity.  Growths  on  Si  at  lower  temperatures 
have  surfaces  that  appear  minor  smooth  and  have  supemon- 
ducting  transition  temperatures  closer  to  10  K.  Figure  2 
shows  the^  vs  Tprofile  for  our  highest  sample  on  SL  This 
film  was  |rown  at  approximately  400  *C  to  a  final  thickness 
of  1000  A.  The  data  indicate  that  the  sample  undergoes  a 
superconducting  tiansidon  at  =  12.5  K  with  a  tiansidon 
width  of  about  0.2  K.  A  more  detailed  presentadon  of  the 
effect  of  the  growth  condidons  on  the  resultant  material  and 
electrical  properties  of  vanadium  silicide  films  on  Si  and 
AI2O3  will  be  given  elsewhere.'^  *~r,: 

Seemann-Bohlin  x-tay  diffracdon  measurements  on 
two  similar,  SOO-A-thick  superconducting  samples,  one 
grown  at  400  *C  and  the  other  grown  at  530  *C,  indicate  in 
both  cases  the  presence  ofapoiycrystalline  film  oD4  15,  ¥3^ 
This  is  confirmed  by  cioss-secdonal  transmission  electron 
microscopy  (TEM)  measurements  which  reveal  a  columnar  , 
polycrystalline  grain  strocture  with  a  grain  size  of  appioxi-  . 
mately  100-200  A.  Composidon  versus  depth  information 
for  these  two  samples  was  obtained  with  the  use  of  Anger 
sputter  depth  proffling.  Hgure  3(a)  shows  the  composition 
profile  for  the  sample  grown  at  400  *C.  As  can  be  seen  Grom 

^  CJ31 
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FIG.  2.  Resistivity  vs  temperature  data  for  a  superconducting  V,Si  film  on 
( 1 1 1  )Sl  The  film  shown  becomes  superconducting  at  7^  =  12J  K,  with  a 
transition  width  (10%-90%)ofapproxiinately  0.2  K.  Resistivity  was  mea¬ 
sured  using  the  van  der  Panw  method  and  calculated  assuming  a  film  thick¬ 
ness  of  1000  A. 

the  figure,  the  composidonal  profile  of  the  film  consists  of  a 
single  VjSi  layer  on  Si,  despite  a  somewhat  broad  interface. 
Assuming  the  VjSi  fihn  thickness  to  be  SOO  A,  we  estimate 
the  interface  width,  after  taking  instrumental  resoludon  into 
account,  to  be  approximately  140  A.  In  contrast.  Fig.  3(b) 
shows  the  compositional  profile  of  the  sample  grown  at 
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FIG.  3.  Auger  composifiBBri  profllcs  ibr  two  iQna,  eKh  aooiiosliy 
500  A  in  thickness.  Tbe  be  were  grown  at  (a)  400*CaDd  (b)  flppnuu* 
natdy  330  *C  Compowiiaaal  cafibntioo  is  adieved  by  defining  the  b  in 
(a)tobeaccnratdyV,Si.Thebioadintei&oeiteaiihprimnrilyfioniinter- 
bdaltonghiKssf— lOOAlmthesetnipies.*  *  “. -  - - 
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slightly  higher  temperature  (530  *C).  If  we  take  the  compo¬ 
sition  of  the  top  layer  to  be  VjSi,  as  suggested  by  our  x-ray 
and  electrical  measurements,  we  see  the  formation  of  an  in¬ 
termediate  layer  having  a  V:Si  composition  of  approximate¬ 
ly  1:1.  This  intermediate  layer  is  evidently  due  to  a  solid 
phase  reaction  between  the  depositing  VjSi  film  and  the  Si 
substrate.  It  is  interesting  to  note  that  no  “VSi”  compound  is 
present  in  the  V-Si  binary  phase  diagram,’  and  whiles  V^ij 
phase  does  occur,  this  phase  is  not  stable  below  1160*C. 
Hence  the  chemical  identity  of  the  intermediate  layer  is  un¬ 
clear  at  present,  though  it  conceivably  could  be  metastable 
V^Si,. 

The  intermediate  layer  is  also  clearly  evident  in  cross- 
sectional  TEM  images,  which  reveal  that  the  boundaries 
between  it  and  the  substrate  on  one  side  and  the  VjSi  over¬ 
layer  on  the  other  are  quite  rough  ( ~  1(X)  A).  This  interfa¬ 
cial  roughness  (rather  than  interdiffusion)  is  probably  the 
origin  of  the  broad  interfacial  transitions  observed  in  the 
Auger  profiles  shown  in  Fig.  3.  Thus,  from  the  above  results, 
we  see  that  a  growth  temperature  of  4(X)  *C  minimizes  the 
reaction  with  the  substrate,  yet  ensures  nucleation  of  the 
superconducting  A  15  VjSi  phase  on  Si. 

In  summary,  with  the  use  of  modem  molecular  beam 
epitaxial  growth  techniques  and  through  a  careful  choice  of 
growth  conditions,  we  have  successfully  fabricated  the  meta- 
stable  A  15  VjSi  phase  on  a  crystalline  Si  substrate.  Our 
study  demonstrates  that  the  1 5  V ^Si  structure  will  nucleate 
at  unusually  low  temperamres  during  molecular  beam  co¬ 
deposition  onto  a  clean,  heated  substrate.  Too  low  a  growth 
temperature  results  in  the  failure  to  form  the  superconduct¬ 
ing/1  15  structure,  while  growth  at  higher  temperatures  re¬ 
sults  in  a  reaction  between  the  depositing  film  and  the  Si 
substrate.  Although  we  have  not  yet  been  able  to  achieve 
epitaxy,  the  results  obtained  here  suggest  that  the  appropri¬ 
ate  growth  template  coupled  with  the  appropriate  substrate 
surface  preparation  may  result  in  epitaxial  growth  of 
(111 )  V,Si  on  ( 1 1 1  )Si  in  the  future. 
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Growth  and  Characterization  of  Superconducting  VaSi  on  Si 
and  AI2O3  by  Molecular  Beam  Epitaxial  Techniques 

E.  T.  CROKE,  R.  J.  HAUENSTEIN/  C.  W.  NIEH  and  T.  C.  McGILL 

California  Institute  of  Technology 
Pasadena,  California  91125 

A  study  of  the  growth  parameters  governing  the  nucleation  of  metastable  supercon¬ 
ducting  A15  VaSi  on  Si  and  AI2O3  is  presented.  Nominally,  500A  films  of  Vi_;,Si,  were 
produced  through  codeposition  of  V  and  Si  onto  heated  (111)  Si  and  (li02)  AI2O3  sub¬ 
strates.  Samples  were  prepared  in  a  custom-built  ultrahigh  vacaiun  (UHV)  ch^ber 
containing  dual  e-beam  evaporation  sources  and  a  high  temperature  substrate  heater. 

V  and  Si  fluxes  were  adjusted  to  result  in  the  desired  average  film  composition.  Vm^tSin.,^ 
films  prepared  at  temperatures  in  excess  of  550°  C  on  Si  show  significant  reaction  with 
the  substrate  and  are  nonsuperconducting  while  similar  films  grown  on  AI2O3  exhibit 
superconducting  transition  temperatures  (T^)  approaching  bulk  values  for  V3Si  (16.6— 

17.1  K).  Codeposition  at  temperatures  between  350  and  550“  C  results  in  supercon¬ 
ducting  films  on  Si  substrates  while  growth  at  lower  temperatures  results  in  nonsu- 
perconducting  films.  Lowering  the  growth  temperature  to  400°  C  has  been  shovra  through 
ex  situ  transmission  electron  microscopy  (TEM)  eind  Auger  compositional  profiling  to 
Tninimi^P  the  reaction  with  the  Si  substrate  while  still  activating  the  surface  migration 
processes  needed  to  nucleate  A 15  V3SL  Variation  of  film  composition  about  x  =  0.25  is 
shown  to  result  in  nonsuperconducting  films  for  high  x  and  superconducting  films  with 
Te  approaching  the  bulk  V  value  (5.4  K)  for  low  x.  Finally,  lowering  the  Vo.75Sioj25  de¬ 
position  rate  is  shown  to  raise  T^. 

Key  words:  Silicide,  superconductor,  VsSi,  codeposition,  epitaxy,  UHV,  MBE 


L  INTRODUCTION 

Recently,  much  excitement  has  been  generated  m 
the  scientific  community  due  to  the  discovery  of 
materiads  which  exhibit  high  superconducting  tran¬ 
sition  temperatures.  It  is  now  well  known  that  these 
materials  possess  a  hi^ily  layered  structure  resem¬ 
bling,  to  a  certain  extent,  a  superconductor/insu¬ 
lator  superlattice.  Our  work  with  the  V3Si/Si  syst^' 
has  been  motivated  by  the  possibility  of  fabricating 
such  a  structure  and,  perhaps,  gaining  control  over 
the  structural  considerations  governing  in  the  new 
superconductors.  In  addition,  superconductor/semi- 
conductor  epitaxial  structures  may  be  of  consider¬ 
able  interest  in  novel  electronic  device  applications. 

Besides  being  a  simple  binary  system,  VjSi/Si  also 
has  been  shown  to  have  a  very  close  lattice  match 
on  the  (111)  surface.  Successful  epitaxial  growth 
of  VsSi  on  Si  would  greatly  enhance  the  chances 
of  growing  layered  superconductor/semiconductor 
structures.  Zur  et  al?  have  shown  that  (111)  VsSi 
matches  to  (111)  Si  to  better  than  0.4%;  however, 
the  two-dimensional  unit  cell  size  on  the  (111)  Si 
surface  is  relatively  large  (three  times  the  area  of 
the  fundamental  unit  cell)  since  the  edges  of  the 
coninect  second  nearest  neighbor  atoms  as  opposed 
to  first  nearest  neighbors.  As  a  result,  the  question 
of  whether  or  not  VsSi  can  be  grown  epitaxially  on 
Si  is  important  from  both  a  physical  as  well  as  a 

.  technological  point  of  view.  Ai. 

*■  ,  .  _ 

*Pennanent  Address:  Hughes  Research  Laboratories,  Mahho, 
CA  90265.  -  r-vii 
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Unfortunately,  previous  work  on  solid  phase  re¬ 
action  (SPR)  of  V-5i  diffusion  couples®’*  and  SPR  of 
thin  V  films  on  Si  substrates®’®  has  demonstrated 
that  V3Si  is  unstable  on  crystalline  Si  as  samples 
have  a  tendency  to  react  fur^er,  resulting  in  growth 
of  the  nonsuperconducting  phase,  VSi2.  In  view  of 
these  results,  our  experiments  have  included  the  use 
of  AI2O3  substrates  in  addition  to  Si  substrates  since 
we  can  safely  assume  that  the  depositing  Vi^JSi,  film 
will  not  react  with  the  AI2O3. 

In  this  paper,  molecular  beam  epitaxial  (MBE) 
growth  techniques  were  used  to  produce  supercon¬ 
ducting  A15  VsSi  films  directly  on  a  crystalline  Si 
substrate.  Simultaneous  growth  on  Si  and  AI2O3 
substrates  at  temperatures  in  excess  of  550°  C  re¬ 
sult  in  nonsuperconducting  films  on  Si  while  the 
corresponding  films  on  AI2O3  exhibit  siq)erconduc- 
tivity  with  Tc  approaching  the  bulk  value  for  V3Si 
(16.6-17.1  K^).  Further  analysis  of  the  growth 
temperature  dependence  of  for  depositions  on 
Si  indicate  the  presence  of  a  plateau  region  of  su¬ 
perconductivity  between  350  and  *550“  C.' Auger 
compositional  profiling  and  transmission  electron 
microscopy  (TEM)  have  been  used  to  characterize 
the  quality  of  the  films  and  identify  the  scqwrcon- 
ducting  materiaL  Studies  of  the  dependence  of 
on  Si  concentration,  x,  reveal  that  a  maximum  oc¬ 
curs  near  z  =  0i!5.  b  addition,  for  hi^  x,  films 
prepared  on  Si  are  nonsuperconducting  and  fiir  low 
X,  Tc  approaches  the  bulk  value  for  pure  V  (5.4  K). 
Finally,  decreasing  the  growth  rate  is  shown  to  re¬ 
sult  in  improved  film  qu^ty  and  the  hi^iest  T.  ob¬ 
served  for  a  500A  V3Si  film  on  Si  (13.5  ED. 
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n.  EXPERIMENTAL 

The  samples  prepared  for  this  study  were  code¬ 
posited  from  a  Themionics  dual  e-beam  source  in 
a  custom-built  iiltrahigh  vacuum  (UHV)  chamber 
with  a  base  pressure  <8  x  10  Torr  and_ typical 
growth  pressures  of  approximately  5  x  10  *  Torr. 
The  growth  system  has  been  adequately  described 
elsewhere.®  Two  Inficon  (model  XTC)  controUers  were 
independently  used  to  maintain  closed-loop  control 
over  the  V  and  Si  deposition  rates.  Rates  were  cal¬ 
ibrated  using  a  Tencor  Instruments  alpha-step  200 
stylus  profiler.  Nominally,  lOOOA  V  or  Si  were  de¬ 
posited  on  AI2O3  substrates  and  the  thicknesses  were 
measured  at  four  different  spots  on  the  wafer.  The 
spread  in  the  measurements  indicated  an  uncer- . 
taiiity  of  about  5  percent  in  the  final  thicknesses. 
The  deposition  rate  ratio,  tJts,,  was  calculated  from 
the  expression 

£v  ^  I  -  X  PsjMv 

r^i  X  pviVfsi 

where  x  is  the  molar  fraction  in  Vi_;iSii,  psi  and  pv 
are  the  densities,  and  Msi  and  My  are  the  molecular 
weights  of  Si  and  V,  respectively.  For  x  =  0.25  (the 
Si  concentration  appropriate  for  VaSi  growth),  the 
ratio,  Ty/rsi  =  2.07,  was  calculated  from  Eq.  (1).  This 
ratio  changed,  of  course,  when  we  varied  the  Si  con¬ 
centration  in  the  depositing  film.  In  practice,  ac¬ 
curate  rate  control  over  a  few  angstroms  was  dif¬ 
ficult  to  achieve,  especially  at  low  deposition  rates 
since  instantaneous  rate  fluctuations  in  the  e-beam 
process  could  be  as  large  as  O.SA/sec;  however,  at 
the  growth  temperatures  used  in  this  study,  the  at¬ 
oms  have  sufficient  thermal  energy  to  blur  the  ef¬ 
fect  of  these  fluctuations.  As  a  result,  we  have  cho¬ 
sen  to  calculate  our  compositions  using  the  final 
thicknesses  of  the  constituents  in  Eq.  (1)  instead  of 
the  rates.  Propagating  the  uncertainty  in  the  thick¬ 
nesses  through  the  ciculation  results  in  an  uncer¬ 
tainty  in  average  Si  concentration  of  about  1%. 

Depositions  were  made  onto  heated  (111)  Si  and 
(ll02)  AI2O3  substrates  which  had  been  first  de¬ 
greased  in  trichloroethylene,  acetone,  methanol,  and 
then  rinsed  in  de-ionized  water.  After  this  degreas¬ 
ing  procedure,  the  Si  substrates  were  etched  in  a 
50%  HF  solution  to  remove  SiOj,  rinsed  once  again 
in  de-ionized  water,  and  immediately  loaded  into  the 
UHV  chamber.  Samples  prepared  for  the  study  of 
growth  temperature  versus  were  outgassed  01 
at  approximately  800"  C.  Growth  temperatures  were 
calibrated  through  observation  of  Au/Si  and  Al/Si 
eutectic  reactions  in  situ,  and  correlated  to  heater 
power  levels,  assuming  a  T*  dependence.  At  the  eu¬ 
tectic  points,  the  calibration  is  exact,  assuming  the 
eutectic  reactions  occurred  at  the  correct  tempera¬ 
tures  (363"  C  for  Au/Si  and  577’  C  for  Al/Si).  Sffl, 
we  estimate  our  substrate  temperature  calibration 
to  be  accurate  only  to  within  25°  C.  Other  factors 
are  also  present  which  may  affect  actual  substrate 
temperatxire.  For  instance,  thermal  drift  due  to 
heating  of  the  substrate  heater  itself  could  cause 


the  wEifer  temperature  to  change  during  deposition. 
We  have  attempted  to  address  these  factors  by  care¬ 
fully  reproducing  heater  power  levels  and  temper¬ 
ature  ramping  times  from  run  to  nm. 

Once  we  determined  that  we  could  consistently 
grow  superconducting  VsSi  on  Si,  we  discontinued 
use  of  the  AI2O3  substrates.  The  study  of  the  effect 
of  Si  concentration  on  used  a  modified  substrate 
cleaning  procedtire  and  employed  the  growth  of  a 
100 A  Si  buffer  layer  at  approximately  760"  C  in¬ 
stead  of  the  high  temperature  outgassing.  Prior  to 
growth  of  the  samples  used  in  the  study  of  growth 
rate  versus  Tg,  the  substrates  were  heated  to  ap¬ 
proximately  860°  C  and  were  exposed  to  a  O.lA/sec 
Si  flux  for  30  sec.*’^“  Subsequently,  a  200A  Si  buffer 
layer  was  grown  at  760"  C.  Changes  in  substrate 
cleaning  procedure  were  intended  to  improve  the  Si 
starting  siirface,  with  the  objective  of  obtaining  an 
atomically  clean  surface.  At  present,  we  do  not  have 
the  capability  of  characterizing  this  surface  in  situ 
aside  from  visual  inspection.  Fortunately,  subse¬ 
quent  TEM  was  able  to  show  that  most  of  the  Si02 
was  indeed  removed.  The  sample  set  described  in 
this  study  is  presented  in  Table  1. 


m.  RESLETS  AND  DISCUSSION 

Our  initial  attempts  at  growing  VjSi  involved  si¬ 
multaneous  deposition  onto  (111)  Si  and  (ll02)  AI2O3 
substrates  which  had  been  first  outgassed  in  situ  at 
approximately  800"  C.  V  and  Si  were  codeposited  in 

Table  I.  Vi-^Six/Si  sample  characteristics. 

Samples  1-9  were  used  in  the  study  of  TV  vs.  7^. 

The  dependence  of  Tg  on  composition  x  was 
determined  through  a  study  of  samples  5  and  10- 
13,  and  samples  14  and  15  were  used  to  investigate 
the  effect  of  deposition  rate  on  TV 


Thick- 


Sam¬ 

ple* 

Sub- 

strate(s) 

ness 

(A) 

x  (%Si) 

r,  r  C)  Tg(K} 

1(a) 

(111)  Si 

5000 

0.25 

25 

-  - 

1(6) 

(li02)AlA 

5000 

0.25 

25 

— 

2 

(111)  Si 

500 

0.25 

305 

3 

(111)  Si 

500 

0.25 

350 

8.4 

4 

(111)  Si 

500 

0.25 

380 

9.8 

5 

(111)  Si 

500 

0.25 

400 

9.8 

6(0) 

(111)  Si 

500 

0.25 

500 

9.6 

6(a) 

(li02)Al20i, 

500 

0.25 

500 

10.6 

7(a) 

(111)  Si 

500 

0.25 

530 

-.9.9 

7(6) 

(1102)  AlA 

500 

0.25 

530 

9.8 

8(a) 

(111)  Si 

500 

0.25 

565 

. ...  — 

8(6) 

(1102)  AlA 

500 

0.25 

565 

:.'l3.1 

9(a) 

(111)  Si 

525 

0.25 

640 

— — 

9(6) 

(1102)  AlA 

525 

0.25 

640 

■714.6 

10 

(111)  Si 

500 

0.20 

400 

5.4 

11 

(111)  Si 

600 

0.22 

400 

■■'7.0 

12 

(111)  Si 

500 

0.27 

400 

-  IIJ 

13 

(111)  Si 

500 

0.30 

400 

— 

14 

(111)  Si 

500 

0.25 

400 

.12.2 

15 

(111)  Si 

600 

0.25 

400 

13.5 

*Sample  pairs  denoted  by  (a)  and  (6)  involve  simultaneons  de¬ 
position  onto  Si  and  substrates. 
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a  3 : 1  flux  ratio  at  various  growth  temperatures  and 
Tc  was  measured  with  a  germanium-diode  ther¬ 
mometer  while  the  sample  was  held  over  liquid  He. 
Growth  rates  and  final  thicknesses  were  chosen  to 
result  in  the  growth  of  a  SOOA  V3Si  film  at  a  rate 
of  about  2.6A/sec.  Figure  1  shows  the  dependence 
of  Tc  on  Tg  observed  in  this  study.  Samples  pre¬ 
pared  at  temperatures  in  excess  of  550°  C  are  non¬ 
superconducting  on  Si  and  have  T*  approaching  the 
bulk  value  for  VsSi  (16.6-17.1  K)  on  AI2O3  sub¬ 
strates.  (Presumably  we  might  obtain  a  higher 
on  AI2O3  if  we  had  chosen  to  investigate  growth  at 
higher  temperatures.  Higher  growth  temperature 
would  give  the  depositing  atoms  more  energy  and 
freedom  to  nucleate  larger  grains.  (Xir  purpose  here, 
however,  is  to  study  growth  on  Si.)  Films  grown  on 
Si  do  not  superconduct  unless  the  growth  temper¬ 
ature  is  between  350  and  550°  C.  Room  tempera¬ 
ture  depositions  do  not  superconduct  on  either  sub¬ 
strate,  presumably  because  the  surface  migration  of 
the  depositing  atoms  is  insxifficient  to  form  the  A 15 
structirre.  Therefore,  we  find  a  plateau  region  be¬ 
tween  350  and  550°  C  where  A15  V3Si  can  be  grown 
directly  on  a  (111)  Si  substrate  with  Tc  near  10  K. 

Transmission  electron  diffraction  (TED)  and  plan 
view  TEM  have  confirmed  the  presence  of  poly¬ 
crystalline  A 15  V3Si  in  our  samples.  Figure  2(a)  is 
a  plan  view  TEM  of  a  500A  film  grown  at  400°  C 
(sample  5).  The  polycrystalline  nature  of  the  sur¬ 
face  layer  is  clearly  evident  in  the  image  and  the 


Growth  Temperature  (*’€) 

Fig.  1  —  Remits  from  T,  measurements  on  SOOA  VaTsSios  Hlms 
grown  on  (111)  Si  and  (1102)  Al»Oj  substrates.  Samples  grown 
on  Al^Os  at  hi^  tempentme  have  T,  approaching  the  bulk  value 
for  V]Si  (16.6-17.1  K).  Growths  on  Si  react  considerably  with 
the  substrate  at  growth  temperatures  in  excess  of  SaOT  C,  yet, 
between  350  and  550*  C,  superconducting  A15  VjSi  is  nucleated 
with  a  T,  typically  near  10  K.  Films  grown  at  room  temperature 
do  not  superconduct  on  either  substrate. 


(b) 

Fig,  2  —  a'  Plan  v.t  TEM  showing  randomly  or. sntsd  grains 
from  a  500.\  V  i.  grown  on  '111’  Si  at  400'  T  sample  5). 
The  grain  size  is  esn.-.ited  to  range  between  lOT  and  200A  in 
the  plane  of  the  intsrz;;,  'b)  Transmission  electrtn  dimaction 
pattern  along  the  111  tone  a.xis  normal  to  the  sunt'ace  of  the 
sample  studied  in  a  .-..ngs  surrounding  the  tra.-jmttted  beam 
icentrai  spot:  were  .ntrited  and  are  due  to  polycr  stalline  A15 
V;Si.  Rings  surrouni.-.g  diffraction  spots  due  to  :ne  substrate 
result  from  multiple  siittering  from  the  film  and  :ne  substrate. 


grain  diameter  i;  e.-timated  to  range  from  100-200A. 
Further  evidence  the  presence  of  polyg— .'smllme 
VjSi  is  presented  m  Fig.2lb).  This  image  is  the  elec¬ 
tron  diffraction  ci.::ern  from  the  samp'.e  analyzed 
for  Fig.  2iai  %v;:h  the  electron  beam  aligned  along 
the  I  lib  zone  ax;;  normal  to  the  surface.  The  rings 
surrounding  the  :ri.-.smitted  beam  icentrai  spot)  have 
been  inde.xed  an;  ere  due  to  VjSi.  Multiple  diffrac¬ 
tion  is  responsibi;  :br  the  rings  observed  around  the 
Si  substrate  dinra.nion  spots.  It  is  apparent  from  this 
pattern  that  the  Si  grains  are  random.iy  oriented 
with  respect  to  substrate. 

Two  oOO.A  S.  samples  grown  on  (111  Si  (sam¬ 
ples  .5  and  7'c  '  sre  selected  for  further  analysis 

using  .Auger  com; .  s.cional  profiling.  Calibration  was 
achieved  by  de:':r..r.g  the  surface  layer  to  be  accu¬ 
rately  \’iSi.  The  -rsults  are  presented  in  Fig.  3.  In 
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Fig.  3  —  Auger  compositional  profiles  &om  two  nominally  aOOA 
VsSi  samples  grown  on  (111)  Si  (samples  5  and  7(a)).  (aj  The 
sample  grown  at  400*  C  has  a  sharp  interface  with  the  Si  sub¬ 
strate  indicating  that  minimal  reaction  has  taken  place.  The  dif¬ 
fuse  nature  of  inter&ce  is  due  entirely  to  lii^tations  asso¬ 
ciated  with  instrumental  resolution,  (b)  A  growth  temperature 
of  530°  C  is  shown  to  result  in  significant  reaction  with  the  sub¬ 
strate  producing  an  intermediate  layer  with  a  1: 1  V:Si  average 
composition. 


Fig.  3(a),  a  film  grown  at  400°  C  is  shown  to  be  com¬ 
posed  entirely  of  VjSL  The  diffuse  nature  of  the  in¬ 
terface  with  the  Si  substrate  is  due  to  instrumental 
broadening.  Figure  3(h)  suggests  that  significant 
reaction  with  the  substrate  has  occurred  for  a  sam¬ 
ple  grown  at  530°  C  producing  an  intermediate  layer 
composed  of  V  and  Si  in  a  1 : 1  average  composition. 
It  is  interesting  to  note  that  the  binary  phase 
diagram*'  for  the  V-Si  system  does  not  indnde  a 
"VSi”  phase  and  that  although  a  VgSig  phase  does 
exist,  it  is  not  stable  below  1160°  C.  U^ortunately, 
a  plan  view  TEM  specimen  was  not  prepared  for  this 
sample  so  it  was  not  possible  to  identify  the  com¬ 
position  of  this  intermediate  layer  although  it  could 
conceivably  be  metastable  VgSij. 

Cross-sectional  TEM  was  also  performed  on  these 
samples  and  the  results  are  presented  in  Fig.  4.  In 
Fig.  4(a),  the  Si  buffer  layer  and  the  polycrystalline 
VaSi  overlayer  firom  the  sample  grown  at  400°  C  are 
clearly  visible.  The  buffer  layer  is  indeed  epitaxial 
despite  the  contrast  present  in  the  image.  Unfor¬ 
tunately,  it  appears  to  contain  a  higher  than  de¬ 
sired  density  of  extended  defects  (spaced  by  several 
thousand  A).  However,  since  the  VaSi  grain  size  is 


Fig.  4  —  Cross- St 
scribed  in  Fig.  I 
cr.'stalline  V.;Si  :■ 
the  sample  grc-vr 
530’  C  reaction  .v 
intermediate  la;  t: 
(b'  '.vith  typical 


(b) 

rZM  performed  on  the  tv.  :  samples  de- 
crystaliine  S:  buffer  lai-'e.-  and  a  poly- 
■T  IS  shown  to  :':rm  a  sha.—  interface  in 
0'  C.  ibi  For  a  growth  te.mperature  of 
Si  substrate  p.'oduced  a  r  ;  i;-‘crystalline 
as  appear  to  be  rolumnar  .r.  coth  (a)  and 
rs  ranging  from  100-'200.- 


oniy  about  lOO.A..  '.e  can  safely  assume  :'r.at  buffer 
layer  quality  is  the  limiting  factor  in  our  at¬ 
tempt  to  obtain  e:::axial  VjSi.  In  addi:i:n,  the  im¬ 
age  shows  r. :  e-.ience  of  reaction  tvith  the  sub¬ 
strate  since  the  ir.ierface  appears  very  sharp.  The 
530"  C  sample  sh;  -vn  in  Fig.  4  b'  has  o'c-nously  re¬ 
acted  with  the  S:  .-ubstrate.  One  can  clearly  distin¬ 
guish  two  poiycr  stalline  layers  in  this  image, 
thereby  corroboriiing  the  result  presented  in  Fig. 
3ibi.  The  crystalhies  appear  coium.nar  in  these  im¬ 
ages  with  tv-pical  grain  diameters  between  100  and 
200.4  in  agreement  with  our  plan  view  image  pre¬ 
sented  in  Fig.  2  a  Therefore,  we  conclude  that,  by 
lowering  the  grtv.nh  temperature  to  act  ut  400°  C, 
a  polycrystalline.  superconducting  VjSi  fiim  can  be 
grown  metastabl;  directly  on  crystalline  Si  with 
minimal  reaction  v-ith  the  substrate. 

In  addition  tc  siudying  the  dependente  of  T,  on 
growth  tempera: -.'e.  we  also  investigated  the  effect 
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of  film  composition  on  Depositions  were  made 
at  our  previously  determined  optimal  growth  tem¬ 
perature  of  400“  C  and  deposition  rates  were  mod¬ 
ified  to  provide  a  V:Si  flux  ratio  that  varied  from 
4:1  to  7:3  (20%  to  30%  Si).  Figure  5  shows  that,  for 
high  V  concentrations,  Tc  drops  to  about  5.4  K,  the 
accepted  value  for  pure  V,  suggesting  that  pure  V 
has  segregated  out  of  the  depositing  layer.  As  the 
Si  concentration  is  increased,  we  see  a  maximum  in 
the  Tg  occurring  near  i  =  0.25  as  might  be  expected. 
For  X  s  0.30,  samples  are  typically  non-supercon¬ 
ducting,  indicating  that  a  non-superconducting 
Vi_^4  phase  has  segregated  out  of  the  Si-rich  films 
in  the  same  way  that  pure  V  segregated  out  of  the 
V-rich  films.  So  far.  TEM  specimens  have  not  been 
prepared  in  order  to  verify  this  conclusion.  This  study 
has  shown  that  a  Si  composition  of  approximately 
X  =  0.25  results  in  the  highest  quality  supercon¬ 
ducting  VjSi  materiaL 

Several  500A  VaSi  samples  have  been  grown  at 
400“  C  in  an  attempt  to  study  the  effect  of  growth 
rate  on  but  unfortunately  have  yielded  wme- 
what  inconclusive  results.  Since  fluctuations  in  de¬ 
position  rate  become  more  significant  at  lower  rates, 
it  becomes  increasingly  difficult  to  maintain  an  ex¬ 
act  3:1  V:Si  stoichiometric  ratio.  Deposition  rates 
for  VaSi  growth  were  varied  between  2.6  and  0.8A/ 
sec.  The  first  few  growths  at  low  VaSi  deposition 
rates  have  produced  samples  with  higher  Te  values. 
Our  best  sample  shows  a  Tg  of  13.5  K  for  a  500A 
VjSi  film  grown  at  400“  C  with  a  deposition  rate  of 


only  about  O.SA/sec.  Although  this  result  is  en¬ 
couraging,  further  work  wrill  be  necessary  if  epitax¬ 
ial  VaSi  on  Si  is  to  be  realized. 


CONCLUSIONS 

Through  the  use  of  modem  molecular  beam  epi¬ 
taxial  growth  techniques,  superconducting  VjSi  films 
have  been  successfully  grown  directly  on  a  crystal¬ 
line  Si  substrate.  A  study  of  the  growth  parameters 
governing  nucleation  of  VaSi  on  (111)  Si  has  also 
been  conducted.  Growth  temperature  has  been  shown 
to  strongly  affect  the  nature  of  the  deposited  film. 
A  temperature  of  400“  C  seems  to  produce  the  high¬ 
est  qusJity  VaSi  while  minimizing  reaction  with  the 
Si  substrate.  Additionally,  for  Si  concentrations  dif¬ 
ferent  from  X  =  0.25,  Tg  begins  to  degrade  due  to 
segregation  of  inferior  materials  fiam  the  A 15  VaSi. 
Lowering  the  growth  rate  has  been  shown  to  in¬ 
crease  Tg  although  strict  stoichiometric  control  has 
been  difficult  to  obtain.  Finally,  though  our  results 
have  not  yet  indicated  epitaxial  growth,  the  use  of 
more  refined  growth  procedures  may  result  in  epi¬ 
taxial  (111)  VaSi  on  (111)  Si  in  the  near  future. 
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